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PREFACE 


ADDITIONS AND CORRECTIONS 


p. 28. Formula (30) is due to Lamd. Cf. A. E. H. Love, The Mathematical 
Theory of Elasticity^ 2nd edition, p. 55. 

p. 101. An asymptotic expression for (.s) wlien n is a large positive integer 
can be derived from a formula given by L. Fej^r in 1909. This 
formula is accessible in a paper by 0. Perron, Arhiv cler Mat. u. 
Phys. (1914). 

p. 118. The factor c in front of the double integrals should be omitted, 
p. 120. Pelete the minus sign in the second of equations (277). 
p. 127. Line 8. This statement is incorrect, the equations are poristic, the 
special case is the only one which can occur, 
p. 132. Line 20. On account of the porism just mentioned, the hope may be 
abandoned, 

p. 150. Ex. 13. For equations (10) of § 5 read equations (2) of § 2. 
p. 154. Ex. 24. The equation should read 


ir_£L 

dcLiJi-r 

+ 


oos (a - f) ^ ^ Sin (a - e) ^ J + |_{t - r) -g^J 

i [<r sin (a - e) I? + (t - r) g - <r (« - r) OOS (a - e) ^] = 0. 


at r em^ '' FaS ' 

computations. There are many parts of the subject indeed to 
which a pure mathematician might make useful additions , in 
' particular, I might direct attention to p. 21, line 2, and p. 101, 
where there are one or two matters which require further 
discussion. 




PREFACE 


r^jlHIS book is intended as an introduction to some recent 
developments of Maxwell’s electromagnetic theory which 
are directly connected with the solution of the partial differential 
equation of wave-motion. The higher developments of the 
theory which are based on the dynamical equations of motion are 
not considered at all. Even with this limitation the subject is a 
vast one, and to bring the work of perusing the literature within 
my power I have omitted an account of the modern theory of 
relativity which has been expounded very clearly in several recent 
publications. 

For a thorough understanding of the present subject a very 
extensive knowledge of mathematics is necessary, but there are 
parts of the subject in which a reader with only a limited 
mathematical equipment may soon feel at home and perhaps do 
useful original work. With the idea of enabling such a reader to 
obtain a quick grasp of the nature of the subject and the results 
obtained, I have thought it advisable to state without proof a 
number of relations of which adequate demonstrations can only 
be obtained by means of complicated and difficult analysis. 
I have also endeavoured to keep the analysis as elementary as 
possible, but in some places where the work is perfectly straight- 
forward a few details are omitted. 

The book is far from being a complete treatise on the subject, 
for I have not given any existence theorems to show that the 
solutions of certain problems exist and are unique, and no 
attempt has been made to enter into the details of numerical 
computations. There are many parts of the subject indeed to 
which a pure mathematician might make useful additions; in 
• particular, I might direct attention to p. 21, line 2, and p. 101, 
where there are one or two matters which require further 
discussion. 
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PKEFACE 


Chapter viii and paragraph 5 contain some of my own con- 
tributions to the subject. At present there se.em to be several 
different directions in which future developments may be made, 
and so it seems unwise to give a hasty judgment concerning the 
physical significance of the results. Ideas which naturally 
present themselves are that the aether can be regarded as built 
up from singular curves of the type considered in § 43, and that 
§§ 41 and 44 may throw some light on the question of the difference 
between positive and negative elementary electric charges. I 
hope to discuss an hypothesis relating to the first idea in a future 
note, but am unable to give any support at present to the 
second idea. 

I gratefully acknowledge my indebtedness to Sir Joseph 
Larmor who read the manuscript before it was revised and made 
some helpful suggestions, to Prof. Ames who read the greater 
portion of the manuscript, to Prof. Morley and Mr Ha8s6 who 
helped me with their advice and vigilance in reading the proof- 
sheets, and to the officers and staff of the University Press for 
their careful work and constant consideration shown in matters 
connected with the printing. For the correctness of the new 
formulae and examples I alone am responsible ; if any errors are 
discovered I shall be grateful if my readers will inform me. 

HARRY BATEMAN. 


October^ 1914. 
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CHAPTEK I 


FUNDAMENTAL IDEAS 

§ 1. The fundamental equations for free aether. 

In Maxwell’s electromagnetic theory the state of the aether 
in the vicinity of a point {x, y, z) at time t is specified by means 
of two vectors E and H which satisfy the circuital relations* 

^ IdE 1 , V 

rotA = 5- (1), 

c ot c ot 

and the solenoidal or sourceless conditions 
div E = 0, div H=0. 

If right-handed rectangular axes are used the symbol f rot H 
denotes the vector whose components are of type 

dy dz ’ 

the three components of H being Hx, Hy^ Hz respectively. 
The symbol div£f denotes the divergence of H, i.e. the 
quantity 

aF. 9^^ 

dx dy dz 

The vector E is called the electric displacement or electric 
force and H the magnetic force. The quantity c represents the 

* The equations are written in the symmetrical form in which they were 
presented by 0. Heaviside, Electrical Papers^ Vol. 1, § 30, and H. Hertz, 
Electric Waves, p. 138. Sir Joseph Larmor points out that a set of equations 
equivalent to these was first used by MacCullagh in 1838 as a scheme 
consistently covering the whole ground of Physical Optics, Collected Works of 
James MacCullagh (1880), p. 145. 

t We use here the units and notation employed in Lorentz’s The Theory of 
Electrons, Ch. i, except that large letters are used to denote vectors and E is 
written in place of D. Many writers use the symbol curl instead of rot. 
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velocity of propagation of homogeneous plane waves and is 
commonly called the velocity of light ; we shall assume it to be 
a constant, although in the most recent speculations it is treated 
as variable^. 

Some of the modern writers on the theory of relativity 
maintain that the introduction of the idea of an aether is 
unnecessary and misleading. Their criticisms are directed 
chiefly against the popular conception of the aether as a kind 
of fluid or elastic solid which can be regarded as practically 
stationary while material and electrified particles move through 
it. This idea has been very helpful as it presents us with 
a vivid picture of the processes which may be supposed to take 
place, it also has the advantage that with its aid we can attach 
a meaning to the term absolute motion, but herein lies its 
weakness. Larmor, Lorentz and Einstein have shown, in fact, 
that the differential equations of the electron theory admit of 
a group of transformations which can be interpreted to mean 
that there is no such thing as absolute motion. 

If this be admitted, the popular idea of the aether must be 
regarded as incorrect, and so if we wish to retain the idea of a 
continuous medium to explain action at a distance we must 
frankly acknowledge that the simplest description we can give 
of the properties of our medium is that embodied in the 
differential equations (1). 

If we abandon the idea of a continuous medium in the 
usual sense only two ways of explaining action at a distance 
readily suggest themselves. We may either think of the 
aether as a collection of tubes or filaments attached to the 
particles of matter as in the form of Faraday's theory which has 
been developed by Sir Joseph Thomson and N. R. Campbell ; 
or we may suppose that some particle or entity which belonged 
to an active body at time t belongs to the body acted upon at a 
later time t+r. From one point of view these two theories are 
the same, for if particles are continually emitted from an active 

* A. Einstein, Ann. d. Phys. Vol. 35 (1911), p. 898 ; Vol. 38 (1912), pp. 355 
and 443. M. Abraham, Phys. Zeitschr. (1912), pp. 1 — 5, 310 — 314, 793—797 ; 
Ann. d. Phys. (1912), pp. 444 and 1056; Fifth International Congress of Mathe- 
maticians, Proceedings^ Vol. 2, p. 256. 
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body they will form a kind of thread attached to it. The first 
form of the theory is, however, more general than the second. 

At present we are unable to form a satisfactory picture of 
the processes that give rise to, or are represented by, the vectors 
E and H, We believe, however, that some points may be made 
clear by studying the properties of solutions of our differential 
equations. 

It will be seen from the investigations of Chapter viii 
that the mathematical analysis connected with these equations 
is suitable for the discussion of three distinct theories of the 
universe, which may be described briefly as follows: — 


Aether 

Continuous medium. 

Discontinuous medium con- 
sisting of a collection of tubes 
or filaments. 

Continuous medium. 


Matter 

Aggregates of discrete par- 
ticles. 

An aggregate of discrete 
particles attached to the tubes. 

An aggregate of discrete 
particles to which tubes are 
attached. 


The last theory may be supposed to include that form of 
the emission theory of light in which small entities are projected 
from the particles of matter under certain circumstances and 
produce waves in the surrounding medium. This theory might 
be justly ascribed to Newton*. 

For other theories of the aether the reader is referred to 
Prof. E. T. Whittaker's recent workf A History of the Theories 
of the Aether, 

In the first part of this book the analysis is adapted almost 
entirely to the first theory, the high development of which we owe 
to the pioneer work of Maxwell, FitzGerald, Hertz, Rayleigh, 
Heaviside, J. J. Thomson, Lorentz and Larmor. The other 
theories have not yet received much attention but it is hoped 

* A form of the theory in which the entities are electric doublets has been 
developed by W. H. Bragg and applied to the X and y rays. British Associatiori 
Reports (1911), p. 340. 

t Dublin Univ. Press ; Longmans, Green and Co. (1910). 
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that the analysis of Chapter viii will lead to further develop- 
ments so that a comparison can be made between the different 
theories. It is quite likely that one theory will be enriched by 
the developments of another. 

§ 2. Electromagnetic fields. 

For many purposes it is convenient to work with a complex 
vector* M^H ± iE, where i = V — 1 and the ambiguous sign 
f is independent of the ambiguity which occurs in the determi- 
nation of V — 1. The differential equations (1) may then be 
replaced by the simpler equations 

rot if = + - div M = 0 (2). 

c dt ^ 

When a solution of these equations has been found a pair 
of vectors E and H satisfying equations (1) may be obtained by 
equating coefficients of the ambiguous sign. In working with 
an ambiguous sign it must be remembered that when two 
ambiguous signs are multiplied together the ambiguity is 
removed. The chief advantage in using the two independent 
ambiguities ± and - 1 is that we can assume that the vectors 
E and H are the real parts of expressions of the form and 
we are at liberty to equate the coefficients of either i or ± in 
any of our equations. 

Definition. A solution of the differential equations (2) or (1), 
which provides us with single-valued vector functions E and H 
for each space-time point (os, y, z, t) belonging to a certain 
domain D, is said to define an electromagnetic field in the 
domain D, 

Since the differential equations are linear the sum of any 
number of solutions is also a solution. The physical meaning 
of this is that when two electromagnetic fields are superposed, 
they are together equivalent to an electromagnetic field. 

Two superposed electromagnetic fields can of course be 
related to one another in some way. When electromagnetic 

* The use of a complex vector H- iJB is recommended by L. Silberstein, Ann, 
d. Pkys. Vols. 22 and 24 (1907) ; Phil. Mag. (6), Vol. 23 (1912), p. 790. He does 
not, however, use the ambiguous sign. 
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waves fall upon an obstacle, a secondary disturbance is produced 
which depends in character upon the nature of both the primary 
waves and the obstacle. 

We shall find that in some cases it is possible to find two 
fields in which the vectors {E, H), (E', H') are connected by the 
two relations embodied in the equation 

{MM') = + MyMy' + M,MJ = 0 (3) 

for all values of {x, y, z, t) belonging to some domain. 

When this is the case the fields are said to be conjugate 
within this domain. 

If we use the notation 

34^+14^ -fir/, 

we may write 

(iro = ~ {E^) ± 2i {EH) = 4 ± L, 

where 4 and 4 are two quantities which we shall call the 
invariants^. It is easy to see that when two conjugate fields 
are superposed the invariant 4 for the total field is the sum of 
the invariants Ji for the two component fields. Similarly for 
the invariant /a- 

When the invariants are zero over a given domain the field 
may be called self-conjugate for this region f. 

§ 3. The flow of energy. 

An entity whose volume density^' p is a function of {x, y, z, t) 
will vary in a manner which can be described as a simple flow 
with component velocities (w, v, w) if the equation of continuity 

w 

is satisfied. This equation implies in fact that there is no 

They are invariants for the group of linear transformations which leave 
the electromagnetic equations unaltered in form. Of. H. Minkowski, 

Nackr. (1908) ; E. Cunningham, Proc. London Math. Soc. (2), Vol. 8 (1910), 
p. 89 ; H. Poinoar4, Rend. Palermo (1906) ; M. Planck, Ann. d. Phys. Vol. 26 
(1908). Other invariants are given by these authors, 
t Silberstein calls it a pure electromagnetic wave. 

X The limitations to which the idea of density is subject and the question of 
the continuity of the function p are discussed by J. G. Leathern, Volume 
integrals and their use in physics,” Cambridge Mathematical Tracts (1905). 
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creation or annihilation of the entity in the neighbourhood of 
(x, y, z, t). 

Now it is easy to see that the equation of continuity is 
satisfied in virtue of equations (1) if we put 

p = + pu^c{EyH,--E,Hy\ 

and two similar equations. We shall regard p in this case as 
the volume density of the energy contained in the electro- 
magnetic field. The vector 2 whose components are of the 
type c (Eyllz — EzHy) can then be supposed to indicate the rate 
at which energy flows through the field. Since • 

(c^ — u- — — w^) = J {E^ — + c- {EHy, 

it appears that energy travels through the field with a velocity 
which is less than the velocity of light. The velocity c is 
attained only in the case of a self-conjugate field. 

The vector 2 was introduced by Prof. Poynting* and is 
usually called Poynting’s vector. The idea of describing the 
transfer of energy in this way also occurred to Prof. Lamb 
before the publication of Poynting’s work. 


Example, Prove that the equation of continuity may be satisfied V>y 
putting 

1 0(9 „ dQ ^ de 1 0^ „ ^ de 

d6 „ d8 „ de 


ay 


Ide „ dd „ ds „ 

pw=- E, - ^ 


= ~ S— JLJLx > pc = — ^ Jbx — 7^ M/y — -TZ-Ji 

0 ot ox ^ cy ^ ^ ox dy ^ dz 

where 6 is an arbitrary fmiction. Obtain a similar solution by replacing 
^by^andiS^by 




§ 4. First solution of the fundamental equations. 

Let us use the symbol Slu to denote the Dalembertianf 
of u, viz. 

^ . _ld^u dhi 1 d^u 

U ^2 0 jj 2 0^2 d " 0^2 0^2 ^2 05^2 ^ 

* Phil Tram, A, Vol. 175 (1884), p. 343. See also H. A. Lorentz, The 
Theory of Electrons, p. 22. 

t This is the name suggested by Lorentz, loc. cit, p. 17. Many writers use 
Cauchy’s symbol □ to denote the Dalembertian, hut I think Q is preferable 
because its form suggests a wave. Murphy’s symbol A is also used here in place 
of the usual symbol V^. E. B. Wilson and G-. N. Lewis use the symbol 0 
to denote the Dalembertian of u. Cf. Proc. Amer. Acad, of Arts and Sciences, 
7ol. 48 (1912), p. 389. 
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and the symbol grad U to denote the vector whose components 

du du dU 

are 

dx oy oz 

respectively. Let us also use HA, where A is a vector with 
components A^?, A^, A^, to denote the vector whose components 
are OA^;, HAjy, OA;^. The equation H 2 c = 0 will be called the 
wave-equation and a solution of this equation a wave-function. 
A vector function A will be said to satisfy the wave-equation 
when each of its components is a wave-function, i.e. if HA = 0. 
We may now satisfy equations (1) and (2) by writing 

Jf = ± rot X = ~ ^ ^ (^)j 

where the scalar pote7itial A=^ + i^ and the vector potential 
X = JS T iA satisfy the equations 

nA = 0, fli = 0, divi + i|^ = 0 (6). 

c ct 


The last three equations may be solved in a general way by 
writing 

L = - ^ ± 'irot (? + grad iT] 


c dt 
A = — div 

c dt 


.(7), 


where the vector G=T + ill and the scalar K satisfy the wave- 
equation 

Hu= 0 (8). 

The solution of equations (1) which is embodied in (5), 
(6) and (7) is a simple extension of Hertz’s solution* and is 
suggested by Whittaker’s solution*}* in terms of two scalar 
potentials. It is clear that the function K drops out when we 
differentiate to find M and so the electric and magnetic forces 
depend only on the vector G. The form of this vector indicates 
that the electromagnetic field can be regarded as the sum of 
two partial fields ; one of these is derived from the vector 11 and 


* Ann. d. Phys. Vol. 36 (18S8), p. 1. The general solution is given by 
Eighi, Bologna Mem. (5), t. 9 (1901), p. 1 ; IlNuovo Cimento (5), t. 2 (1901), p. 2. 
He finds suitable expressions for the vectors II and T in a number of cases, 
t Proc. London Math. Soc. Ser. 2, Yol. 1 (1903). 
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will be called a field of electric type, the other is derived fr 
the function T and will be called a field of magnetic type. 

This resolution of an electromagnetic field into two pari 
fields is analogous to the one used by H. M. Macdonald ^ in • 
study of the effect of an obstacle on a train of electric wa’\ 
The component fields are then of such a type that in one c, 
the magnetic force normal to the obstacle vanishes over ■ 
surface of the latter, in the other case it is the electric fo 
normal to the obstacle that vanishes. The same idea has b( 
used recently by Mief and Debye J in the treatment of the c 
of a spherical obstacle. 

In Hertz's solution we have T == 0, A" = 0 and 11 ’ 
components (0, 0, S). 

The components of E and H are consequently given by ' 
formulae 

^■^~^'c'dydt 



Hy = 


1 dhS [ 
c dxdt I 


( 9 ) 


„ _d^-s 1 a--S' 

df’ 




Hertz uses Eulers wave-function § 

A? = ” sin K{r — ct), = X“ + y- + 


and obtains in this way a theory of his oscillator||. The elec 
and magnetic forces become infinite at the origin which is the 
fore a singularity of the electromagnetic field. A singularity 
this type is called a vibrating electric doublet and is regarded 
the simplest model of a source of light or electromagnetic wa^ 

* EUctric Waves, Ch. vi. 

t Ann. d. Phys. Vol. 25 (1908), p. 382. t Ibid. Vol. 30 (1909), p. 6' 

§ Periodio solutions representing a disturbance sent out from n-fold p 
had been used previously by H. A. Rowland and applied to the elucidatio; 
optical phenomena. Amer. Journal of Mathematics, Vol. 6, p. 359 ; Phil. A 
Vol. 17 (1884), p. 423. Of. also Stokes, Camhr. Phil. Tram. (1849). 

II To deal with the ease in which the vibrations are damped we assi 

sin K (r- ci). Of. K. Pearson and A. Lee, Phil. Trans. A, Vol. 
(1900), p. 159. 
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The solutions of equations (1) which are obtained by superposing 
elementary solutions of this type are of great importance in 
physical optics. 

When r is very great the most important terms in the 
expressions (9) are 


E. = - 


HTXZS 





tc'yzs 




r- 






where 5 = sin /c (r ~ c^). All the other terms are of order I/?-- or 
These expressions give 

Hence at a very great distance from the origin the field is 
practically a self-conjugate field and so the energy travels with 
a velocity very nearly equal to the velocity of light. The 
expressions indicate that Poynting^s vector is ultimately along 
the radius from the origin; now the electric and magnetic 
forces are at right angles to Poynting s vector and so the 
vibrations of the light-vector, whether we take it to be the 
electric or magnetic force, are at right angles to the radius. 
The waves sent out from the source have, then, the character 
of monochromatic light at a great distance from the origin*. 
The amplitudes of the vibrations at points on the same radius 
are proportional to the quantities 1/r when r is large, and so if 
the intensity of the light be measured by the square of the 
amplitude the inverse square law is fulfilled. 

Since the electric force is ultimately at right angles to the 
radius there is no total charge associated with the singularity, 
for the charge is equal to the surface-integral of the normal 
electric force over a large sphere”" concentric with the origin 
and this integral is evidently zero. We are consequently 
justified in regarding the singularity as a doublet and in fact 


For a fuller diaoussion see Larmor, Phil. Mag. (5), Vol. 44 (1897), p. 503 ; 
Aether and Matter, Chap, xiv, where it is shown that energy is radiated from a 
moving charge only when the velocity of the charge alters in either magnitude 
or direction. 
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as a simple electric doublet of varying moment as is indicated 
by the way in which the electric and magnetic forces become 
infinite*. The axis of the doublet is along the axis of 

The electric lines of force due to a vibrating electric doublet 
have been drawn by Hertzf for various stages of the motion. 
The general character of the lines of force is indicated in Fig. 1. 



It will be noticed that the lines are all at right angles to a 
plane perpendicular to the axis of the doublet. M. Abraham § 
has used a Hertzian doublet to obtain a model of the electro- 
magnetic field produced by the oscillations in a vertical antenna, 
the plane just mentioned being supposed to represent the earth 
which is regarded as a perfect conductor. Zenneck|| has, 
however, pointed out that when the imperfect conductivity of 
the earth is taken into account the circumstances of the 
* See § 42. 

t Ann, d. Phys, Vol. 36 (1888), p. 1. The case of damped vibrations is 
considered by K. Pearson and A. Lee, loc. ciU 

t I am indebted to the Macmillan Company and A. Gray, Esg., for permis- 
sion to reproduce this diagram. 

§ Phys. Zeitschr. Vol. 2 (1901), p. 329; Theorie der Elehtrizitdt^ Vol. 2, 
§ 34 ; Encyklop, d. Math. Wiss. Band 5, § 18. 

II Ann. d. Phys. Vol. 23 (1907), p. 846 ; Phys. Zeitschr. Vol. 9 (1908), p. 50 ; 
lUd. p. 553. 
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propagation are somewhat different. The spreading of electro- 
magnetic waves over the earth’s surface has been investigated 
thoroughly by A. Sommerfeld* * * § * * * §§ and his pupil H. v. Hoerschle- 
mannf, and their results seem to indicate that the imperfect 
conductivity of the earth is an important factor in directing 
electric waves and in enabling their effects to be detected at 
great distances. The ionisation of the air by sunlight is also 
an important factor, as has been pointed out by J. J. Thomson, 
W. H. EcclesJ and J. A. Fleming^ Marconi’s experiments have 
indicated that the circumstances of propagation are not yet 
thoroughly understood. No good reason has been given to 
explain why communications by means of electric waves can be 
made more easily when the receiving station is in a north or 
south direction than when the direction is east or west. The 
curious contrasts in the results obtained with waves of different 
frequencies in day and night communications are also un- 
explained ||, 

The use of the vector n instead of the scalar S was 
recommended by Abraham IT. Von Hoerschlemann has obtained 
in this way a model of Marconi’s bent antenna which gives 
a directed effect to the radiation. A number of arrangements 
of Hertzian doublets that can be used to imitate the action 
of antennae have been described by Fleming^*, Larmor’|’*|-, 
Sommerfeld and Macdonald 1^. 

In the theory of FitzGerald’s magnetic oscillator§§ we have 

n = o, r = (o,o,iv^), 

N being Euler’s wave-function. Whittaker’s solution is 
obtained by adding the solutions of Hertz and FitzGerald. 

* A7m. d. Phys. Vol. 28 (1909), p. 665. 

t Jahrh. d. draht Teleg. Vol. 5 (1912). 

X Proc, Boy. Soc. A, Vol. 87, p. 79. 

§ British Associatio7i Beports, Dundee (1912). See also 0. J. Lodge, 
Phil. Mag. Vol. 25 (1913), p. 775. 

11 See Marconi’s address to the Royal Institution, June, 1913. 

H Theorie der Elektrizitat, Vol. 2, Ch. i. See also Bighi, loc. cit. 

Proc. Boy. Soc. A, Vol. 78, p. 1. 

ft Ibid, in a footnote to Fleming’s paper. 

XX Proc. Boy. Soc. A, Vol. 81, p. 394. 

§§ Trans. Boy. Enhlin Soc. Vol, 3 (1883) ; Scientific Writings, p. 122. 
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§ 5. Second solution of the fundamental equations. 

It is easy to see that equations (2) will be satisfied if we 
can find two functions (a, /3) such that 

djoL, i d (a, 13) 

^ d (y, z) - c a {oc., tj 

i d {a, y3) 

^ d{z, x) cd ( 2 /, t) 

M - ^ = 4-1 

d{xry)-- Cd{z,t) 

An electromagnetic field that is specified in this way is 
necessarily a self-conjugate field, for if we multiply together 
the two expressions for and do the same for My, Mz, we find 
that = 0. A particular pair of functions a, /S is obtained by 
putting 

« = ^ cos 0 -h y sin 0 == ii? sin 0 — 2 / cos O-’Ct (11), 

where 6 is an arbitrary constant. To generalise this field we 
multiply the expressions for My, Mz by an arbitrary 
function^ of a, ^8, ff and integrate with regard to 0 ; we thus 
obtain a very general electromagnetic field in which 

r2ir 

Msc —Ttj f{a, /3, 6) cos 6 dd | 

My^ + i ry(a, j3, ff) sin 0 del (12). 

J 0 

r27r 

= - f{oL,^,e)de j 

J 0 

The components of the electric and magnetic forces are 
obtained by equating the ambiguous and unambiguous parts in 
these equations; it is easy to verify that they are all wave- 
functions. 

It should be remarked that these definite integrals may 
give a representation of the electromagnetic field, required for 
the solution of a problem, only in a certain limited domain of 

* When -we speak of an arbitrary function it must be understood that the 
function may be subject to certain limitations which render the integration and 
differentiation under the integral sign intelligible operations. 
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the variables x, y, z,t\ the integrals may in fact represent 
discontinuous functions. 

The limits of integration could have been taken to be any 
other constants instead of 0 and 27r ; they can also be taken to 
be functions of x, y, z, t of the type a>, where o) is defined by an 
equation of the form 

X sin <» — ?/ cos — F {<*>), 

F being an arbitrary function. 

A suitable pair of functions a, ^ is also obtained by putting 

“ = i8 = r-ci (13). 


and in this case Poynting’s vector is along the radius from the 
origin. A more general type of electromagnetic field in which 
this is true is obtained by multiplying the above expressions 
for the components of M by an arbitrary function of a and /S. 

Other pairs of functions a, y8 of a very general nature are 
obtained in Chap. viii. It should be remarked that in all cases 
the functions (a, /3) are of such a nature that if F (a, /3) is an 
arbitrary function of a and F satisfies the partial differential 
equation 



which is of fundamental importance in geometrical optics* and 
may be called Hamilton's eqvMion, It is found that this 
equation is also satisfied in many cases by the functions of 
X, y, t which are the limits of a definite integral representing 
a wave-function, when the function under the integral sign is 
a wave-function for all values of the parameter with regard to 
which we are integrating. Thus the function (w just defined and 


the function t — (r + n) which will be used later are solutions 
c 


of this equation. 


* For another connection between this equation and the electromagnetic 
equations see A. Sommerfeld and J. Bunge, “Grundlagen der geometrischen 
Optik,” Ann, d. Phys, Vol. 35 (1911), p. 277. 
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I 6. The fundamental equations for a material medium. 
For a material medium which is stationary relative to the 
axes of coordinates, the equations (1) must be replaced by the 
more general equations* 


rot jE = — - ^ , div JT = 0 
c ot 


( 15 ), 


where D is the electric displacement, E the electric force or 
field strength, H the magnetic force and B the magnetic indue- 
dD 

tion. The quantity ^ represents the total current which 

is made up of a conduction-current C, a displacement-current 
dD 

^ and a convection-current pv, p being the volume density of 


electricity. 

Various notations have been used for the different vectors 
of an electromagnetic field. Most English writers use (a, 6, c) 
for the components of the magnetic induction, (a, y) for those 
of the magnetic force, (/, g, h) for the components of the electric 
displacement and (P, Q, B) or (Z, F, Z) for those of the 
electromotive intensity or electric force This is not to be 
confused with the mechanical force F of electromagnetic origin, 
whose components are sometimes denoted by (Z, F, Z). 


* Lorentz (1892—1895) and Larmor (1895) have derived these equations and 
a corresponding set of eq.Tiations for moving bodies by a process of averaging, 
starting from the fundamental equations of the theory of electrons in which we 
have D=E, J-pv. Cf. H. A. Lorentz, Afead. van Wetemchappen te 

Amsterdam (1902), p. 305 ; Encykl. d. Math. Wise. Bd. 5, § 14, pp. 200—210. This 
method of averaging has been developed so as to give results in accordance with 
the Theory of Belativity by M. Born, Math. Ann. Bd. 68 (1910) andE. Cunning- 
ham, Proc. London Math. Soc. Sex. 2, Vol. 10 (1911), p. 116. The Born- 
Minkowski equations differ slightly from those of Lorentz and indicate the 
existence of an electrostatic field due to the motion of a magnetised body. 

It has been realised by the foregoing writers and others that the principle of 
relativity alone is not sufficient to determine a complete set of equations for 
moving bodies, a theory of the constitution of matter is needed. Cf. H. B. Hass6, 
Phil. Mag. Jan. (1914). 

t Clerk Maxwell, Electricity and Magnetism^ 3rd edition (1892), Vol. 2, 
p. 267. 
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In any material medium there are certain constitutive 
relations connecting the vectors D, E, B, IT, E In moving 
media, crystalline media and ferromagnetic bodies the relations 
are rather complicated, but for an isotropic medium in which 
p = 0 the relations can be represented to a good degree of 
approximation by the simple equations 

D^eE, J^aE, (16), 

where e, cr, p. are scalar quantities which are generally regarded 
as constants; they will be regarded in fact as the optical 
constants of the medium. The quantity cr is called the con- 
ductivity^ jjb the pernieahility ^ and e the dielectric-inductive 
capacity. 

The units that are used here are the so-called modified 
units in which Heaviside’s suggestion of eliminating a factor 
4 < 7 r has been adopted. We can pass to electrostatic units or 
electromagnetic units by replacing our quantities J?, if, etc. by 
(xE, /3JE[, etc., where a, /3 are certain factors which are given in 
the following table: 



0, D, J 

E 

B 

E 

m 

Electrostatic system 

s / 4 it 

sj4:rr 

c 

^/47^ 

1 

C fJ^TT 

C ^/47^ 

Electromagnetic system 

C >/47r 

1 

C 

1 

/v/47r 

1 

\firiT 

tJ^TT 


We use e here to denote a quantity of electricity, and m a 
quantity of magnetism. 


§ 7. The energy equation for a material medium. 

If we use 2 as before to denote the vector whose components 
are of type o {EyHz - we find that 


' . as.v 


_L _L. ^ ^ W JT I TT t T^\ O 


* Of. H. A. Lorentz, EncyUopddie der Math. Wiss. Bd. 5, § 18, pp, 83 — 87. 
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If /i, is a constant and ^eE'^ be regarded as the 

energy per unit volume, the change in the distribution of the 
energy can be described by means of a flow a and a loss per 
unit volume of magnitude ctE'^ due to the transformation of 
electric energy into heat (Joule’s heat)^. If B does not depend 
on the instantaneous value of H so that is not a constant 
there is a loss of energy due to hysteresis. B may depend 
upon H alone but not be a single-valued function of //, con- 
sequently in a cycle of changes ^{H.dB) is not zero and may be 
taken as the heat per unit volume developed during the 
description of the cycle. Notice that 

\{H,dB) - - 5{B . dH) round a cycle, 
and is always positive since the value of B for a given value of 
H is greater when H is increasing. The experimental analysis 
and the accompanying theoiy are due to E. Warburgf and 
independently in much greater development to J. A. Ewing J by 
whom the name hysteresis was applied to such phenomena. 

§ 8. Solution of the fundamental equations for a 
material medium. 

Let us assume that cr, jjb, e are constants and that H 
the real parts of expressions of the form where A is 

a complex quantity independent of t Then if we write 

(IJ), 

and regard E, H now as the complex vectors of which they wcrc^ 
formerly the real parts, the differential equations to be satisfied 


by M are 

rot ibf == ± kMy div ilf = 0 (18). 

These may be solved by putting 

AT = rot n ± ~ grad div 11 ±kTl (19), 


where 11 is a solution of the equation Ati + and may b(j 

of the form U ±iV. 

* For a fuller discussion see R. Gans, MnfUhrung in die Theorie der 
Magnetismus; H. A. Lorentz, Encyhl. d. Math. Wise. Bd. 5, § 14, p. 240, Heft; 1 
(1903) ; Heaviside, Electrical Pajgers, Vol. 1, pp. 437 — 450. 

+ Ann. Fhys. Chem. (3), Yol. 13 (1881), p. 141. 

X Phil. Trans. A, Yol. 176 (1885), pp. 523 — 640; Proc. Roy. 8oc. Yol. 34 
(1883), p. 39; Phil. Mag. Yol. 16 (1883), p. 381; Magnetic Induction in Iron 
and other Metals, London (1892). 
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Problems which deal with the effect of small obstacles upon 
light or electric waves require the taking into account of the 
properties of the material of which the obstacle is composed : it 
is only for long waves and conducting media that the obstacle 
can be treated as a perfect reflector. This is illustrated by the 
work of Maxwell Garnett^ and Mief on the optical properties 
of colloidal suspensions of metals and in the work of Sommerfeld 
to which we have already referred. 

Unfortunately, however, the analytical difficulties are very 
great when imperfect conductivity is taken into account The 
simple-looking problem of the reflection of the disturbance 
produced by a moving charge, when the obstacle is an infinite 
plane sheet of metal or other conducting substance, has not yet 
been solved accurately § and there are many similar problems 
that have completely baffled mathematicians. 

Much more progress has been made with problems dealing 
with perfect reflectors. These problems are to some extent 
ideal but some of the characteristics of actual physical problems 
are often preserved ||. Apart from this, such problems are of 
considerable mathematical interest, and have been studied by 
some writers simply on this account. 

§ 9. Boundary-Conditions. 

The conditions to be satisfied at a surface separating two 
different media are obtained by integrating equations (15) across 
a thin layer of transition 11. Taking the axis of z along the 

* PM. Tram. A, Vol. 203 (1904), p. 385 ; Vol. 205 (1905), p. 237. 

t Ann. d. Phys. VoL 25 (1908), p. 377. 

t An important solution of the equations has been given by M. Brillouin, 
“Propagation dans les milieux conducteurs,” Comptes Eendus, t. 136 (1903), 
pp. 667, 746. (See Ex. 20, Ch. ix.) For other references see Ex. 5, p. 23. 

§ An approximate solution was suggested by Maxwell and has been developed 
by Larmor, Proc. London Math. Soc., Ser. 2, Vol. 8, p. 1. An accurate solution 
for the case in which the sheet is treated as infinitely thin and the charge moves 
with uniform velocity parallel to the plane has been given by G. Picciati, Bom. 
Acc. Line. Bend. (5), II 2 (1902), p. 221. 

II A metal behaves practically as a perfect conductor to electric waves when 
the oscillations are rapid but slow compared with waves of light. Cf. J. Larmor, 
“Electric vibrations in condensing systems,” Proc. London Math. Soc. Ser. 1, 
Vol. 26, p. 119. 

IT Eayleigh, Scientific Papers, Vol. 1; Phil. Mag. Vol. 12 (1881), p. 81; 
H. Hertz, Electric Waves, pp. 207, 238; Larmor, Phil. Tram. (1895), 


B. 


2 
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normal to the surface we shall assume that E, H, D, B a; 
their derivatives with regard to sc, y, t are finite within t 
layer and that the conductivity or is also finite. The equatio 

dz dy 0 dt ^ dy dz g\^ dt J ’ 

then show that ^ and /a? + ^ are finite and so th( 

integrals with regard to z across a thin layer of thickness 6 a 
less than ad where a is a finite positive quantity independent of 
This means that the tangential components of the electric for( 
magnetic force and electric current are continuous in crossii 
the surface. 

Again, the equation 


dBg^ ^ dBy , dBz _^ 

dx dy dz 


dB, 


shows that — is finite and so the normal component of tl 

magnetic induction is continuous. This result may, howeve 
be regarded as a consequence of the previous one. 

The normal component of the electric displacement may 1 
discontinuous, for the equation div D = p gives 


/ ^dD 

dz==J pdz + terms of order 6 , 

where dz is the discontinuity in the electric displacement. ' Hen 
if o is the surface-charge of electricity per unit area we have dz= 
When the media are both conducting we have 5= = 0 and tl 
normal component of the electric displacement is continuous. 

It is known that in the case of a good conducting boc 
rapidly alternating currents are confined within a very th; 
layer close to the surface*. In the ideal case of a perfe 
conductor or perfect reflector the field-vectors are zero with; 


p. 733; H. M. Macdonald, Electric Waves, p. 14; H. A. Lorentz, Zeitschr 
f. Math, u, Phys. Bd. 22 (1877). 

* Of. H. Lamb, Phil. Trans. A (1883) ; 0. Heaviside, Electrical Papers, VoL 
p. 168 ; J. J. Thomson, Becent Besearches, p. 281. For some recent work ' 
the subject see a paper by E. F. Northrup and J. E. Carson, Jov/rn. of i 
Franklin Institute, Feb. (1914), p. 125. The results of many other research 
on the skin-effect and alternating current resistance are given in J. A. Fleming 
The principles of Electric Wave Telegraphy. 
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the body of the conductor. This case is characterised by or = oo 
and the tangential components of the magnetic force are no 
longer continuous as a point moves across the surface, we have 
in fact if h is the discontinuity in the magnetic force 

h = / — crEsodz 4 * terms of order 0. 

^ Jo dz cJo 

Hence if K is the surface-current and h the discontinuity 
of the magnetic force, we have 

hy=-lK^, K=-Ky. 
c c 

At the surface of a perfect conductor the tangential com- 
ponents of the electric force must vanish and as a consequence 
of this we can say that the normal component of the magnetic 
induction must also vanish. If the medium outside the con- 
ductor is free aether the surface-conditions are simplified on 
account of the relations D = E, 

In the case of a very thin conducting sheet it is convenient 
to treat the thickness of the sheet as negligible and regard the 
tangential components of the magnetic force as discontinuous 
when a point moves along the normal from one side of the 
sheet to the other. 

The boundary-condition is then* 

-hy = ^E^, h^ = ^Ey, 

C 0 

where adz is the conductivity of the sheet. 

If we wish to extend the idea of Green's equivalent layer to 
electrod3m.amics we must consider electromagnetic fields in free 
aether with surfaces at which the tangential components of the 
electric and magnetic forces are discontinuous ; this requires an 
electric current sheet and a magnetic current sheet on the 

* Of. T. Levi-CivM, Rend. Lmeei (5), llg (1902), p. 75. These conditions 
are used by Picciati in his solution of the problem of an electric charge moving 
parallel to a conducting sheet. In some previous papers, Rend. Lmeei (5), lli 
(1902), pp. 163, 191, 228, Levi-Civit^ had used the electromagnetic potentials to 
determine the effect of a conducting sheet on an alternating current flowing 
along a straight wire parallel to the sheet ; the boundary-conditions are then 
determined by the discontinuities of the potentials due to the induced current. 
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surface*. For a complete generalisation we ought to consider 
the cases when the normal components are also discontinuous 
and when the surface is in motion. The last circumstance 
alters matters to some extent and must now be discussed. 

In t"he case of electric waves in free aether, the vectors 
E and H may be discontinuous at a wave-boundary. If this 
can be regarded as the limit of a thin layer of transition within 
which equations (1) are satisfied and the vectors E, H are 
finitef, the values of these vectors on the two sides of the 
boundary must satisfy certain conditions which may be found 
as follows. 

Let the equation of the moving boundary be expressed in 
the form • 

( 20 ). 

If now we apply Green's theorem to the integral 

I [[(if, + icM, I ± tcM, |) dyd^ 

+ (ify + icMg ~ + fcif, dzdx 

+ {m, + icM^ ^ + ioMy dxdy^ ( 21 ), 

which is supposed to be taken over a closed surface, we find 
that it vanishes on account of the equations (2) provided t is 
supposed to be expressed in terms of x, y, z according to some 
definite law which we shall take to be that expressed by (20). 

We now apply this theorem to a disc-shaped surface whose 
two faces very nearly coincide. We shall suppose that on one 
side of the disc the vector M represents the field of the advancing 
waves and that on the other side it represents the field obtain- 
ing just before the arrival of the waves. We shall also suppose 

* Of. J. Larmor, Froc, London Math. Soc. (2), Yol. 1 (1903), p. 11; 
H. M. Macdonald, Electric Waves, p. 16 ; Froc. London Math. Soc. (2), Yol. 10 
(1911), p. 91. 

t The idea is practically due to Stokes, Math, and Fhys. Fapers, Yol. 2, 
p. 275, but was not worked out in detail. The different possible types of dis- 
continuity are discussed with some care by Love. The case in which E and H 
are continuous but some of their derivatives are discontinuous at the moving 
boundary may be discussed more simply by analysis analogous to that given fn 
Hadamard’s Lemons sur la Fropagation des Ondes, Paris (1903), Ch. 2. See also 
Ex. 2, p. 23, and the references to Duhem and Silberstein on the next page. 
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that the derivatives of M are finite or behave in such a way 
that an application of Green’s theorem is justifiable. Now 
let M be the discontinuity of M, i.e. the difference in the values 
of M at two neighbouring points on opposite sides of the wave- 
boundary. Then when the two faces of the disc coincide we 
find that a certain surface-integral over one face of the disc is 
zero. The surface-integral is of the same type as (21) except 
that M is written in place of M, Since the face of the disc can 
be chosen arbitrarily the integrand must vanish and so we 
obtain three equations of the type* 

#* + ic J’j, £ ± ic J, ^ = 0 (22). 

These equations give 





1 


Hence the wave-front advances with the velocity of light 
and the difference between the two electromagnetic fields at 
the wave-boundary behaves as a self-conjugate field in which 
Pojmting’s vector is along the normal to the moving 
boundary. 

If the equation of the wave-boundary be expressed in the 
form 

F{x, = 


we find on calculating the values of 


dx' dy^ ' dz' 


that 




( 14 ). 


This is the differential equation of the characteristics, it 
expresses that the moving boundary moves normally to itself 
with the velocity of light. According to the theory of Stokes f 


* Equations equivalent to these are obtained in a different manner by 
0. Heaviside, Electrical Papers^ Vol. 2, p. 405 ; A. E. H. Love, Proc. London 
Math. Soc. Ser. 2, Vol. 1 (1903), p. 37 ; L. Silberstein, A7in. d. Physik, Vol. 26 
• (1908), p. 751 ; P. Duhem, Comptes Bendm^ t. 131 (1900), p. 1171. 

t Proc. Camb. Phil. Soc. (1896) ; Manchester Memoirs (1897) ; Math, and 
Phys, Papers^ Vol. 4. 
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and Wiechert*, Eontgenrays consist of pulses travelling through 
the aether, the energy in a pulse being confined within a thin 
shell. The above theory indicates that the front and rear 
surfaces of the shell move forward with the velocity of light. 

A slight modification of the preceding method can be used 
to find the conditions to be satisfied at a moving surface which 
is the boundary between two difierent media. If we write 
X = J5 ± iS and assume that the surface-charge and 

surface-current can be neglected, the six boundary-conditions 
can be expressed by saying that the three quantities of type 

must be continuous as the boundary is crossed. The equation 
of the moving boundary is expressed as before in the form (20). 
When the moving boundary is the surface of a perfect conductor 
or perfect reflector, the boundary-conditions are simply that the 
three quantities of type 

should vanish f. 

In the last two cases the boundary-conditions do not imply 
that the boundary moves normally to itself with the velocity of 
light ; in fact, the motion of the boundary can be quite arbitrary. 


EXAMPLES. 


1. The surface of discontinuity is the sphere r=^ct and the electro- 
magnetic field within this surface is expressed by the equations 


IS- -S 

o), 




where n= sin p r 4- e). 


* Abh. d. Physyohon. Ges. zu Konigsherg^ 1 Pr. (1896), p. 1 ; Ann, Phys. 
Chem. Bd. 59 (1896), p. 283. See also J. J. Thomson, Phil. Mag, (5), Vol. 45 
(1898), p, 172. 

t These conditions are equivalent to the condition that the mechanical force, 
which would act on a charge moving with the normal velocity of the surface, 
must be along the normal to the surface. Cf. Heaviside, Electrical Papers ^ 
Vol. 2, p. 514 ; Electromagnetic Theory^ Vol. 1, p. 273 ; Hertz, Electric Waves^ 
p. 257. 
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The constants 'p and q being known determine the constants A and e 
in order that the field outside the surface r=^ct may be the electrostatic 

field for which the potential is — ~ . 

(Cambr. Math. Tripos, Part II, 1904.) 


2. If the vectors E and B. are known for all points (^, y, s, i) of a 
moving surface 

F(x, 0=0 

the values of all the derivatives of E and and consequently values of E 
and jET at points not on the moving surface can generally be found pro- 
vided F does not satisfy the difterential equation 



(Havelock, Proc. London Math. JSoc. Ser. 2, VoL 2, p. 297.) 

An electromagnetic field is conjugate to an electrostatic field. 
Prove that the flow of energy in the electromagnetic field takes place along 
the lines of electric force in the electrostatic field. 


4. Let the line 00 of length c be drawn in the direction of Poynting’s 
vector at each space-time point 0 of a self-conjugate electromagnetic 
field and let OV represent a velocity v associated with the point 0. 
Prove that if this electromagnetic field is conjugate to another field 
{E, H) in which cE is the vector product of v and E^ the direction of E is 
parallel to VO. 


5. Prove that when o*, e, fM are constants, the vectors E^ E of § 6 
satisfy the differential equations 


AE- 


ejidjE 


era dE 1 , 

3^^=-gradp, 


aE- 


€/JL d^E (TfJL dE 

c2 dt 


(Maxwell.) 


Solutions of these equations for the case in which p = 0 have been given by 

0. Heaviside, Phil. Mag.^ Jan. (1889), p. 30 ; Electrical Papers^ VoL 2, 
p. 478. 

H. Poincar^, Comptes Rendus (1893), p. 1030 ; Theorie analytique de 
la propagation de la chaleur^ Oh. 8. 

J. Boussinesq, Comptes Rendus (1894), pp. 162 — 223 ; Theorie 
analytiqu^e de la chaleur, t. 2 (1903), p- 538. 

Kr. Birkeland, Archives des sciences physiques., Geneva (1895), p. 5. 

0. Tedone, Rend. Lined, Mar. 31st (1913), Jan. 18th (1914). 

M. L Pupin, Trans. Amer. Math. Roc., I'ol. 1 (1900), p. 269. 


6. Prove that the function 
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satisfies the equation 


0^24 _ 0 % 02 ^ 

0^2 0^2 ^ 0 ^ * 


7. If 24 satisfies the partial differential equation 

/024\2 /024\2 /02^\2 1 /32iY 

\0;r/ \92/ \3i/ ’ 

and we make the transformation 

"“="+-^97’ 

-024 . / 9^ 

yo=y+/^, 

where / is an arbitrary function, then 24 also satisfies the partial differential 
equation 

/024Y /024Y /024Y_1 /024Y 

w/ ^ \Wo) \^o) * 

8. Plane electromagnetic waves fall on the convex surface of an infinite 
paraboloid of revolution 4i7=a—2*, whose surface is a perfect reflector. If 
the incident waves are given by expressions of the type 




9(°) /3) 

0(^, t) ’ 


where a=y+ 22 ;, (S—x-hct, the boundary-conditions at the surface of the 
paraboloid may be satisfied by subtracting from the primary field a 
secondary field represented by expressions of a similar type but with 
7/ — {z 

a=a ^— — , 8^a-r + ct 


If waves represented by the above expressions with a—y~-iz^ cif, 

fall on the concave surface of the paraboloid the boundary-conditions at 
the surface of the paraboloid may be satisfied by supposing that the 
secondary disturbance is of the form - M' - M” M'” ^ where the fields 
M'\ M"' are represented by expressions of the above type, where «, 8 
have the values 

,, y-\-iz 

“ ^=a-r-ct, 

a"'=a?^, ^"=a+r-ct. 


respectively. With these suppositions the forces are finite at the focus, 
when / is independent of a. 


CHAPTER II 


GENERAL SURVEY OF THE DIFFERENT METHODS OF 
SOLVING THE WAVE-EQUATION 

§ 10. The object to be attained. 

It has been shown in Chapter I that the solution of 
Maxwell’s equations can be made to depend upon the solution 
of a single partial differential equation which is either the 
wave-equation flu = 0 or the equation Ait + kPu = 0 which 
is satisfied by wave-functions of the form u = z). 

The properties of functions satisfying these equations must 
accordingly be studied at some length. It is desirable, 
also, that all types of such functions should be studied 
and not merely those which admit readily of application to 
physical problems. If certain solutions of the fundamental 
equations must be rejected in the treatment of the boundary 
problems of mathematical physics, a knowledge of their 
behaviour is at any rate useful as it gives a clear indication of 
the reason why such solutions must be rejected. There is, 
however, another reason why the scope of the inquiry should 
not be restricted. The theory of wave-functions forms a 
natural extension of the theory of functions of a complex 
variable* and may consequently lead to results of great value 
for the general theory of functions. 

* This point of view is adopted, for instance, by Yolterra, Bend. Lmcei (4), 
m2,pp. 225-330,274-287 (1887); iVi,pp. 107-115, 196-202 (1889); Vi,pp. 158-165, 
291-299, 599-611, 630-640 (1889); Bend. Palermo, Yol. 3, pp. 260-272. See 
also Appell, Acta Math. t. 4, p. 313 (1884) ; Painlev4, Toulouse Ann., t. 2b 
(1888) ; Bocher, Bull. Amer. Math. Soc. Yol. 9, p. 455 (1903). The theory of 
, functions of two complex variables is closely connected with the theory of wave- 
functions. Of. H. Poincar^, Acta Math. Yol. 2 (1883) ; Yol. 22 (1898) ; H. F. 
Baker, Camh. Phil. Trans. Yol. 18 (1899), p. 431; Proc. London Math. Soc. 
Ser. 2, Yol. 1 (1903), p. 14. 
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We shall now describe briefly some of the principal methods 
of solving the wave-equation. 

§ 11. Eeduction to ordinary differential equations. 

The aim of this method is to determine elementary solutions 
of the form 

^^=/x(«)/.(^)/b(7)/4(S) (23) 

where /i,/2,/3,/4 are particular functions of their arguments and 
a, /3; 7, B are particular functions of x, y, z, t This method was 
used by D. Bernoulli in 1732 in the treatment of the vibrations 
of a hanging chain, the partial differential equation being how- 
ever in this case different. 

The general theory of elementary solutions is due to Lame* 
who transformed Laplace’s equation into curvilinear coordinates. 
For a historical account of the development of the theory we 
may refer to Prof. B6cher s book Die Beihenentwickelungen 
der Fotentialtheorie, Leipzig (1894) and to Byerly’s Fourier 
Series and Spherical Harmonics. 

A simple elementary solution of the wave-equation is 
obtained by putting /9 = y, 7=5, S = ^ ; we can then 

take 

u = 0lx-^my + nz^pt (24) 

where the constants I, m, n, p satisfy the relation 

p^ = c^{P + m^ + n^) (25) 

and can be either real or complex quantities. 

When p is a purely imaginary quantity and Z,m, n unre- 
stricted, the solution is periodic and more general periodic 
solutions may be derived from this one by summation, m, n 
being regarded as variable parameters subject to the relation 
(25). When ?, m, n are purely imaginary the solution (24) is 
appropriate for the representation of plane waves of mono- 
chromatic light, the intensity and phase of which are the same 
at all points of any plane perpendicular to the direction of 
propagation. 

* Liouville's Journal, t. 2, pp. 147-183 ; Legons sur les coordonnm 
ou/rviligneSf Paris (1859). See also E. MatMeu, Gours de ^physique matMmatique 
(1873). 
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Taking m — n=^0,l = — ~ — vso that the axis of os is in the 
direction of propagation, we may write 

Ex — Ez~ Hcc — Hy — O, Ey = a cos v (x — 

Hz — Cl cos V (x — ct) (26) 

where a is a constant. These waves will be said to be linearly 
polarised in a direction parallel to the axis of y and will be 
called homogeneous because E and H do not depend on y and z. 
It will be noticed that the constant c represents the velocity of 
propagation of a phase of the disturbance. 

A wave-function of the type 

== sin vx cos vet 


is appropriate for the representation of standing waves. Ex- 
pressions for E and H may be written down by analogy with 
the above. To obtain a representation of plane waves in a 
conducting medium, we must use a solution of 


d^u 

dx^ 


• kHc = 0 , 


where k^ has the complex value given in § 8. Putting F = 0, 
U = (0, 0, that 

E^ = Ey = H^ = H, = 0, Ey = '^ 

c 

Hz = ik ~ (27) 

where the real parts of the quantities are retained. If 

k = ^-\-i'q 

where rj is positive, the oscillations of the vector E are damped 
owing to the exponential factor 

The elementary electromagnetic fields that have just been 
found are fundamental in the theory of the reflection and 
refraction of light at a plane surface. This theory is given in 
the text-books on Physical Optics* and need not be reproduced 
here. Various attempts have been made*f to prove that any 


* See for instance, Wood’s Physical Ojptics (1911), Ch. 13 ; Jeans, Elec- 
tricity and Magnetism (1911), Ch. 18. 

t See for instance a series of papers by G. Johnstone Stoney, Phil. Mag, 
(5), Vol. 43 (1897), pp. 139, 273, 368 ; Vol. 44, pp. 98, 206 ; Brit. Assn. Reports 
(1902), p. 539 ; Phil. Mag. Feb. 1903. The idea is probably due to Stokes. 
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electromagnetic disturbance in the aether can be represented as 
the sum of a finite or infinite number of elementary disturbances 
of the character of plane-waves travelling in various directions. 
Such a representation is generally only suitable within a 
restricted domain of the variables x, y, z,t] nevertheless, it may 
sometimes be employed with advantage. 

When waves of all directions and frequencies are considered, 
the method of summation leads to Whittaker’s formula* 


rvr r^ir 

I /[a; sin a cos /3 4 - 3 / sin a sin /3 + 5 cos a — ct, a, yS] 

J oJ 0 

xsina^docdfi (28) 

for a wave-function. The case when 


/31 = ^<6 <Tr 

-0 ^>6 

has been used by Debye f in a discussion of the behaviour of 
waves of light in the vicinity of a focus. In order that an 
integral of the type (28) may represent a wave-function it is 
not necessary for the limits of integration to be those chosen. 
The limits for a may, for instance, be 0 and 0 where 0 is a 
root of an equation of type 

X sin 0 cos /3 4- y sin 9 sin ^ z cos 6 — ct = F {6). 


In order to obtain other types of elementary solutions it is 
necessary to transform our differential equations to a system of 
orthogonal coordinates {u, v, w) for which the linear element is 
given by 

, du^ dv^ dw^ 

ds-- + ( 29 ). 


If Hu, Hy} are the three components of a vector H in 
directions normal to the surfaces u — const., v = const., w = const, 
through a point {x, y, z\ the corresponding components of rot H 
are of the typej 


VW 


\W ) dw 


V 


.(30), 


* Math. Ann. (1903). See also G. N, Watson, Mess, of Math. Vol. 36 
(1906), p. 98. 

t Ann. d. Phys. Vol. 30 (1909), p. 735. 

t See for instance, H. M. Macdonald, Electric Waves, Ch, 6 ; M. Abraham, 
Math. Ann. Bd. 52, p. 81. Some very general transformation-formulae are 
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the new expression for div H is 


UVW 



and the wave-equation becomes* 


UVW 


'd f u d4>\ d f V d(i>\ d / w d<i>\ 

du\VW duj^ Zv\WU dvj ^ dw\UV dwl 


1 ^ 
& df 


(32). 


It is also sometimes advantageous to transform the wave* 
equation to a system of coordinates for which 

+ dy^ 4- dz-^ — c^dt^ = A^dp *4 B^dy^ -h C^d^^ — D^dr^, 
the wave-equation then becomes 
a fBCDdcj>\ d (CD A dcl>\ d fDABd<l>\ 

3^-4 d^j'^dvK B d^K c'di;) 

__ 9 ( ABG d^\ 

~dT\ D dr) 


§ 12. The generalisation of wave-functions, 

^ When a solution of the wave-equation has been obtained 
other solutions may be derived from it in various ways. For 
instance, the function obtained by differentiating the given 
wave- function any number of times with regard to the coordinates 
X, y, z, t, is also a wave-function. By adding together arbitrary 
constant multiples of all the wave-functions obtained in this 
way we may obtain a very general type of wave-function. 

Another method of generalisation is to make an arbitrary 
change of rectangular axes. The wave-equation is a covariant 
for such a transformation and so is a wave-function of the new 
coordinates. A number of arbitrary constants can be introduced 
into the solution in this way. We can also make a linear 
transformation of coordinates for which the expression 

dx^ -j- dy- + dz^ — c^dt^ 

contained in the papers of V. Volterra, Rend. Lincei, Ser. 4, Vol, 5, pp. 599, 
630 (1889), and J. Larmor, Camhr. Phil. Trwns. Vol. 14 (1885), p. 121. 

* Lam4, Journ. de VJ^cole Poly technique^ Cah. 23 (1833), p. 215; Legom 
sur les coordonnees curviUgnes, t. 2. A simplified proof was published by Lord 
Kelvin, Camhr. Math. Journ. Vol. 4 (1843). 
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remains unaltered in form and the preceding remarks still hold 
good. 

To illustrate this let us first of all add together two particular 
cases of Euler’s wave-function ^ y (r + ct), viz. 


1 ^ 1 

r<y-ct^ r\r + cty 


we then see that (r^ — cV) ‘ is a wave-function. Generalising 
this by writing le y —yo, z — Zo, t- to in place of x, y, z, t, 
we obtain the wave-function 


{x - XoY + (y - Vof -f- (a - ZoY - c“ (t - toY 


.(34). 


When we have obtained a wave-function involving one or 
more arbitrary parameters we noay obtain others from it by 
differentiating with regard to the parameters, or by integrating 
with regard to them after having naultiplied the expression by 
an arbitrary function of the parameters. For instance, from 
the above wave-function we may derive the more general wave- 
function 

/(t)cZt 


k 


{a> - + (y ~ -c^{t- rf 


.(35), 


where the integration is between constant limits. It is not, 
however, really necessary for the limits to be constant, we may 


for instance take them to be — oo and t — (r 4* n), where 

c 


-i- j/o^ + V- resulting integral is then a wave-function 

provided /(r) behaves in a suitable manner. If we take /(r) = 1, 
we obtain a wave-function 


u = 


:log 


R — r — Vo 


where 


R ^R+r+u 

R‘={x- a!o)“ + (y — yoY + k- ^oY- 


.(36), 


This function is independent of t and so must be a solution of 
Laplace’s equation Am = 0. It is closely connected with the 
function used on p. 3 of Basset’s Hydrodynamics, Vol. 2. 
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§ 13. Transformations. 

In addition to the linear transformations that have already 
been mentioned there are certain other transformations which 
enable us to pass from one wave-function to another^. 

The first transformation, which is analogous to inversion, is 
defined by the equations 


where 



if ^2 



It is easy to verify that if f{x, y, z, t) is a wave-function, 
then 


IKi- 1 . 5 - 1 ) 


is also a wave-function f. 

as 

y 




I - cH^ 4- 1 

z — ct' ^ z—ct' ^ 2{z’-ct) ' 2c(z — ct) 


.(39). 


It is easy to verify that if f(x, y, z, t) is a wave-function, 
then 

CO y 1 


' — cV \_z 




^z-'Ct' z — ct' 2{z — ct)' 2c (z ■ 
is also a wave-function. Since is a wave-function we 

may deduce in this way that 

1 

e z-ct 

z — ct 

is a wave-function. 

It should be noticed that if we put 4 t = ^ — c^, a — z-\-ct, 
a function of the type 

u = F(x, y, t) e"®* (41) 


* For a general account of these transformations see a paper by the author, 
Proc. London Math, Soc. Ser. 2, Vol. 7 (1908). 

t This of course is a simple generalisation of Kelvin’s theorem for Laplace’s 
equation. The generalisation to the corresponding equation in n variables is 
mentioned by Bdeher, Bull, of the Amer, Math. Soc. Vol. 9 (1903), p. 459. 

t Proc. Londm Math. Soc. Ser. 2, Vol. 7 (1909). This transformation is 
equivalent to a conformal transformation of a space of four dimensions which 
was discovered by Cremona. Cf. Darboux, Logons sur les systemes orthogonaux 
et les coordonnees curviligneSj Paris (1910). 
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is a wave-function if F satisfies the differential equation 


d^F , d^F dF 

dy‘^ dr 

The wave-function we have just obtained indicates that 


1 a-g+i/ 

te-~i;r^F (43) 


is a solution of the above equation. This solution is fundamental 
in the theory of the conduction of heat. It is evident that any 
solution of the equation of the conduction of heat in two 
dimensions can be used to construct a wave-function. 

Our second transformation theorem for the wave-equation 
also tells us that if F{Xy y, r) is a solution of the equation (42) 
the function 


1 - 
- e 


r 



y 



,(44) 


is also a solution. This result is due to J. Brill* and Appellf. 

The corresponding theorems for the two-dimensional wave- 
equation 

da? ^ df & d? '' '' 

are first that if fix, y, t) is a wave-function and 


the function 


s^ = x^+y^-cH\ 

y i\ 

s-^ {?’ sV 


( 46 ) 


is also a wave-function. This is equivalent to Lord Kelvin's 
theorem for Laplace's equation if we simply replace ict by z. 

The second theorem is that if f{x, y, t) is a wave-function, 
then 






2(2/- cO’ 


? + l ~ 
2c{y-ct)^ 


( 47 ) 


is also a wave-function. Writing 4!T — y — ct, cr — y + ct as 
before we find that u = e~''F (co, r) is a wave-function if 


da? dr 


( 48 ). 


* Messenger of Mathematics (1891). 
t LiouvilWs Journal, Ser. 4, t. 8 (1894). 
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The second theorem can now he used to show that 

= e ^ (49) 

is a solution of this equation and that if F {x, r) is one solution 
the function 

T-^e-TrF(^, -i) (50) 

is also a solution. This theorem is likewise due to Brill and 
Appell, it can evidently be generalised to the equation in 
n variables. 


EXAMPLES. 

1. Prove that it is ]30ssible for a train of plane electric waves 
to travel along an infinite isolated slab of dielectric material without 
being dissipated by spreading out into the adjacent empty space. Show 
that if 2a is the thickness and K the inductivity of the slab, the velocity 
of propagation of such waves of length X along the slab, when polarised so 
that their magnetic vector is parallel to it, is 




where 6 is the lowest real or the pure imaginary root of the equation 

■ 


(Larmor, Cambr. Math. Tripos, Part II, 1906.) 


2. Prove that with the notation of § 13, the function 
1 ^4-^0 * 

{x — xqY + (y -yo )^+ {z - - iof' Jssq + ^^0 +yyo + 

is a wave-function, zq, to being arbitrary constants. 


3. 


If F{x^ y, 2 , t) is a solution of the wave-equation, the function 



cr^T tret 
__ 


^''0. ?. i. I) 


dr 


is, under suitable conditions, a solution of the equation 




erdV 

C 0 ^ ' 


The quantity has the meaning assigned to it in § 13 and is supposed in 
the present case to be negative. 


B. 


3 
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4. The function 

F=- 4 = I e-^'‘d'K 

/s/tT J QQ 

satisfies the equation 

?!r 

3^ 3y2 0;! ’ 

of the conduction of heat in two dimensions : it is zero over a semi-infinite 
line which covers part of the axis of x and moves in the positive direction 
with uniform velocity -v. The isothermal lines at any given instant are 
confocal parabolas. 


5. Prove that, under suitable limitations^ the function 




t 


<u 


is a solution of 




Obtain in this way the particular solution 

- {x-^-iyY^ . 


p2 = .q72 


6. Prove that if the integral 


/: 


-^^da 


^2 4-^2 _j. 4^2 _ 4^^ 


= F 


satisfies 
and if the integral 


^ dV 

dx^ dy^ 0^ ’ 




^ da 


0 iJx^+Aa^’-Aat 


V 


satisfies 


0F 

dsc^ ~ 0{2 3f • 


7. Prove that if p^=x^+y^, the integral 

F= f °° sin^(g+a)(fg 
J p+p» {y-y^f 


0*7 027 


satisfies the equation 
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§ 14. The elementary solutions. 

If we make Laplace’s transformation 

x = rsm6 cos(f>, y = r sin 5 sin 0 , ...(51), 

the equation Au + khi = 0 becomes 


d^u 2du 1 ^ A* p 9^\ 1 

0^2 r sin ^ ^ ^W)'^ r^sin^i? 

This is satisfied by a function of the form 
u = It (r) @ (6) <I> (<f)) 


+ kht = 0...(52). 


if 


ocp- 


.(53), 


dr^ 


2dR \k-‘ , ,/ 


sin’^ 


0 = 0... (54), 


R = 0 


.(55). 


The first equation is satisfied by = cos (m(f> + a), the 
second by P^’"(cos^) and Qn^{cos6), where these are the 
associated Legendre functions. The third equation may be 
written in Bessel’s form 


d^w 1 dw 
dr^ r dr 



w = 0 


(56), 


where w = r^R, and is satisfied by Jn+i and J_ (kr). 

In these solutions m and n can have any constant values. 
It should be noticed that when n H- ^ is an integer the Bessel 
. functions that have just been written down are not independent 
and the second solution (kr) of Bessel’s equation must be 

used. 
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In dealing with a problem such as the effect of an obstacle 
on a train of electric waves, the secondary waves sent out from 
the obstacle must have the character of diverging waves at a 
great distance from the obstacle. In the case n=0 the diSeren- 
tial equation for R is satisfied by 

(57), 

and if the real part of k has the same sign as co when the 
electric and magnetic forces are the real parts of expressions 
of the form we can obtain a solution appropriate for the 

representation of a diverging wave by taking the positive sign, 
for then we have a function of the form 

r 

Neither of the given solutions of BesseFs equation has the 
required form in fact 

/j (^) = sin kr, J-:^(kr) = ^~ cos kr. 

We may, however, obtain solutions of the form (57) by taking a 
suitable combination of the preceding solutions. 

In the case of electromagnetic fields in the free aether the 
physical interpretation of the elementary wave-functions when 
n is zero is as follows^ : 


- cos k (r — ct) 
r 

- cos fc (r + ct) 
r 

i cos kr . cos kct 

T 


Progressive divergent waves. 
Progressive convergent waves. 

Standing forced waves, source at origin. 


- sin Air. cos kct Standing free waves. 

To obtain the solution of (56) appropriate for divergent 
waves when n has any value we writef 


* A fuller discussion is given by A. Sommerfeld, Jahresbericht der deutsch. 
math. Verein, Bd. 21 (1913). 

t The theory is due to Stokes, Phil. Trans, Vol. 158 (1868), p. 447; 
Collected Papers, Yol. 4, p. 321. See also Rayleigh’s Sound, Vol. 2, p. 304.* 
It should be mentioned that different notations are used by different writers. 
This is the notation used by Debye, Ann. d. Phys. Vol. 30 (1909), p. 57. 
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Vn (OO) = ‘^p'n («) - *%« («) 




'JTX 

'2 


L.(58). 


{oc) — (^) + i%n (^) 

These new functions fn are connected with HankeFs 
cylindrical functions'^ by the relations 


Vn 




r, 






When the real part of x is large and positive we have the 
asymptotic expansion 


%(«;) = (- 


1 + 




2x 1 ! 

1 {n — 1) n (n + 1) (n -}- 2) 


-f 


...(59). 


{2x)- 2 ! 

The series terminates when n is an integer and then gives 
a true representation of the function. To get tn(^) we change 
the sign of L Various other notations have been used for the 
solution of BesseFs equation that is suitable for the representa- 
tion of diverging waves. Lamb uses Du(x), {v—n + ^) to denote 
the solution of equation (56) which has the asymptotic value 

D,im) = f— (-+j) , 

yirxj 

while many other English writers use Kp{ix) to denote the 
solution with the asymptotic form 

The following foi-mulae will be found useful : 

(/»')”+’■ 


■^n{kr)-- 


X 


1 - 


1.3...(2n+l) 
(kry- ^ 


{kr)* 


2(2n+3) ^ 2 . 4 . (2w + 3) (2n, + 5) 

See Nielsen’s JSaiidhuch dev Cylinderfunktiomn^ p. 16. 


••] (60), 



38 


POLAR COORDINATES 


[CH. 




1 - 


i n (n + 1) 
2kr 1 
(n—l)n (n + 1) (w + 2) 




{2n)l 


...(61). 


{2kr'f' 21 , ... V ■(2fe")”.n!j 

We have written down the last term in the series on the 
supposition that n is an integer. In this case 

■fn (kr) = i + (- e*’’], | r \ large . . .(62), 


(2n 4- 1) = (w + 1) {x) -Ti'fn+i (^) ■ ■ .(63), 

(2n + = as («) + ■^m+i (<»)] (64f), 

Ky(x)=f e“*“*“coshi/a.da (65), 

JO 


"" 2>'r(i)r(y+ F ) J 0 i 

( 66 ), 


/: 




{n+f) 

{n = p) 


l?^>0,j9>0 (67). 


2?i + 1 

In the last formula n and p are supposed to be integers. 
For further properties of Bessel functions the reader should 
consult Gray and Mathews’ Treatise on Bessel Functions, 
Nielsen’s Handbuch der Cylinderfunktionen, and Whittaker’s 
Analysis. Tables are given in the first work and in Jahnke 
and Emde’s Fwiktionentafeln, Leipzig (Teubner). A fev?- ad- 
ditions to the tables have been made recently by J. W. Nicholson, 
Proc. London Math. Soc. Ser. 2, Vol. 11, p. 104 ; Dinnik, Archiv 
der Mathematik und Physik (3), Bd. 20, Heft 3, 1912 ; J. R. Airey, 
Phil. Mag. Vol. 22 (1911), p. 85, Brit. Ass. Reports (1911); 
A. Lodge, Brit. Ass. Reports (1909) ; J. G. Isherwood, Manchester 
Memoirs (1904). 

The best definitions of the generalised Legendre functions 
for unrestricted values of m and n are those given by Hobson*. 

* FMl. Trans. A, Vol. 187 (1896), pp. 443 — 531. E. W. Barnes has recently 
given new definitions of the functions as integrals involving Gamma Punctions 
which make it possible for the principal formulae to be proved very quickly. 
His definition of (x) differs from that of Hobson by a numerical factor which 
becomes rather troublesome when n is an integer and m is not. 
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We are interested here in the case -when the variable is cos0 
and 0 is a real angle, such that 0 < 0 < tt. We may then put, 
with the usual notation of the Gamma and hypergeometric 
functions, 


e 


Fn^(cQS 9) = cot^^ rz^ + l; 1-m; sin^^ 

( 68 ), 

Qn^ (cos 6) = 


2 sin (n -h m) TT F (1 — m) 

|cos {n 4- m)7r cot"^ | f{^ n-^1] 1 — m ; 

9 f 

— tan’” ^ jP ?i, n + 1, 1 — m, cos^ 


• ^ 
sm-5 


.(69), 


P„-“(COS 0 ) = 


cosec™ 0 


_r- 

-f (2 cos <j{> — 2 cos 


cos (n + I) ^ 


i2 (w 4- -1) > 0 

When m is a positive integer, we have 

^ r7m 

Qrr(<c) = (i-c<^y (a:) 


.(70). 


•(71). 


and when is a positive integer 


Pn(«!) = 


1.3...(2n-l) 

1.2. ..n 


L 2(2 


(ra-1) 


(2w-l) 


w(w-l)(n-2)(TO-3) 

2.4(2n-l)(2n-3) 

Pn (cos 0) = cos”0 — — — COS’‘-- 0 . sm^0 


-•••] (72), 


+ 


2= 

w (n — 1) (n. — 2) (n — 3) 
22.4= 


cos""'*^ . sin'*^ (73), 


1 D / a\ (- 1)“ /'I' 


■(M). 


05 ). 
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The last two formulae illustrate the method of deriving 
more complicated wave-functions from simple ones by differen- 
tiation. The functions 


- and 
r 



r -f ^ 


are in fact solutions of (52) when k = 0. 
The formula 


/ 


+1 


Pv 


' (x) dx=^0 {n^v) 

__ 2 (?H- m ) ! 

2n + 1' (71 — m) 1 


.(76), 


in which m, n, v are positive integers or zero, enables the 
coefficients in an expansion of a function in series of functions 
Pn ^ (^) to be determined by a simple integration. 

For further properties of the Legendre functions we must 
refer to Heine’s Kugelfunktionen, Byerly’s Fourier Series and 
Spherical Harmonics, Whittaker’s Aiialysis, Nielsen’s Handhuch 
der Gylinderfunktionen and Theorie desfonctions mdtaspheriques, 
and to memoirs by E. W. Hobson^ and E. W. Barnes f. 

Tables of the Legendre functions have been published by 
J. W. L. GlaisherJ, J. Perry§ and A. Lodge||; some tables of 
the functions P-^F (p) have been given by H. TallquistH'. 

For the history of the functions of Legendre and Bessel 
the reader should consult the article by A. Wangerin in the 
Encyklopadie der Mathematischen Wissenschaften, Bd. II. 1, 
Heft 5 (1904), p. 695. 


§ 15. Relations between various solutions. 

We have already remarked that when a wave-function or 
a solution of equation (52) involving arbitrary constants has 
been found, other solutions may be derived from it by the 
method of summation or integration. By choosing our sum or 
integral so that it represents certain simple solutions of the 

* Proc. London Math. Soc. Ser. 1, Vol. 22 (1891), p- 431, and op. cit. 

t Quarterly Journal, Vol. 39 (1908), p. 97. 

X Prit, Ass. Peport (1879). 

§ Phil. Mag., Dec. (1891). See also Byerly, loc. cit. 

jl Phil. Trans. A, Vol. 203 (1904). 

IT Acta Societatis Fennicae, Vols. 32, 33. 
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fundamental equation, a number of important identities may 
be obtained. A few formulae will be written down to illustrate 
this^ 

If R- = — 2rri cos 9, 

R Vu^i) »=o 

(Heine and Hobson) (77), 

4 e"**-® =4“ 2 (2n + 1) (Ar) (M) Pn (cos 0) (n > r) 

R ^W. = 0 /Hrn\ 

(Heine and Macdonald) (78), 


(fo-) Pn (COS 0) = ^ + J (*n) ~ ^ 

(79). 

We may illustrate the peculiar behaviour of certain definite 
integral solutions of our fundamental equation by the following 
example, in which k is supposed to be real and positive. 


Let 


rm 

fm{a;)= I 
Jo 


cos xt X (t) dt, 

rm 

where % (t) is a function such that I \xO')\dt is convergent, 

J 0 

then it may be proved by means of Fourier’s double integral 
theorem thatf 

1 f (n) dn if) {k > m) 

" J — 00 * • 1 

=/* (’") (^ m). 


Hence 


u 


_1 r dnkR 
'ttJ-oc B U{ri)dH 


is a solution of (52) which reduces to either (r) or (r) 
when ^ = 0. Solutions of (52) which are derived from the 
elementary solutions by integrating with regard to n have 
been employed by H. W. MarchJ. He makes use of an inversion- 
formula § 

* Some very general formulae are given by L. Gegenbauer, Monatsh, f. Math. 
Bd. 10 (1899), p. 189. 

t This equation is obtained in a different manner by G. H. Hardy, Proc. 
London Math. Soc. (2), Vol. 7 (1909), p. 445. 

:|: Ann. d. PhysiJc, Bd. 37 (1912). See also H. Poincar^, Comptes Eendusy 
t. 154 (1912), p. 795 ; W. v. Bybcyhski, Am. d. Phys. Bd. 41 (1913). 

§ This is in some respects analogous to the inversion-formula given by 
F. Or. Mehler, Math. Ann. Bd. 18 (1881), p. 161. 
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/((?) = [ P„ _ I (cos 6) ^}r (a) ada 

r I 

’/"(«)=] (cos 7)/(7) Sin 7^7 j 

sin aj3 . d/3 


^ I . cos aTT 2 
S^_^{cos 7 )= + - 


ttTT COS 


7 TT j y {2 (cos 7 — COS /?)}^ 


[CH. 

■(80), 

.(81). 


It is sometimes instructive to find how a wave-function, 
depending on an arbitrary function, can be expressed in terms 
of elementary wave-functions. Now in the second example 
of § 5 the electric and magnetic forces are all of the form 


* 70 ] («*). 


or are the sums of terms of this form. Consequently, a function 
of this type may be expected to be a wave-function and it is 
easy to verify that this is the case*. We may now deduce that 


(83) 


ff 

is a solution of (52) and consequently it follows that tan’" 2 

Q 

and cof" ^ solutions of equation (54) when n = 0. We 
have in fact 

1 


Pc’" (cos 6) = 




cos miT . cof" 2 2 ) • 

In the last formula m must not be zero or a negative integer. 


§ 16. The convergence of series of elementary solutions. 
When h = Q our fundamental equation (52) reduces to 
Laplace’s equation = 0 and we have the familiar elementary 
solutions 


* The theorem also follows immediately from a result given by A. B. 
Porsyth, Messenger of Mathematics (1898), p. 114, and B. W. Hobson, loc. eit. 
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r'^Pn”" (cos 6) cos — (j>c), r” Qn’” (cos 6) cos m((f}— (^o) 

^ (cos 6) cos ?92 (<^ - <f>o), (cos 6) cos m (<^ - </)o) 

(84). 

A pair of series of the type 

© n /fy \ 71+1 

2(“) p„(0,<^) 

which converge when r=a are suitable for representing harmonic 
functions inside and outside the sphere r = a because the first 
converges absolutely when r< a and the second when r > a. 

The case in which 0 is very similar. When is large 
the function ^jpn (kr) may be replaced by 

(^.^)n+l 

1.3,..(2n 4-1) 
and so a series of the form 


2 yjr^ (kr)fn {0, 

n-O 

converges like a power series. Again, when kr is real, we have 

I ^nQcr )\ 2 = 1 + - ^^(9^ + 1) ~p^ 

+ |^(»^-l)n(«+l)(« + 2).^, + ... 




n being a positive integer. It is clear from this equation that 
I i decreases as r increases, hence if a series of the form 


2 Ukr)fn(0.<l>) 

n-O 


converges absolutely for any value of r it converges absolutely 
for all greater values of r. 

For a discussion of the convei'gence of series of spherical 
harmonics we may refer to 0. Neumann’s book Ueher die nach 
Kreis-, Kugel-und Cylinder ’•Funktionenfortsclireitenden Entwicke- 
lungen, Leipzig (1881); to Heine’s Kugelfiinktionen, Bd. 1, p. 435, 
Bd. 2, p. 361, and to papers by U. Dini, Ann. di Mat. (2), t. 6 
(1874); H. Poincard, Oomptes Rendus, t. 118 (1894), p. 497; 
S. Chapman, Quarterly Journal, Vol. 43 (1912), p. 1 ; T. H. 
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Gronwall, Math. Ami. VoL 74 (1913), VoL 75 (1914), Gomptes 
Mendus Amer. Trans. Jan. (1914), VoL 15; C. Jordan, 

Gouts d' Analyse, 2nd ed., Vol. 2, p. 252 ; and B. H. Camp, Bull, 
of the Amer. Math. Soc. Vol. 18 (1912), p. 236. The con- 
vergence of series of Legendre polynomials has been discussed 
very thoroughly by G. Darboux, Liouvilles Jotirnal (2), t. 19 
(1874), p. 1; (3), t. 4, p. 393; 0. Blumenthal, Dissertation, 
Gottingen (1898); E. W. Hobson, Ptog. London Math. Soc. (2), 
Vol. 7 (1909) ; L. Fejer, 3fath. Ann. Bd. 67, p. 76; D. Jackson, 
Amer. Trails. Vol. 13 (1912). 

§ 17. The scattering of plane homogeneous electro- 
magnetic waves by a spherical obstacle. 

The effect of small particles in scattering incident radiation 
has been discussed very thoroughly by Lord Rayleigh^ who has 
used it as the basis of a mathematical theory of the blue colour 
of the sky. The action of a single spherical particle is of 
fundamental importance and so the electromagnetic theory 
of the scattering of light by a dielectric sphere has been worked 
out by Lord Rayleigh f, Prof. Lovet and other writers. This 
theory can also be developed so as to cover the mathematical 
theory of the rainbow. 

The more general theory of the scattering of incident radia- 
tion by a spherical obstacle § with arbitrary optical properties || 
admits of some very interesting applications in the study of the 
colours exhibited by metal glasses, metallic films and colloidal 
solutions or suspensions of metals. The electromagnetic theory 
of these colours has been developed by J. Maxwell Garnett IT, 
G. Mie* * * § **, R. Gansf-|- and Happelff, who have considered 

* mu. Mag. Vol. 41 (1871), pp. 107, 274, 447; Vol. 12 (1881), p. 81; 
Collected Papers, Vol. 1, pp. 87, 104, 518. 

t Phil. Mag. Vol. 44 (1897), pp. 28 — 52; Collected Papers, Vol. 4, p. 321; 
Proc. Poy. Soc. Vol. 84 (1910), p. 25 ; Vol. 90 (1914), p. 219. 

X Proc. London Math. Soc. Vol. 30 (1899), p. 308. 

§ The work of Stokes, Camb. Trans. Vol. 9 (1849), p. 1, with later appli- 
cations, Collected Papers, Vol. 4, and of L. Lorenz, Wied. Ann. Vol. 2 (1880), 
p. 70, opened up the subject. 

|( The case of small conductivity was discussed by G-. W. Walker, Quart. 
Journ. Vol. 30 (1899), p. 204; Vol. 31 (1900), p. 36. 

IF Phil. Trans. A, Vol. 203 (1904), p. 385; Vol. 205 (1905), p. 237. 

Ann. d. Phys. Vol. 25 (1908), p. 377. 

ft Ibid. Vol. 29 (1909), p. 280 ; Vol. 37 (1912), p. 881. 
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the cases of spheres and ellipsoids endowed with the oj 
constants e, o* * * § ** . 

The particular case of a perfectly conducting sphere 
worked out by J. J. Thomson^ and has been discusse 
greater detail by J. W. Nicholson f. 

The problem is also of importance in connection with 
theory of comets’ tails which has been developed by E 
FitzGerald and AiThenius§. The pressure of light ( 
perfectly conducting spherical obstacle has accordingly 
calculated by K. Schwarzschild|| and J. W. Nicholson 
more general case of a sphere with the optical constants €, 
has been treated very fully by P. Debye 

Let us assume that the electric and magnetic forces E 
at any point of space are the real parts of vectors E, H oi 
form where -d is a complex quantity independent 

We then \vrite 


H±ivE=.Me^\ where = 

fJLCO 

and we find that the differential equations satisfied by M 
medium whose optical constants are e, a can be writt( 
the form 


+ hMy = 

+ hl^I^ = 


sin 6 
1 


~ {r sin e . Mi) - ~ {rMe) 


dMr d , . d ,, s' 

- 5 - (r sin 0 . Mi) 


r sin 6 \_d^ dr 
1 


— (rM ) — 


where k = or if = 
c 




.(S 


* Receiit Researches^ p. 437. 

t Froc. London Math, Soc. (2), Yol. 9 (1910), p. 67; Yol. 11 (1912), ] 

X Scientific Writings, pp. 108, 531. 

§ Fhys. Zeitschr. Yol. 2{1901),pp. 81—97 ; Das Werden der Welten, I 
(1907), p. 85. 

11 Sitzungsher. d, Kgl. Bayer, Akad, d. Wiss. Yol. 31 (1901), p. 293. 

IF British Association Reports (1910), p. 544; Monthly Notices of the 
Astronomical Society, Yol. 70, p. 544. See also J. Proudman, lUd. Yol. 73 \ 
p. 535. 

** Ann. d, Physik, Yol. 30 (1909), p. 57. 
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These equations may be satisfied by putting 

^ r drd0 r sin 6 
1 02 (rXl) ^ k i 

rsin^ oro^ r 

where the function fl = Z7 ± iF is a solution, of (52). 

The electric and magnetic forces may be derived at once 
from these expressions by equating the ambiguous and un- 
ambiguous parts as explained in § 2. 

We shall now assume that the incident wave of plane 
homogeneous monochromatic polarised light is represented by 
{Mr, Mq, M^) = e^^rcos6±i4> cos 0, ±i) (89), 

the electric vector being parallel to the axis of y. 

The corresponding function flo is obtained by solving the 
equation 

^ (rUo) + = sin ^ ^ , 

We easily find that 

rflo = [cosec 0 . +/i (0, </>) -h f, {0, <f>) 

(90). 

Choosing the unknown functions so that rOo is finite for ^ = 0 
and 0 = TT we obtain finally 

r 0 0 ~ 

rDo = 2 cosec 0 . — cot ^ — tan ^ 

A A 

( 91 > 

We may now assume an expansion for rllo of the form 
rOo = k'~^ 2 dn'^n (kr) (cos 0). 

n=l 

To determine the coefficients we multiply by sin 6 and 
differentiate both sides of the equation. Then since 

j 

^ [sin 6 . P„i (cos ^)] = - n (m + 1) sin ^ . Pn (cos 6) ( 92 ), 
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the coefficients may be determined at once with the aid of 
Lord Rayleigh's expansion* 

ikr sin 6 . = 2 (2u +1) (k'^) Pn (cos 6) sin 6 

w=0 

(93). 

We thus obtain 

rOo - k-^ 2 fn (kr) (cos (94). 

n=i n{n± i.) 


Now let Oi, fla be the functions from which the electric 
and magnetic forces in the scattered light and transmitted 
light may be derived respectively. The appropriate forms are 
given by equations of the type 


r f7i = ^ An^n (kr) Pn^ (cos 6) cos ^ 

71=1 

rFi = 2 Bn^n Qcr) P„* (cos 9) sin <j) 

7i = l 


.(95), 


rJJ^— ^ Cnirn(hr)Pn^(Q0Sd)C0B<j) 

71 = 1 

r Fa = 2 Dn ‘^nikr) Pn^ (cOS d) sin (j) 

71 = 1 

where h, t] are the values of k, v respectively within the sphere. 

It is easy to deduce from (88) that the tangential components 
of the electric and magnetic forces are continuous in crossing 
the sphere r = a, if when r = a 

-[C^o + f7->] = -(7. 



A(F„ + F)=AF 

* Theory of Sounds Vol. 2, p. 272. The expansion was also obtained 
independently by Heine, Kugelfunhtioneji (1878), Bd. 1, p. 82. 
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(*«) 

kAn^n' (ka) - hOnfn (ha) = | 

1 2^7 -4- 1 

(ka) - hDnfn (ha) = t fn (ka) 


[CH. 


1 .. 


92 (?^ + 1) 

2 n 4- 1 


Solving these we get 


...( 98 ). 


'I^n (ka) 


2ra+l Fn 


A _ «n.+l /g -2 ^ (7 ^ «»+l 

n(n + !)(?„’ A* n(n+ !)(?„’ 


5„ = i»+i 


2w + l 


where 


n (ti + 1) ^ 


-^n — 7 7 ^ 


____ 2n + l 

hk n(7i-\-l)Ln' 


Hn = r)y}rn (ha) (ha) — i/i/r,/ (ka) (hay 
En = (ka) yfr^ (ha) — 9?'v/r,/ (fo) yjpn (ha) 

Fn = (^a) i/rn (ka) yjtn (ka) ^ . . .(99). 

Gn = V^n (ha) (ka) - z^5'n (ka) yfrn (ha) 

Ln = vy^n (hcL) ^n (ka) (ha) 

We can prove that our series all converge absolutely and 
uniformly at about the same rate as a power series of the form 


O)'’' 


1.3...(29^+l)’ 


kcv^ 

where x is either &r or — . The proof depends on the fact 
that when n is large we have approximately 


rpTl’^-l 


, / \ f / / \ (n’¥l)x^ 


(a>) 1 . 3 ... (2n “ 1) 

tv 


?« (a:) = - 1 . 3 . . . (2n - 1) niaT’^-^e-^. 
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§ 18. Free damped vibrations for the space outside the 
sphere. 

It should be noticed that some of the terms of our series 
become infinite when either = 0 or X^ = 0: fortunately, 
however, the roots of these equations turn out to be complex 
and so when k is real no values of k need be excluded from the 
discussion. The damped vibrations determined by the roots of 
the equations = 0 may be distinguished as the electric 
vibrations, those determined by equations of type = 0 as the 
magnetic vibrations. Some of the roots of the equations have 
been calculated for the case of a totally reflecting sphere by 
Sir J. J. Thomson*, who finds that the roots are all complex. 
The vibrations for the space inside a totally reflecting sphere 
have been discussed by Prof. J. W. Nicholson f, those for the 
space between two concentric spheres by Sir J. J. Thomson]:, 
Sir Joseph Larmor§, Prof. H, M. Macdonald ll and A. LampalT. 

P. Debye, who has calculated some of the roots for a case of 
a dielectric sphere, finds that the roots are complex and of two 
types. When the index of refraction is large, the imaginary 
part of a root of the first type varies very little with the index 
of refraction N and approaches a limit different from zero when 
N oo . If on the other hand p is a root of the second type, 
Np tends to a finite real limit, viz. a root of (-^p) = 0, when 
and so the imaginary part of a root of this type must 
be very small when N is large. 

The vibrations belonging to the space outside a sphere 
must be in all cases damped on account of the loss of energy by 
radiation; when the refractive indices of the outside medium 
and sphere are very nearly equal, they are clearly very strongly 
damped ; thus it is only when the refractive index is large that 
some of them are durable. It is doubtful whether a substance 
exists which has a large refractive index and does not absorb 
light to a marked extent. 

* Proc. London Math. Soc. Ser. 1, VoL 15 (1884), p. 197 ; Recent Researches, 
p. 361. 

t Phil. Mag. 1906, p- 703. 

:j; Recent Researches, p. 373. 

§ Proc. London Math. Soc. Ser. 1, Vol. 26 (1894), p. 119. 

jl Electric Waves, Chapters 6-7. IT Wien. Ber. 112 (1903), p. 37. 
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It will be seen later that the characteristic vibrations play 
an important part in determining the size of the sphere on 
which the pressure of a given type of incident radiation has 
a maximum value. 

Prof. Love* has used the solutions corresponding to the 
characteristic vibrations to discuss the mode of decay of an 
arbitrary initial disturbance. He makes use, in fact, of the 
functions where k is one of the roots of one of the 

equations = 0, = 0 ; only the sphere is treated as a perfect 
conductor. 

This method can easily be extended so as to provide us 
with a method of discussing the problem of the scattering of an 
arbitrary primary disturbance by a spherical obstacle. In this 
method we assume that the total disturbance outside the sphere 
can be represented by 

r?7 = 2 2 m,p^n Pn^ (cOS 6) COS m {<f> - 4>o) 1 

w=0 w-Ojp I 

rV=’t t t Bn, (VO (cos 0) sin m - <f>o) I 

(«). 

where the A?/s are roots of one of the equations of the type 
(ka) ~ v^n {ko) 'f,, (ha) = 0) 
vyjrn {ha) ^Jrn (ka) — ri'y^n {ka) (ha) = Oj * * 

The coefficients must then be chosen so that this total disturb- 
ance has the same character as the primary disturbance at its 
singularities outside the sphere and at an infinite distance, 
taking into account of course the presence of diverging waves 
from the spherical obstacle the effect of which is, however, 
negligible at infinity. 

The field inside the sphere is represented by equations 
similar to (co), one with hp written in place of kp. The boundary 
conditions are satisfied in virtue of (100). 

If a series of type (co) should fail to represent the disturbance 
outside the sphere, it may be necessary to add terms cor- 
responding to the fi:ee characteristic vibrations. These are of 
the form (co) with the function written in place of and 
numbers kn determined by equations of type Ln = 0, (?n = 0. 

* Proc. London Math. Soc. Ser. 2, Vol. 2 (1904), p. 88. 
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§ 19. The case of a very small obstacle. 

When a, the radius of the sphere, is very small compared 
with the wave-length A of the incident radiation, we may treat 
ka and ha as small quantities. We may then obtain some idea 
of the relative magnitudes of the different coefficients in our 
series by using the expansions (60) and (61). 

It is easy to see that the values of An, Bn decrease very 
rapidly in absolute magnitude as n increases. The disturbance 
radiated from the sphere can consequently be represented 
approximately by superposing a small number of partial waves, 
the effect of the others being negligible. 

Eemembering that vh = rjk, we find that when H is finite 


rr - A^) 

^'^(2^+3).P.3^..(2?2 + 1)2’ 

3 ^ {n + 1 ) (^^" - 

P.3^..(2n + 1)2 



Ln ivh^k-^' 


,_2 ('^+ l)k‘^-{-nh^ 
+ 1 


•ika 


) 


and it is easy to see that all our series converge. 

It appears from these expressions that the ? 2 th magnetic 
wave, i.e. the disturbance due to the 72th term in the expansion 
for Ui, is of the same order of magnitude as the (n + l)th 
electric wave, i.e. the disturbance due to the (n + l)th term in 
the expansion of V^. This is in sharp contrast with the result 
obtained by Sir J. J. Thomson for the case of the totally 
reflecting sphere wherein the 7^th electric wave and the nth. 
magnetic wave are of the same order of magnitude. 

The first electric wave is clearly of chief importance and 
Mie has proposed to call this Rayleigh's radiation. We easily 
find that 






-1 k^{k^-hf) 
2.32 3A-"4-2A- 


...( 102 ). 


. The following diagrams, which are taken from Mie’s paper, 
indicate the character of the electric lines of force for the first 
four partial ^dbrations of each type. For the magnetic waves 

4—2 
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Er = 0 and so the electric lines of force are spherical curves. 
In the case of the electric waves the lines of force lie on 
certain cones and the diagrams represent the intersections of 
a sphere with these cones, the vertices- of the cones being at 
the centre of the sphere*. 

§ 20. Polarisation of the scattered light. 

Let us now look for cases when the light scattered by the 
sphere is linearly polarised. It is easy to see that Eq and M4, 

These conditions are 


both vanish when = 0 and = 0. 
ou 


9 ^ 

both satisfied by <jf> = 0, i.e. when the observer looks in a direction 
at right angles to the electric vibration in the incident wave 
(Fig. 3). 

It appears from the figure that the component of the electric 
vector of the scattered light, which 
is at right angles to the direction 
in which the observer is looking, 
is parallel to the electric vector 
in the incident wave. In a similar 
way it is found that E^, and Mq 

both vanish when </> = ± ? , i.e. 

when the observer looks in a 
direction at right angles to the 
magnetic vibration in the incident 
wave. The magnetic vibrations in the incident and scattered 
waves are now found to be parallel. 

The experiments of Steubingf with different kinds of col- 
loidal gold solutions have shown that when the solution is 
illuminated with polarised light and viewed in the manner 
described, there is always a small quantity of unpolarised light 
sent out from the particles, but the greater portion of the 
scattered light is polarised in the way the theory requires. 
The slight disagreement between the theory and observations 
is attributed to the fact that the metallic particles are probably 

* The author is indebted to the publisher of the Annalender P/ij/stfc, Herr Johann 
Ambrosius Barth, for permission to reproduce the figures on pp. 52, 59 and 64. 
t Dusertation, Greifswald (1908) ; Ann. d. Phys. Vol. 26 (1908), p. 329. 
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not all spheres*, that some may have developed into crystals 
perhaps of octahedral form. The mathematical theory of the 
scattering of waves has not yet been fully developed. The 
problem is, however, one of great importance in meteorological 
optics. The case of a regular distribution of atoms or molecules 
has recently been brought into prominence f by experimental 
work on the scattering of Rontgen rays by a crystal Approxi- 
mate mathematical theories have been given by several writers §. 


§ 21. Intensity of the scattered light. 

When Rayleigh’s radiation alone is considered, we have 






0 ^ 


^^{TV) + lc^rV 


_ ih djrV) „ _ 1 9^ (r7) 

* rsin^ 3^ ’ * vr drdd 

„ ik d(j-V) p 1 d^jrV) 

^ r ^ VI' sin 6 drd<l> 

■where 


...(103), 


W -7 


sin^sin ^...(104). 


At a great distance from the origin the radial electric force 
is of order ^ while the transverse electric and magnetic forces 

are of order Hence the intensity of the scattered light 

diminishes ultimately according to the inverse square law when 
points on the same radius are considered. It also varies as the 
square of the volume of the particle. 


* TMs remark is made by both Maxwell Garnett and Mie. 
t It had previously been considered by Lord Rayleigh, “ On the influence of 
obstacles arranged in rectangular order on the properties of a medium,’’ PML 
Mag. (5), Vol. 34 (1892), p. 481 ; Scientijic Papers^ Vol. 3, p. 19 ; and by 
T. H. Havelock, Proc. Boy. Soc. A, Vol. 77 (1906), p. 170. 

J Laue, Eriedrich, und Knipping, Sitzungsher. der Konigl. Bayerischen 
Akad. d. Wiss. Juhe 1912. 

§ See, for instance, W. L. Bragg, Proc. Camb. Phil. Soc. Vol. 17 (1913), 
p. 43 ; Proc. Boy. Soc. A, Vol. 88 (1913), p. 428 ; M. Laue, Milnchener Ber 
(1912), p. 363; Ann. d. Phys. Bd. 41 (1913), p. 989, Bd. 42 (1913), p. 397*; 
P. P. Ewald, Phyt. Zeitschr. (1913), p. 465 ; L. S. Ornstein, Amsterdam Proc. 
(1913) ; M. Bom u. T. v. Karman, Phys. Zeitschr. (1912), p. 297. 
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An approximate formula for the intensity is 

I = (eos^ 6 sin^ (/> -f cos^ <jE)) . . .(105), 

the intensity of the incident radiation being 

If cr=:0 for the medium outside the sphere, the quantity 
k is inversely proportional to the wave-length X of the incident 
radiation. Hence when h is large compared with k, we have 
Lord Rayleigh's result that the intensity of the scattered light 
varies inversely as the fourth power of the wave-length. The 
short waves are on this account scattered far more profusely 
than the long ones and so we have an explanation of the blue 
colour of the sky. 

The above formula for the intensity of Rayleigh's radiation 
indicates that there is no light of this type in a direction for 

which 0 = = i.e. when the observer is looking in a direction 

parallel to the direction of the electric vibration in the incident 
wave. To obtain an expression for the intensity of the light 
sent out in this direction we must take into account the second 
electric wave and the first magnetic wave. Referring back to 
the expressions for and Ai, we find that when cr is neglected 
both inside and outside the sphere, the intensity of the scattered 
radiation varies inversely as the eighth power of the wave- 
length*. This corresponds to Tyndall's residual blue" which 
is purer than the blue seen under other conditions. 

§ 22. The absorption of light by a spherical obstacle. 

The total energy absorbed from the incident radiation by 
a single particle consists of two parts; first of all the energy 
scattered and secondly the energy which flows into the particle 
and is transformed into heat or chemical energy. Both these 
quantities may be calculated by the followung method which is 
a simplification of the one given originally by Mie. 

The flow of energy across unit area of a very large sphere 
concentric with the spherical particle takes place at a rate 
measured by the radial component of Poynting's vector, i.e. 

* Lord Rayleigh, mil. Mag. VoL 12 (1881), p. 81. 
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Now if E, E and H, H are conjugate complex quantities, 
this component is equal to 

Integrating with regard to t so as to obtain the mean value 
of this quantity over a period, we obtain an expression which 
may be written in the form f 

S = ~ {M,M^ - MeM^, - 

olv 

where M, M* are the values of M corresponding to the signs 
+, — respectively in (86) and AT, M* are derived from them by 
changing the sign of i. 

The function fl for the outer space is the sum of £1q and ili, 
hence we may write 

rfl — S Pn^ (cos 6) [une^'"^ + 

n^l 

_ 00 

rfl = 2 Pn^ (cos 6) [une~^ + Vne~^ + w„e^] 

n=l 

rll* = 2 Pn (cos 0) [Une’~^ + 'Vne^'^ + Wne'^] 

n=l 
00 

rfl* = 2 Pn^ (cos 6) [Une^^ + 4- 

«=1 

where 

(^n - Pn) (^)- 

When the expression S is integrated over the spherical 
surface, the result can be expressed in the form 

Iq’\' Ii Iqi) 

where Jq depends only on the incident radiation, Jj depends 
only on the scattered radiation and represents the amount of 
energy absorbed by scattering, Joi depends on both types of 
radiation and represents the total absorption of energy from 



t We assume now that (r=0 outside the sphere. 
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the incident wave. The sum of the three terms with its sign 
changed represents the amount of energy that flows into the 
sphere and is truly absorbed by it. 

Performing the integration with regard to (f> and collecting 
the different terms together, we find that the surface integral 
can be written in the form 


— OO 00 

t 

wi, = l 


/o 


+ sm= 5 


dPJ dFm^\ dd 


dd dO J sin 9 


r p 1 , dPm p i'\ yg 


— lew. 


dWfy^ 

'dr ] 

UdUn. dVn \ / du,n dVm \ _ d^ d^ 


[V dr dr j \ dr dr ] dr dr 

+ kk (Un + Vn) (Um + ^m) - kTcW^Wr, 

This expression can be simplified with the aid of the relations 


/:( 


1 (]p 1 

p 1 p I 

m i T/i ^ 


dti 


de 


^d9 = 0 


.(107), 


/:( 


Pn Pm + sin^ d ■ 


n m 


de dd 


de ^ 

/sin e 
271- (n + 1)- 
2^4- 1 


0 m^n 


m^n 


...(108). 


The first of these is evident since the integrand is the 
differential coefficient of a function which vanishes at both 
limits. To prove the second we make use of the formulae 

d ^ 

sin 6 (cos ^) = “” cos 0 (cos 6) 

del 

— n{n + 1) sin d . Pn (cos 9) 


Pn^ (cos 9) = --^Pn (cos 9) 
1 + cos® 9 


sin0 


PnPm^ + n (n + 1) cos 0 . PnPm 


+ m (?n + 1) cos 9 . PmPn 


.( 109 ). 


= - ^ [cos 9 . 
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The integral is thus equivalent to 

n(n+ l)m(m + l)J Pn(cos 0) Pm (cos 6) sin 6d6, 

and so has the value we have assigned to it*. 

Our surface integral now reduces to 

\ + 'kWn 


[CH. 


TT I, , .(dun , dv^\ . , dw, 

'>n + 1 3^ 


S \k(u. 


\dr'^ dr 


dr 


dv. 




— kWn 


dw^ 9 ^ 2 ( 7 ^+ ly 
dr 


Since cr = 0 outside the sphere, we have h=^h, 
r is very large we may use the approximations 

2^^ + l 


2?2< + 1 
Also when 


Un^Q>o^\kr 


Un'^oos \kr — n + l ■ 


2n + l 


y...(iio). 


’’l I.-2. 

2) ‘ n{n + l) 

t)„ ~ i + S„) 

Wn'^\ {An - Bn) 

Vn'»^{An + Bn) 6^^ (- f)”+l 

~ i (2n - P«) e*" (- 
We thus find that Iq = 0 and 

7„ = £ n(n + l) [i» ( J„ + 5„) - (- i)» + Bn)] . • .(HI), 

C12)- 

When v = l these expressions are just half those given by 
Mie on p. 436 of his memoir. The reason why they must be 
doubled is that our expressions for the incident light give an 
intensity 

In the case of a solution containing A” spherical particles 
per unit volume the total absorption coefficient is thus 

_ "KT 00 

A=^ t n(n + l)[f«(^„ + P„)-(-i)"(.l„ + P„)]...(113). 

n^l 

It is interesting to study the variation of this quantity with 
the wave-length for particles of different sizes. Mie has drawn 
* Cf. L. Lorenz, Oeuvres scientifiqueSf p. 526. 
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in] 

the absorption curves for particles of gold varying in size from 
20 fjbiJL to 180 fji/jL, 



It will be seen that for ruby red gold solutions containing 
very fine particles there is an absorption maximum in the 
green corresponding to a wave-length of about 525 ya/x. 

Mie has also drawn curves showing the pure absorption in 
colloidal gold solutions. 

The colours of silver particles in colloidal silver solutions 
have been discussed with the aid of the mathematical theory 
by E. Mtiller*. The particles of a silver solution show beautiful 
colour phenomena, all colours of the spectrum from the extreme 


Ann. d. Phys. Bd. 35 (1911). 
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blue to the extreme red being present* For other applications 
of the theory we must refer to the papers of Lord Eayleigh, 
Garnett and Mie and to Prof. Wood’s Physical Optics. 

The mathematical theory for a large number of particles 
has been developed further on certain simplifying assumptions 
by F. Hasenohrlf, A. SchusterJ, W. H. Jackson§, L. V. King||, 
A. Einstein IT and M. v. Smoluchowski**. It has been used 
recently by W. J. Humphreys ff in a study of the effect on 
climate of large quantities of volcanic dust in the upper 
atmosphere. 


§ 23. The pressure of radiation on a spherical obstacle. 

We shall now calculate the pressure of radiation on a spherical 
obstacle, following the work of Debye except in some of the 
details. The pressure is calculated on the assumption that the 
force exerted by an electromagnetic field on a unit charge 
moving with velocity v has components of typej| 

= + (114), 

c c 


and that the equations of the field can be derived by a process 
of averaging from the electron equations 


Now if 


^ ^ 1 (dE \ 

div E — p, 


rot ^ — 


1^ 

c dt 


div = 0 


(115). 


X, = J - E/ - Ef) + i H/) 

iKy == ExEy + HxSy == Yx 

Xz = ExEz + HxHz == Zx 

Sx = c{EyHz-EzHy) 

* A coloured reproduction of an ultramicroscopic picture of a silver solution 
is given by H. Siedentopf, Ber. d. JDeutsch. Phys. Ges. (1910), p. 6. 

t Wim. Berichte (1902), % Astrophysical Journal^ Vol. 21 (1908). 

§ Bull, of the Amer. Math. Soc. Yol. 16 (1910), p. 473. 

11 Phil. Tra7i$. A, Yol. 212 (1913), p. 375. 

IT Ann. d. Phys. Bd. 33 (1910), p. 1275. 

** Boltzmann Festschrift (1904), p. 626 ; Ann. d. Phys. Bd. 25 (1908), p. 205 ; 
Phil. Mag. Yol. 23 (1912). 

tt Bull, of the Mount Weather Observatory (1913) ; Jourruil of the Franklin 
Institute, Aug. (1913). 

tt We assume now that a-=0, €=/t=l for the external medium, so that i^ = l. 
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etc., we have identically on account of (115) 

ia;S. 


P^X 


^Jx dXy dX, 

dx dy dz dt 


.(117). 


Hence by Green’s theorem, if (I, m, n) are the direction 
cosines of the outward drawn normal to a surface cr enclosing 
all the charges in the field, 

19^ 

dt 

The i:c-component of the force exerted by the electromagnetic 
field on the obstacle is thus the same as if there were tractions 
X, F, Z at each point of the surface and a volume force 
idS 


JI(lX^ + mXy + nX,)d<r^jll + dxdydz...(118). 


dt ‘ 


similarly for the other components. 


Now when we integrate with regard to t so as to obtain the 

9/5. 


mean value of the pressure over a period, the term 


dt 


contributes nothing on account of the periodicity of S. The 
pressure may consequently be calculated from Maxwell's tractions 
X, Yy Z. In the present case the mean value of the pressure 
in the direction of the axis of z is derived from the quantity* 
(IZ^ Hh wiZy -f- TiZ^ = — cos d \_E/q^ -j- 4" 

^ e;^ ^ h;^'\ ^ sin e ie;e' + E'E'\ 
The surface integral of the mean value of this quantity over 
a period is accordingly 

- ? rr [cos ^ {EtEe+E^E^+HeBe+H^B4.-ErEr-H,.Hr) 
^JoJo _ _ _ 

+ sin 6 (ErEe + ErEg + ErEg + ErE^)] sin 6 . dd dcp. 

When the incident radiation only is taken into account, this 

integral becomes simply 

<v,2 frr r2rr 

— ” I ( cos ^ sin 0 d6d^ 

^ J 0 Jo 

and vanishes completely. Again, it is easy to see that when r 
is large the components Er, Er for the scattered field are small 

* We assume now that the surface cr is a sphere whose centre is at the origin. 
Our expression is easily obtained by writing down the expressions for j&y, Ey, E/ 
in terms of E/, Ed% E<p\ 
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compared with the tranavoi-se components, a product such as 
Hr3g being of order may be neglected when r is large, au<l 
so we have only to consider the integral 

riT fitr 

Pjj = — j I I (Effli/Q 4* IIqMq + H^l[^ sin 6 cos 0 d$d<pj 

4 J 0 j 0 

where JS= JilQ + JHu + The torins which depend 

only on the incident field may, moreover, be diarcgardcd. 

To evaluate we necul the values of the integrals 

/WV dI\A 




+ sin ’ d 


(le (W 


-fdPn' 

V dd 


eos0. 


cot 6 . ddf 
(id. 


By using the relations (109) we may transform the first of 
these into 

n (n + 1) m (m + 1) | i\ (cos 6) (cos 6) sin $ cos d dO 
Jo 

~ fo dd ^ ^ ^ 

Now 

(2n + 1) coa d . Pn (coa d) = (n + 1) i »+i (cos d) + wP„_i (coa d) 

( 119 ), 

and 

(2n + 1) cos d . P„’ (coa d) = wP>„,^i (coa d) + (« + 1) P*»_, (cos d) 

( 120 ); 

hence when the second integral ia integrated by parts wtj 
obtain two integrals of the typo (7(i) and ao wc have finally 

P = 0 m^n ± .1. , 

.(n-iyn(n + iy 
“ ^ (2n- 1) (2ft + 1 ) - 1 


nft='(ft+l)(ft + 2)* 

■^(2ft + l)(2ft + 3) + 


.( 121 ). 


When the second integral P is integrated by parts it becoinos 


n 

P ” P 9i^P m sin d dd “ 0 7ft at ft 

jo ^ 


ft(ra + l) 
" 2 ft +1 


m = n 


.( 122 ). 


Ill] 


PRESSUEE OF RADIATIOIT 


63 


Writing Pj in the form 

Pz-’-gJojo (MsMe + Me^Me* 

4- sin 0 cos 0d6d<^ 

and making use of equations (121), (122) and (106), we obtain 

p S 4- 1) (n + 2)^ r f dunj,i dvn+i \ f du^ dvn\ 

^ 01=1 (2ti + 1) ( 2^1 4- 3) dr dr J \ dr dr) 

dVJ*fiJi~-Y d'lDfl , 7 n / V /•— _ \ 7 n 

+ (Woi+I + y-ft+i) (Wn 4- v,») 4- k^Wn+iWn 


dUn 

dr 


dVfi 

dr 


^^n+l \ _ 1 _ dw,^ dw, 
dr dr ) 


dr dr 


+ Jd‘{Un + IJn) (Mn+l + ^n+i) + k-W^Wn+i 

, . " w(n+ 1) ' 

-Aaw 2 , , I 

„=1 2n, + l 


'n+i 


-[ \( dU<yi 


•(^fcw4■ v«)(^ 


dr dr 


(^i7i4-^ox)(^' 
dvn 


4" Wf, 


( dUn ^ dVn 
dr dr 


dWn 


dr 


■) 


We now use the asymptotic expressions for Un, when 

r is large and omit the terms that depend only on the u*s. 

We also write 


jLyi 


*71—1 


2n 4- 1 




n(n+l) 

and obtain after some simplifications 




{n-i + ^ Q 


Pz==— "n 2 (2w + 1) [On, + «» + /3n. + ^»] 

^n-1 

00 2^ 4-1 — 

+ Ji jr(jr+i) 

00 7jr72,-j-2) _ — 

-j- Tt/c^ ^ \P-n^n+i + ^n^n+i + + Bn^n+i\ 

n=l ^ 4- 

(123). 


This expression turns out to be negative, i.e. the pressure 
acts in the direction in which the incident light is moving. If 
the constants in the incident light are chosen so that its 
intensity is unity, the above expression must be doubled. The 
numerical value of the pressure has been calculated from the 
above formula by Debye in a number of cases. When the 
radius of the spherical obstacle is small compared with the 
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wave-lengtli of the incident light, the functions yjrn {ka), {ka), 

'\}rn{ha), ^n{ka) occurring in the expressions for can be 

expanded in ascending powers of a. This method, however, 
fails when ka approaches unity and numerical values of the 
functions must be used. 

In the case of a totally reflecting sphere 


„_ i^n(ka) 

(ka) ’ 


'’"Im 


and Debye finds that if 


)=:ak=^ 


27ra 


where X is the wave-length of the incident light, L denotes the 
light-pressure and W — ^Tra^ is the energy of the incident train 
of waves per unit length of a cylinder circumscribing the 
spherical obstacle and having its axis parallel to the direction 
of motion of the waves ; then 

— = .^p4[i (125). 


The first term of the series was given by Schwarzschild. 
The convergence of the series is slow as p approaches unity and 
several terms of the series (123) must be taken into account. 
By using numerical values of the functions ^n{p\ ^n(p) aD.d 
their derivatives for n = l, 2... 5, Debye has succeeded in 

drawing a curve for 
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It will be seen that the pressure has a maximum value for 
a certain value of p, approximately equal to 1. When p is 

large the ratio approaches asymptotically the value 1. 


Debye compares the light pressure so obtained with the 
gravitational attraction for a spherical particle of specific gravity 
s under the influence of the sun’s radiation. He finds that if Q 
is the gravitational attraction, 


L_mo L 

G ■“ \sp * F’ 

and to get a numerical estimate he takes X = 600 yayu, s = l. 

It appears that the ratio vanishes both for small and large 
values of p : it has a maximum value of about 20 for p = 1. 

In the case of a dielectric sphere with refractive index n, 
the expansion corresponding to (125) is 

+ 15' (n.2 + 2)(27i“ + 3) 

( 126 ), 

and is suitable for calculations only when np is small. 

Debye has drawn curves for ^ in the cases n = oo , n = 2, 

n = 1'6 and n = 1-33. When n = 2 the curve appears to have 
three maxima and two minima between the values p = 1 and 



p = 3. 

The greatest value of 


Z 

w 


is now about 2*6; the following 


table indicates when the light pressure exceeds the gravitational 
attraction, the numbers po and pi give the extreme values of p 
belonging to the range in which this is the case. 


n 



Po 


Pi 


00 20 

2 13 

1*5 3 


*3 8 

*6 5 

•8 — 


The maximum light pressure is just balanced by the 
gravitational action when n is about equal to 1*33, the value 
for water : for smaller values of n gravitation prevails. 


B. 


5 
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In the case of an absorbing spherical particle, the equation 
which takes the place of (126) is 



(/c 


cr 

Q) 


+ 2)2 + 



(127). 


When a is small the light pressure and gravitational action 
are both of order and their ratio tends to a finite limit, hence 
for certain types of absorbing material there is no lower limit 
in the size of a particle below which gravitation exceeds the 

light pressure. Debye has drawn a curve for ~ for the case of 

a gold particle and finds that there is a maximum value for 
p = 1*5 nearly. 

L 

The existence of a maximum value for ^ in the cases that 

W 

have been discussed appears to be due to the fact that the 
value of p for which the maximum occurs is very nearly equal 
to the real part of the complex value of p corresponding to one 
of the free damped vibrations*. The first electric vibration 
seems to be of chief importance in determining the position of 
the maximum. 

The determination of the limiting value of for very small 

wave-lengths, i.e. for large values of p, is a matter of some 
difficulty, it depends on some expressions giving the behaviour 
of the Bessel functions for large values of n and p. These have 
been found by J. W. Nicholsonf and P. Debye 
If p is real and n + ^ < p, we have when p oo 


^n(p) = 


(sin t„)4 




co8(p/o-^) 
(sin To)^ 


( 128 ). 


* Cf. Debye, loc. cit.j and the similar remarks for tbe case of optical 
resonance by E, Pockels, Physik, Zeitschr. Bd. 5 (1904), p. 152. 

t British Association Re^ports, Dublin, 1908, p. 595 ; Phil. Mag. YoL 13 (1906) ,• 
p. 195 ; Vol. 14 (1907), p. 697 ; Yol. 16 (1908), p. 271 ; YoL 18 (1909), p. 6. 

X Math. Ann. Bd. 67 (1909), p. 535. 
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where to is an angle lying between 0 and -j for which 

n + ^ 


COS To = - 


P 

and /o = sin to — tq cos tq . 

When n + ^ > p and p oo , we have 


(p) = i 


■fn(p)-- 


-ipA \ 


(i sin To)^ 
Jpfo 


(i sin To)^ 

where To is now the root of the equation 


(129), 


COS To = 


P 


whose imaginary part has a negative sign. 

When n and s are very nearly equal the values of (p) and 
(p) can be made to depend on Airy's integral and are much 
more complicated ; for these we must refer the reader to the 
original memoirs. 


24. Other problems which may be treated with the 
aid of polar coordinates. 

The diffraction of electric waves travelling round the earth 
is a problem of some importance which has been discussed by 
H. M. Macdonald*, Lord Rayleigh f, H. Poincard;!:, J. W. 
Nicholson § and other writers. 

The calculations are very long and depend on the use of the 
formulae to which we have just referred. Rybcynski|| has 
recently treated the problem by a method due to March and 
has taken into account the finite conductivity of the earth. As 
we have already mentioned this was done by Zenneck and 
Sommerfeld for the case in which the earth's surface is treated 


* Proc, Poy, Soc. Vol 71 (1903), p. 251 ; VoL 72 (1904), p. 59 ; Yol. 90 (1914), 
p. 50; Phil. Trans. A, Vol. 210 (1909), p. 113. 
t Proc. Roy. Soc. Vol. 72 (1904), p. 40. 

J Rend. Palermo (1910) ; Proc. Roy. Soc. Vol. 72 (1904), p. 42. 

• § Phil. Mag. Vol. 19 (1910), pp. 276, 435, 516, 757 ; Vol. 20 (1911), p. 157 ; 

Vol. 21 (1911), pp. 62, 281 ; Jahrb. d. draht. Telegr. Bd. 4 (1910), p. 20. 

11 Ann. d. Phys. Bd. 41 (1913). 
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as a plane. The results of Nicholson and Poincar^ indicate 
that diffraction round a perfectly conducting surface is not 
sufficient to explain the apparent bending of the electric waves 
round the earth’s surface. A generally accepted opinion is that 
the ionisation of the atmosphere by the sun’s rays is a very 
important factor in producing the observed effects*. 

The diffraction of a solitary wave or pulse by a spherical 
obstacle might be discussed with advantage. The evaluation 
of certain definite integrals involving Bessel functions, however, 
presents some formidable difficulties which probably account 
for the fact that the problem does not appear to have been 
solved. 

The scattering of electric waves by a perfectly conducting 
conical obstacle has been treated very briefly by H. S. Carslawf. 


EXAMPLES. 

1. Prove that 

1+2 Pn (cos $) Fn (cOS a) = log ^Sec^ ^ COSeC^ 0 , 0 < a < ^ ^ TT. 

(C. Neumann.) 

2. . Prove that 

“ Jo (h sin a) = 2 (271+1) A (cos S) P« (cos a). 

(E. W. Hobson.) 

* Cf. the discussion at the British Association meeting, Dundee (1912), and 
an article by W. H. Eccles in the Year Book of Wireless Telegraphy (1913). 
Some quantitative experiments on long distance telegraphy have been made 
recently by L. W. Austin, who obtains an empirical relation between the 
magnitude of the current received and the distance between the two stations, 
Bvlletin of the Bureau of Standards, Vol. 7 (1911). 
t PhiL Mag. Vol. 20 (1910), p. 690. 
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§ 25. The wave-equation in Cylindrical Coordinates. 

If we put X p Gos <f), y = p sin <j>, the wave-equation becomes 


d^u Idu 1 d^u d“U 1 d^u 
0p^ p dp p^ dz- c- dt^ 


(130). 


Two particular solutions of this equation are suggested at 
once by the general solution of § 5 : they are* 



r2ir 

u=^\ F 
Jo 

z + ip cos a, 

5^ — ^ sin a J 

dec .. 

....(131), 

and 

r+<f> 

F 

J 

z + ip cos a. 

t — - sin a 
c 

doc .. 

....(132), 


respectively. The first of these represents a wave-function 
which is symmetrical round the axis of z and which reduces to 
27rF (z, t) when p == 0. It gives us at once the formulae 


1 

r”Pn — (^ + V cos a)” da. 

II 

...(133), 

1 p / N_ 1 T" da 

^n +1 ^ W 27rJ 0 (z-h ip cos 

II 

...(134), 


where n is an integer f, and many other interesting formulae 
may be written down by simply choosing different wave- 
functions that are symmetrical round the axis. 


. * The first of these is an obvious generalisation of a formula given by 
D. Edwardes, Educational TimeSy Oct. (1904). 

f The formula is also true under certain limitations when n is not an integer. 
See Hobson, Phil, Tracis. A, Vol. 187 (1896). 
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For instance, 

C”- Iog(£r + ^•/)cos 0 ) 


-r 

TT jo 


1 , 2r*- 

i _ lo^ 

T ^ r z 


.(135), 


+ ip COS 6 

1 log (z + ip cos 6) — log Zq + Zq-^R 

TT Jo z-Zo + ip COS 6 R * ?' Hh ^0 — -K 

(R^ = p^^-h(z^Zo)\ z>z,y,.(lS6), 


27r 


sin k(z + ip cos a) ikpsina i _ sinfer 


e ...(137). 


z-\-ip cos a r 

There is another general formula for a wave-function sym- 
metrical about an axis, viz. 


10 = 


r p 

= Ft--- cosh OL, z — p sinh a 

j — 00 ^ 


da ...(138), 


where the function F is of such a nature that 


1 . , dF ^ . dF 

- smh a ^ + cosh a 
c dt oz 


vanishes when a = + oo . As a particular instance of this we 
have the function* 

F(t-^ cosh a) tZa (139), 

which may be regarded as the cylindrical wave-potential for a 
line source of strength F (t) along the axis of z. 

A peculiarityf of the two-dimensional propagation of waves 
is the existence of a tail ” to the disturbance when F (t) is zerc 

for i5 < 0 and ^ > r, for if ^ > t + - we have 

c 

— ^ cosh da (140), 

where pcoshs = c(t—T). It is clear that this expression fo 
n does not generally vanish. The wave-function (139) is thus 
essentially different from Euler's wave-potential for a poini 
source, viz. 



* Cf. H. Lamb, Hydrodynamics, pp. 281, 500 ; V. Volterra, Acta Math 
t. 18 ; Levi-Civit^, Nuovo Gimento (4), t. 6 (1897). The formula is a particula; 
case of a more general one given by Dr Hobson in 1891. 

+ This was discussed by 0. Heaviside, PhiL Mag. (5), t. 26 (1888); Elec 
trical Papers, Vol. 2. See also Lamb, Hydrodynamics, p. 282. 
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lor in this case fl is zero for ^ > r + — , provided the source is 
only active when 0<t<r, 

V. Volterra* has obtained a number of elementary wave- 
functions of the form 

u = t”'F(^'^ = t”‘F(s) (141). 

He finds that F must satisfy the differential equation 
(1 - ^ (2^ - ^ ^ - 1) = 0 . . .(142). 


If we try to solve (130) by means of a function of the form 


u = Wp'^ cos m ((j> — (po), 

where W is independent of <f>, we find that W must satisfy the 
equation 


d^W 2m + l dW 
dp^ p dp dz^ & df 


0 . 


.(143). 


Solutions of this which are independent of z may be derived 
from the following formulaef, in which m > — I-, 


W = J — ^ cos a^sin^oc.da (l^^), 

W=J F (t±^ cosh 7}^ sinh^”^ tj ,di] (145). 


There are, of course, certain limitations concerning the 
behaviour of F (t) at infinity. The first formula enables us to 
determine the value of W when its value is known for p = 0. 


§ 26. Elementary solutions of Au 4* k-u = 0. 
The differential equation 


d^u 1 aw 1 d‘^u d^u ,2 

dp^ p dp p^ d<j)^ dz^ ^ 

possesses elementary solutions of the types 
Jm ip e cos m (^ 

Km ip — ¥) cos m (</> 


■ ^o) 

-4>o) 


.(146) 

< 147 ), 

(148), 


* Acta Math. Vol. 18. 

t E. W. Hobson, Proc. London Math. Soc. Vol. 22 (1891), p. 431. The 
first solution is due to Poisson. 
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where \ h and m are real or complex arbitrary constants. The 
first solution may of course be generalised into 

f (p^/¥+h^) e**V(A ) hdk (149), 

J a 


and a similar remark applies to the second. If we wish to 
express a given wave-function in the first of these forms, the 
following inversion formula due to Hankel is particularly 
useful* 


If 

then 


F (a:) = f (xt) f (t) t dt 

Jo 

/ (i) = f (^0 ^ (^) 


(150). 


Let us use this formula to express ~ in the form (149) 

when m = 0. Since the representation should be valid for ^ = 0, 
we find on putting A® 4, ^2 _ ^2^ a=^k, that 

ikp fao 

Jo(Xp)f(h)Xd\; 

p Jo 

I‘>0- 


Hence we obtain Sommerfeld’s formulaf 


1 

r 





XdX 


(151), 


the upper or lower sign being taken according as ^ 5 0. A 
more complete proof is obtained by applying HankeFs inversion 
formula to the equation 







whicli is established in Prof. Lamb’s paper. 


X^>I(P 

X=<ifc“ 


.(152), 


* See Gray and Mathews, Treatise on Bessel Functions (1895), p. 80 ; 

N. Nielsen, Handbuch der Theorie der Cylinderfunktionen (1904), p. 366. 

t Ann, d, Physik, Bd. 28 (1909), p. 683. Some analogous formulae are 
given by H. Lamb, Proc. London Math. Soc. Ser. 2, Vol. 7 (1909), p. 140 ; 

O, Heaviside, Electrical Papers, Vol. 2, p. 478 ; N. Sonin, Math.' Ann. Bd. 16. 
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A few more formulae will now be written down to illustrate 
the method of generalisation by integration with regard to 
a variable parameter 

, pikR 

Kq (p V — P) cos X b) — ij cos X(a-~b) da 
If m is zero or a positive integer and n is a positive integer, 

,^1 Pi (<=os e) = r(>,-i+ - i ) fo ^ ^ 

JflQl(^os 6) = l) /o z>0. 

(Hobson.) 

I 27. The propagation of electric waves on a semi- 
infinite solid bounded by a plane surface^. 

In this problem the surface of the earth is regarded as an 
infinite plane and the waves are supposed to be generated by an 
antenna, of which one portion is vertical and the other horizontal. 

Let us assume that the electric and magnetic forces are the 
real parts of vectors of the form respectively, 

then if if = jy ± ivE, where v- = , we may satisfy Max- 

fL€d 

well’s equations by putting 

if = rot n ± ~ grad div II ± AjII (153), 

where AH 4- k^U = 0, e/uLco^ -h i/icoa (154). 

To imitate the action of the antenna we shall place two 
vibrating doublets at a point at distance a from the plane. If 
one of these vibrates vertically and the other horizontally, we 
may put for the primary radiation IIo = (Pa;, 0, Pg), where 

P. = -B^, Pz==A^, m = p‘ + {0-ay...{155), 

and the axis of ^ is vertical. We write JB with a negative sign 
to signify that the horizontal branch of the antenna is drawn 
in the negative oj-direction. 

^ A. Sommerfeld, Ann. d. Rhys. Bd. 28 (1909), p. 665 ; H. v. Hoerschlemann, 
Jahrb, d. draht. Teleg. Bd. 5 (1912), pp. 14, 188; Dissertation, Munich (1911). 
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A convenient expression for IIo is obtained by using Sommer- 
feld’s equation 


R 




Jo 


XdX 

~r^ 


0<z<a. 


.( 156 ), 


where l=^'^X^—kK Appropriate functions ITi, Ila for the re- 
flected and transmitted disturbances are obtained by putting 


oo 

ni= S GOSTKp I Jn(Xp)e~^^^’^^^ Fn(X)dX 
n=0 Jo 


/•oo 

1 I t/jj, (Xp^ 

Jo 


^m(z-a) a^(\)dX 


z<a.., (157), 


112 = S cos 91<3 

n=^0 

where — and h is the value of h in the second 

medium. The functions (X), Qn (X) are vectors with com- 
ponents [/n(A^), 0, 0, %n(X')] respectively. 

We can satisfy the condition that the tangential components 
of the electric and magnetic forces should be continuous at the 
surface of separation of the two media by putting for ^ = 0 

Hq -h Hj = 1 I 2 
rot ( Oo 4* Hi) = rot Da 


^ div(no+no=-^divn, 




.( 158 ), 


where fi has been taken to be unity for both media. 

Substituting the integral expressions for Ho, Hi, Ha in these 
equations and equating to 2 ero the coefficients of functions 
of type Jn(Xp) in the resulting integral equation, we obtain 
the system of equations 


T) 

/o e~^ — j- e~’^ = gn e~”“, 
+ -j = Xo 

Z/o - Be~^ = rnffo 


■fne~^=Xn^~^> 
lfne-^ = mgner^, 


m 


{Ae e" 


e la ^ ^ ^ - Xg.e-^), 
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where the last equation has been simplified with the aid of the 
relations fn—gn = 0(n>0) which are a consequence of the 
previous equations. 

Solving these equations we eventually find that if 
= (QiB, 0, Qz)} Il2 = (-Ra;j 0, 


Qz — A j ^J q(Xp) 


— k^mXdX 

g— Us+a) 

hH + Ic^m I 


+ 2SCOS 

R..Uk- fMX, 

+ 2£ »s ^ j; (Xp) 


The directed effect ” depends on the presence of the terms 
involving cos (f> in the expressions for Qz and J?z. Now when 
cr = 00 for the second medium, h = oo , and these terms vanish 
altogether ; hence the possibility of directing the energy of the 
radiation sent out from the bent antenna is due to the im- 
perfect conductivity of the earth. Von Hoerschlemann has 
given a numerical discussion of the above formulae but the 
investigation is too long to be inserted here 


§ 28. Propagation of electromagnetic waves along a 
straight wire of circular cross-section*. 

Let us consider the symmetrical case when the electric 
force at any point is in a plane through the axis of the wire 


* H. Hertz, Electric Waves ; J. J, Thomson, Proc. London Math. Soc. 
Vol. 17 (1886), p. 310; Becent Besewrclies, § 259; A. Sommerfeld, Ann. d. 
JPhysiJc, Bd. 67 (1899), p. 233 ; Gray and Mathews, Bessel Functions^ Ch. 13; 
M. Abraham, Encyhl. d. Math. Wise. Band V. 2, Heft 3 (1910), p. 526; 
J. Larmor, Proc. of the 5th Int. Congress of Mathematicians., Vol. 1 (1912), p. 206. 
The problem eonsidered in this last paper is chiefly that of alternating currents, 
viz. the forced alternations of flow produced by a uniform periodic electric force. 
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0? 


(159). 


and the magnetic force is in circles at right angles to this 
plane. The field equations are then of the types 

€ O’ « 10^ jj- \ € 0Jip . ^ TVT __ 9 ) 

_^9^*^0/fp 
c dt dz dp 

These mjiy be aatisfiet.! by putting 

1 a / ari\ 








vdpJ^ 

^dm an 

c dtdp c dp 

where H satisfies the equation 

a^n / am 


a^^iix 
dpdz I 


€pLdm , cTfjLdu 

' + 


c* dt^ ^ dt ^ dz^ p dp V / 


:AII 


.( 160 ), 


.( 161 ). 


Putting n = u, we find that Au + T(^u = 0, where 


, ^ — ifioxr 


We now assume that for points outside the wire* 

u = A /vo (p Vx® -- &*), 

and that for points inside the wire 

=s Jq (ip — A®), 

where A is the value of ^ inside the wire. These asBumptions 
are made for the purpose of determining the periods and rates 
of decay of electric waves that can travel along the wire and 
maintain their own field. The first solutitm is chosen so m 
to make the flow of energy negligible at an infinite distance 
from the wire so that th(‘. system is self-contained ; and to 
ensure this it is nccoasary to suppose that the real part of 
Vx® — A® is positive. The second solution is chosen s<^ as to 
make the electric and magnetic forces finite on the axis of 
the wire (p =» 0). 

Let p « a be the equation of the surface of the wirc^, then 

* We follow bore tlie work of Bommerfekl in which ifc is suppoaod that diaro 
are no oondnetorB outside the wire. Bir J. J. Thomson allows for the prosenoo 
of external conductors by supposing the dielectric surrounding the wire to be 
bounded by a cylindri cal’ conductor having the same axis as the wire. 
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the continuity of the tangential components of the electric and 
magnetic forces requires that when p = a 


A 


dp 


pKo(p'\/x^ 


■/c») 


= B 


dp 




I 




(162), 


^dp ' J 

V and Vi being the values of the quantity cr + icoe at points 
outside and inside the wire. The elimination of A and B gives 
rise to a transcendental equation for the determination of X. 
The total current flowing along the wire is 


t/ = 27r [ crEzpdp 

Jo 


= — ^irae^^pB ^jQ{ip^X^— h?) for p = a. 

On the other hand the electric force at the surface 
of the wire is 


p dp 


pBe ^0 (ip V^X^ — h^) 

= - jB (\® - h^) e-^ Jo (ia Vx"-A=>). 
Hence, if we put 


for p = a 


E, = J 




where JR and L are the resistance and self-induction of the wire, 
we get 


27rcr 


a?=wpVA*-^® 


The roots of the transcendental equation have been discussed by 
Sommerfeld who used a method of successive approximations. 
The values of X are complex as the waves which travel along 
the wire are damped owing to the imperfect conductance of the 
wire. It appears that when the disturbance does not penetrate 
far into the wire, the damping is small and so the velocity of 
propagation is very nearly equal to the velocity of light. 
When, however, the field does soak into the wire to some 
extent, the damping is of course considerable and so the wave 
travels with a velocity a li ttle le ss than that of light. In the 
first case the real part of Vx^ — is large, in the second case it 
is small. 
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In Lecher’s arrangement^ there are two conjugate parallel 
wires between which the waves travel, consequently the field is 
not symmetrical round the axis of one of the wires. This case 
has been discussed by G. Mief with the aid of bi-polar 
coordinates, and the lines of force have been studied by W. B. 
Morton The latter also considers the case of n parallel wires 
passing through the comers of a regular polygon §. 

The mathematical analysis for the case of a curved or twisted 
wire has not yet been fully developed. The important case of 
a spiral wire has, however, been discussed by H. C. Pock- 
lingtonjl and J. W. Nicholson IT. The latter gives numerous 
references to the literature of the subject. D. Hondros* * * § ** * * §§ ^ has 
recently discussed the propagation of some types of unsym- 
metrical waves along a single wire. The electromagnetic 
theory of an electric cable has been given by Sir Joseph 
Thomsonff and F. HarmsJJ. The latter considers the case 
when the outer conductor of a cable is replaced by air. 

§ 29. Other problems which may be treated with 
cylindrical coordinates. 

The diffraction of electromagnetic waves by a cylindrical 
obstacle has been discussed by Lord Eayleigh§§, W. Seitz||||, 
W. IgnatowskylTU, P. Debye^f, C. Schaefer^J, and J. W. Nichol- 
son*§. Schaefer has made an extensive study of the case of a 

* Ann. d. Phys. Bd. 41 (1890), p. 850. See also D. Mazzotto, II Nuovo 
Cimento (4), t. 6 (1897), p. 172. 

t A?in. d. Phys. Bd. 2 (1900), p. 201. The case in which the capacity 
is small is discussed at length by J. W. Nicholson, Phil. Mag. Feb.— Sept. (1909). 
The current is supposed to flow along one wire and return along the other. 

J Phil. Mag. Vol. 50 (1900), p. 605; Vol. 4 (1902), p. 302. 

§ Phil. Mag. Yol. 1 (1901), p. 563. 

11 Proc. Gamh. Phil. Soc. Vol. 9 (1897), p. 324. 

H Phil. Mag. Vol. 19 (1910), p. 77. 

** Alin. d. Phys. Bd. 30 (1909), p. 905; Dissertation, Munich (1909). 

tt Droc. Boy. Soc. Vol. 46 (1889), p. 1 ; Recent Researches, p. 262. 

Xt ^nn. d. Phys. Bd. 23 (1907), p. 44. 

§§ Phil. Mag. Vol. 12 (1881), p. 81. 

1111 Ann. d. Phys. Bd. 16 (1905), p. 746. 

ITIF Ann. d. Phys. Bd. 18 (1906), p. 495. 

Phys. Zdtschr. (1908), p. 775. 

Ann. d. Phys. Bd. 31 (1910), p. 462. 

*§ Proc. London Math. Soc. Ser. 2, Vol. 11, p. 104. 
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dielectric cylinder and his results have been tested experimentally 
by Grossmann. This is one of the few cases in which the 
influence of the material properties of the obstacle has been 
taken into account in the mathematical treatment of a diffrac- 
tion problem;, the necessity of doing this has been clearly 
indicated by some of the experimental results. 

P. Debye has studied the diffraction problem with reference 
to the theory of the rainbow. This had been done for the case 
of the sphere by L. Lorenz long ago^. 

Electrical vibrations in regions bounded by cylinders of 
various shapes have been studied by Sir J. J. Thomson-f*, Lord 
KayleighJ, Sir J. Larmor§, R. H. Weber |j, A. Kalahnef, and 
J. W. Nicholson*^. The latter has also calculated the pressure 
exerted by a train of plane electromagnetic waves on a perfectly 
conducting cylinderff. 


EXAMPLES. 


1. An infinitely long metal cylinder of specific conductivity <r 
and permeability /x, bounded by the surface is surrounded by a 

dielectric of specific inductive capacity e. A train of waves, in which the 
electric force is perpendicular to the cylinder and to magnetic force and 
if undisturbed would be represented by the real part of + 

is passing in the dielectric. Prove that the magnetic force J3‘2 inside 
the cylinder and the part Hi of the magnetic force outside representing 
the scattered wave, are given by the real parts of 

2 cos mcj>, (hr) cos m(p, 

m in 


where A2= —4:7riJLai(o, and 


47rcr 








* Oeuvres scientij^ques, pp. 405 — 502. See also G-ans and Happel, loc. cit, 
(p. 44). The most recent paper on the rainbow is by W. Mobius, Arm. d. Rhys. 
Bd. 33 (1910). 

t Recent Researches , p. 344. 
t Phil. Mag. (5),' Vol. 43, p. 125. 

§ Proc. London Math. Soc. (1), Vol. 26, p. 119. 

11 Habilitationsschriftj Heidelberg (1902) ; Ann. d. Phys. (4), Bd. 8 (1902), 
p. 721. 

IT Ann. d. Phys. (4), Bd. 19 (1906), pp. 80, 879 ; Bd. 18 (1905), p. 92. 

** Phil. Mag. Aug. (1905), May (1906). 

ft Proc. London Math. Soc. Ser. 2, Vol. 11, p, 104. 
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for positive integral values of m. The constants here have reference to the 
electromagnetic system of units. 

(Cambr. Math. Tripos, Part II, 1905.) 

2. Plane electromagnetic waves represented by 

fall upon the perfectly conducting cylinder Prove that in the 

scattered field 

(P. Debye.) 

3. The wave-potential for a circular ring of point sources is given by 

Q=J“ 0>O. 

(A- G. Webster.) 

4. The wave-function 

a=e^^^J^e~^Jo (p VF+p) cos aX c^X 

is zero at points on the plane 2=0 which lie inside the circle p^<a^; for 
z =0, p^>a^ its value is (p ^ — cos {ks/ p^ ~ a^). (Sonin.) 

5. If Ii^=^p^+(z+a8mh.uf, the integral 

JVo {\p)-Kd\ j" Jo j^X ^ 

represents the function outside a paraboloid of revolution 

whose focus is at the singularity 2 = —a sinhw, /)=0, and which passes 
through the circle p=*a, 2=0. It is zero inside the paraboloid. 


6. The circuital relations in cylindrical coordinates are 


p dEp dllis dJET^ 


c dt dtp ^ 

p dEz ^ , rr s 


dz 
’ d(f> ’ 


p dB<i> 
c dt ' 


dE^ 


dE, 

■p-W’ 


dz 




=0 


and similar equations in which E is replaced by —if and H by E. 


7. If jr= ^ cos - y sin 6)^5 T=^os sin ©^-by cos ©^, Z=^ 2 - 'ot^ where v 
and © are constants, the function Q—F{X, Yy Z) is a wave-function if F 
satisfies the partial differential equation 

j3X2“^V J dY^’^'^1^^^ dXdY'^c^\^ dX'^ ^ dY) * 

2 \ 02^ 


27 - 
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Obtain particular solutions of the form 

where p^=X^ ( ^ m, A, B are arbitrary constants. 

8. An oscillatory current is induced on a circular wire of radius a 
excited by a uniform electric force acting on its surface from the 

surrounding medium. Obtain expressions for the inductance and resistance 
of the wire per unit length when the wire is regarded as straight and no 
disturbing conductor is near. 


B. 
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CHAPTEK V 


THE PBOBLEM OF DIFFRACTION 


§ 30. Multiform, solutions of the wave-equation*. 

The wave-functions required to solve many of the boundary 
problems of Mathematical Physics are not single- valued func- 
tions of X, y, z, t in an ordinary space. We may, however, 
regard them as single-valued functions in a Riemann's space. 
This is a simple generalisation of the Riemann’s surface of the 
theory of functions of a complex variable f ; every plane 
section of the Riemann’s space is in fact a Riemann’s surface. 
Instead of branch lines and branch points we have branch 
membranes and branch curves. Thus in the physical problem 
of the diffraction of light through a circular hole in a screen, 
the boundary of the shadow of the screen is the branch mem- 
brane and the edge of the hole the branch curve. 

We shall commence by finding a multiform solution of 
the equation 


d^u d^u 
doc^ dy^ 


-f k^u = 0 


,(163). 


The fundamental solution u = e** cos « + 2^ sin a) ^ ^ikp cos (<^ - a) jg 
of period 2-^ and can be expanded in the form 

^hp cos(^-a) _ J^Qcp) -p 2 2 i'^Jnikp) cos .(164). 

1 


* This theory is due to A. Sommerfeld, Math. Ann. Bd. 45 (1894), Bd. 47 
(1896) ; Zeitschr. filr Math. u. Phys. Bd. 46 (1901) ; Proc. London Math. Soc. 
(1), Yol. 28 (1897), p. 417. It has been developed by H. S. Carslaw, Proc. 
London Math. Soc. (1), Vol. 30, p. 121 : (2), Vol. 8, p. 365 ; Phil. Mag. Vol. 5 
(1903), p. 374; Vol. 20 (1910), p. 690 ; Fourier"* s Series and Integrals^ Ch. 18; 
W. Yoigt, Gott. Nachr. (1899) ; E. W. Hobson, Camb. Phil. Tram, Vol. 18 
(1900), p- 277. Different methods have been used by H. M. Macdonald, Electric 
Waves, Appendix D ; K. Schwarzschild, Math. Ann. Bd. 55 (1902), p. 177; 
Proc. London Math. Soc. Vol. 26 (1895), p. 156 ; C. W. Oseen, Arkiv for mat. 
Bd. 1 (1904), Bd. 2 (1905). 

t See Harkness and Morley’s Theory of Functiom (1893), Ch. 6- 
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A solution of period 2m7r may evidently be constructed by 
writing down a series of the form 

S anJn(kp) cos — (<jb — a), 

» = ra 

where the a^’s are suitable constants. The solution that seems 
the most natural extension of (164) is 

F^{p, cl>, 4>,) = Jo(/cp) + 2 21-^ cos ^(<f> - <^„) . ..(165). 

To sum this series when m = 2, we transform the terms for 
which n is odd by means of the equation 

n 

Jn (kp) sin" a . da, 

" 25r(i)r("-±l)^“ 

Summing the two series separately we find that 

F, (p, <^o) = e*" +/(<!> - <f>o) +/(<^>o - 4>), 


where 



hence we may write 

{p, 4>, i>o) = 


where t = ^2 kp cos |(<^> - <^o)- 

With the aid of the function F. we can solve some problems 
on the diffraction of plane electromagnetic waves by a semi- 
infinite plane bounded by a straight edge. Let us consider 

6—2 
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2 


the case of a totally reflecting screen*. If the electric force 
in the incident wave is parallel to the edge of the screen, the 
electric force u — Ez for the total disturbance must vanish over 
both faces of the screen and must satisfy the ditferential equa- 
tion (163). These conditions are fulfilled by taking 


Ez = F, (p, <!>, </>o) ^F,(p,<f>,^cf>,) (167). 

This value of E^ also satisfies the right conditions at infinity. 
To prove this we must find an asymptotic expression for F^ 
when p is large. 

Now when r > 0, we have the asymptotic expansion f 
ie" 


J r 


dX — • 


1 1-3 


while when r<0, f has a similar asymptotic expansion with the 

J -00 


sign changed. This means that when cos ^ (<p — <f>o)> e> 0 we 
have 






^—ikp , 


ZV 27rkp cos -I* (<^ — <^o) 

while when cos ^ (<^ — ^o) < ^ < 0 there is a similar asymptotic 
expansion in which the first term is missing. It thus appears 
that the electric force in the geometric shadow vanishes at 
infinity to the order 

If the magnetic force in the incident wave is parallel to the 
edge of the screen, the magnetic force ib = must satisfy the 


du 


differential equation (163) and be such that ^=0 over both 


faces of the screen. The conditions are fulfilled by putting 

Hz = F, {p, cf>, </>o) + F, (p, (^, ~ <!>,) (168). 

Prof. H. M. Macdonald has shown that the solution of 
a problem concerning a perfectly absorbing body can be made 
to depend on the solution of two allied problems^. “ A per- 


* The incident waves are supposed to conje in a direction for which 
^=7r + ^o- approximate solution of this problem was given by H. Poincar^, 
Acta Mathematical Bd. 16, p. 297 ; Bd. 20, p. 313. 
t Bromwich’s Infinite Series, p. 328. 

X Phil. Trans. A, Vol. 212 (1912), p. 337 ; Proc. London Math. Soc. (2), 
Vol. 12 (1913). 
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fectly absorbing body may be regarded as a body which is 
incapable of supporting either electric or magnetic force ; 
hence if G is the electric current distribution on the surface 
of the body when it is supposed to be perfectly conducting, 
and O' is the magnetic current distribution on the surface of 
the body when it is supposed to be incapable of supporting 
magnetic force, the superposition of these two distributions 
gives the electric and magnetic current distributions on the 
surface of the body when it is perfectly absorbing and the 
amplitude of the incident waves is doubled.'* 

Now if we suppose our screen to be incapable of supporting 
magnetic force, the boundary condition is that the tangential 
component of the magnetic force should vanish. When the 


electric force is parallel to the axis of z, 
the screen. Hence 


dE, 

defy 


must vanish over 


E, = F, (p, <jy, (^o) + F, (p, cjy, - (169). 


The solution for a totally absorbing screen is thus simply * 

E,^F,(p.<ly,<f>o) (170). 

Similarly, it can be shown that when the magnetic force 
is parallel to the axis of the solution for a perfectly 
absorbing screen is 

H,^F,(p,<ly,<f>o) (171). 

Prof. Lambt has discussed the case of perpendicular incidence 
with the aid of the parabolic substitution 

^ = p^cos^(j>, ri- sin 

The curves ^ = const., tj = const, are confocal parabolas, 
is the screen. A solution of Maxwell's equations is obtained 
by writing 



where 





11 's. 


dEx du 



du 

Q 

^2 = 0 


dt °dy' 

dt 

d^u dht 

9a!’ 

Id^u 


....(174). 


daf^ dy^ 

~c'^W 



W. Voigt, Gott. Nachr. (1899), p. 1, discusses the case of an absorbing 
screen. 

t Proc. London Math. Soc. (2), Vol. 4 (1907), p. 190 ; Vol. 8 (1910), p. 422. 
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Starting with Poisson’s wave-function* cos ^ ^f{ct — p), 

let us put 

du 

Transforming to the coordinates rj, we obtain 

= (175). 

Solving this partial differential equation by Lagrange’s method 
and adjusting the complementary function so that the boundary 

condition ^ ^ 9 ; = 0 is satisfied, we obtain 

+^F{ct+y)+ iF{ct- y). . .(176), 
where F is an arbitrary function. If the boundary condition is 
24- = 0 for 2 ; = 0 the sign of the second term must be changed. 
It is easy to verify that each of the integrals represents a wave- 
function. 


2 /'* 

Let us now put /(cc) = ~ j F' (x — v^) dv, 

and make the substitution f — cr cos a, 2 ; = <7 sin a, then after 
a little reduction we find that-f* 

TT 

1 r'2 

u =F(ct + y) F[ct-^y — (p-\-y) sec^ a] dec 

J 0 

w 

1 

+ -j^-F’[ct — y — (p — ^)aec^cc]da x<0, 

ir 

u = F (ct+'i/) + F (ct —y) — ^ j~ F [ct + y - (p + y) sec’‘ a] do 

rr 

- ^j^^Lot-y-Cp-y) sec^ a]da ® J ^ 

TT 

1 [T 

^ ~ TT j ^ (p + y) doi 

TT 

+ ^j^^F[ct-y-{p-y)sec^cL]da 


.(177). 


* Journal de VJ^cole Poly technique, Cah. 19, t. 12 (1823). See also V. ' 
Yolterra, Acta Math. t. 18 ; Hantzschel, deduction der Potentialgleichung , Ch. 1. 
t These formulae are not given in Prof. Lamb’s paper. 
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^3 


Each of the integrals represents a wave-function provided 
J^(± oo ) = 0 and 

TT 

tB>naF'[ct±y — (p±y)seo^a]-^0 as a— 

In these circumstances we have the solution of the diffraction 
problem for the case when the initial disturbance is represented 
by u=F(ct + y\ the magnetic force being parallel to the 
axis of z. By suitably choosing F we can deal with the case 
of a solitary wave. 

A new method of solving the problem of diffraction by 
a straight edge has been given recently by Oseen^. 

Problems connected with a wedge have been treated success- 
fully by Sommerfeld and other writers by using a certain type 
of contour integral. The fundamental solution of (163) is now 

w = ^. f (178), 

zmj da V 

where v = and the path of integration is a simple con- 

tour which starts from oo i -h 7 and goes to 00 i -h 7 ' without 
crossing the real axis. The quantities 7 , y are subject to the 
inequalities 

27r > 7 > TT, 0 > 7 ' > — TT. 

This function u is multiform and of period 27i7r, but on an 
n-sheeted Riemann’s surface with the origin as branch-point 
and the line — (tt — ^ 0 ) as branch-section, it is uniform. 

With the aid of this function a number of diffraction 
problems may be solved. 

Thus in the case of a perfectly conducting prism of angle 
27r “ a if the electric force in the incident waves is parallel to 
the edge of the screen and is represented by the real part 
of the expression 

Qtk [ci ^-p cos ■— ^ 0)3 

it can be shown by an extension of the method of images, that 
for the total disturbance the electric force is the real part of 
the expression -f 

= 

* Arkiv for inatematik, a&tronomi ochfysik (1912). 
t Macdonald, Electric Waves, p. 192 (1902). 
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Wi = COS — — cos — (<56 + <f>o), 
a a ^ ^ 

and the path of integi’ation is the same as before. 

In the associated problem when the magnetic force in 
the incident wave is parallel to the axis of z, the magnetic 
force for the total disturbance is the real part of the expression 

” 2^‘ f ^ ^ (180). 

These solutions and the solutions of analogous problems have 
been discussed by W. H. Jackson H. M. Macdonald f , F. 
Reiche A. Wiegrefe §, and other writers. 

§ 31. Elliptic coordinates ||. 

If we put 

^ == cosh o) cos %, y = sinh o) sin % (181), 

the differential equation (163) becomes 

^2/1 y ^2/1 y 

0^ + ^ (cosh^ oj — cos^ %) '^ = 0 ( 182 ). 

The elementary solutions are now of the form 

u=:E (a))F{x), 

where E and F satisfy the equations of the elliptic cylinder IT 

(183). 

solutions of these differential equations have been 
obtained recently by Prof. Whittaker 

* Proc. London Math. Soc. Ser. 2, Vol. 1 (1904), p. 393. 
t Ihid. Vol. 12 (1913), p. 430. 
t Ann. d. Phys. Bd. 37 (1912), p. 131. 

§ Ibid. Vol. 39 (1912), p. 449. 

11 H. Weber, Math. Ann. Vol. 1 (1869) ; Mathieu, Louville^s Journal, Ser. 2, 
Vol. 13 (1868) ; Hartenstein, Hopped Archiv (2), t. 14, p. 170 ; R. C. Maclaurin, 
Cambr. Phil. Trans. Vol. 17 (1898), p. 41. 

H For this equation see Heine, Handbuch der KugelfunUionen ; Lindemann, 
Math. Ann. Bd. 22 ; Hantzschel, Zeitschr. f. Math. u. Phys. Vol. 31, p. 25 ' 
(1883) ; Mathieu, Liouville^s Journal, Ser. 2, Vol. 13 (1868). 

** Math. Congress, Cambridge (1912). 


cosh^o) +p)E= o| 

^2-(^CO3»%+Jp)^=0j 


where 

2^ = cos — — cos - (6 — 6X 
a a ^ ^ 
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Elliptic coordinates are appropriate for the solution of 
problems connected with the scattering of electromagnetic 
waves by an elliptic cylinder*. 

W. Wien has suggested that the problem of the diffraction 
of light through a straight slit in a screen J may be treated 
with the aid of elliptic coordinates by regarding the screen as a 
limiting case of a hyperbolic cylinder. 

H. Weber § has shown that wdien the elliptic and 

parabolic substitutions are the only transformations which lead 
to elementary solutions of the equation (163). For further 
properties of this differential equation we may refer to Pockels, 
Uhey^ die partielle Differentialgleichung Au 4 - k-u = 0, Teubner, 
Leipzig (1891), and to Lord Rayleighs Theory of Sound, 

§ 32. Other diffraction problems. 

The diffraction of light and electric waves by a grating of 
wires is a problem of importance, but the mathematical treat- 
ment is very difficult and the theories that have been given so 
far are of an approximate character. Sir J. J. Thomson || has 
discussed the theory of Hertz's grating IT which consists of 
a number of parallel equidistant metal wires. When electric 
waves whose wave-length is large compared with the distance 
between the wires fall normally on the grating, they pass 
through if the electric force is at right angles to the wires but 
are reflected if the electric force is parallel to the wires. Prof. 
Lamb ** has considered the case of a grating which consists of 
parallel strips of metal ; his theory has been supported by the 

* See, for instance, K. Aichi, Proc. Tokyo Math. Phy». Soc. (2), 4, j). 266 
(1908) ; B. Sieger, Ann. d. Fhysik (4), Bd. 27 (1908), p. 626. 

t Jahresbericht d. deutsch. Math. Verein, Bd. 15 (1906), p. 42. 

t For this problem see K. Schwarzschild, loc. cit. ; Lord Kayleigh, PhiL 
Mag. YoL 43 (1897), p. 259 ; Scientific Papers, Vol. 4, p. *283 ; Proc. Hoy. Soc. 
A, Yol. 89 (1913), p. 194. An interesting experimental result has been obtained 
recently by P. Zeeman, Amsterdam Proc., Nov. 28 (1912), p. 599. 

§ Math. Ann. Bd. 1. See n\mPLnntz^Qh.o\, Reduction der Potentialgleichmuj, 
p. 137. 

I! Recent Researches (1893), p. 425. 

If Collected Works, Yol. 2, p. 190. For some recent experimental work see 
H. du Bois and H. Bubens, A7m, d. Phys. Bd. 35 (1911), p. 243 ; A. D. Colo, 
Phys. Review, Jan. (1913). 

** Proc. Lo7idon Math. Soc. Ser. 1, Yol. 29 (1898), p. 523. 
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experimental work of Cl. Schaefer*, J. Langwitz*, and G. H. 
Thomson f. 

Lord EayleighJ has given an approximate electromagnetic 
theory of the action of a grating on waves of light and this 
theory has been extended by W. Voigt § so as to take into 
account the properties of the material of which the grating 
is made. Voigt's theory has been tested experimentally by 
B. Poganyll who gives an account of previous experimental 
work on the subject. 

The diffraction of light through a circular hole in a screen 
is a problem of interest to mathematicians which has yet to be 
solved IT. A promising method of attack is to regard the screen 
as the limiting case of a hyperboloid of revolution of one sheet. 


Examples. 1. Prove that 




\TT / J 


ihcr J 
e dr 


'00 ^c{t - t) -h pcoa{(j>- fo) 


2. If x+i^=a cosh (co-f-^x)? r)=^kae'‘*^i 

y (kp) ii) (,) 

” 3 2 ‘ 

where the summation extends over all even integral values of n if v is even 
and over all odd integral values of if v is odd. 

(J. H. Hartenstein, Qrunerfs Archiv (2), 1. 14, p. 170.) 


§32a. The introduction and elimination of discontinuities. 

Wave-functions with singular lines or with singularities 
travelling along straight lines with the velocity of light may 
sometimes be employed with advantage in the solution of 
diffraction problems. To illustrate the method to be adopted 
we shall consider the diffraction of waves of sound by an 

* Ann, d. Phys. Bd. 21 (1906), p. 587. The theory is developed from a new 
point of view by C. Schaefer and F. Eeiche, Ibid. Bd. 32 (1910), p. 577 ; Bd. 35 
(1911), p. 817. 

t Ibid. Vol. 22 (1907), p. 365. 

t Proc. Roy. Soc. A, Yol. 79 (1907), p. 532. 

§ Gott. Nackr. (1911). I| Ann. d. Phys. Bd. 37 (1912), p. 257. 

H Approximate solutions have been given by G. G. Stokes, Camh. Phil. ^ 
Trans. (1849); H. Lorenz, Vidensk. Selsk. Skr., Copenhagen (1890); H. A. 
Rowland, Amer. Journ. Vol. 6; A. Grimpen, Diss. Kiel (1890) ; A. E. H. Love, 
Phil. Trans. A, Vol. 197 (1901). 
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infinitely thin semi-infinite plane bounded by a straight edge^. 
We shall suppose that the waves are sent out from a stationary 
source and that the screen acts as a perfect reflectorf. 

Let the axis of z be taken along the edge of the screen and 
the axis of x in the plane of the screen at right angles to the 
edge. Let P be the source of sound and Q an arbitrary point 
on the edge of the screen. If the disturbance issuing fi^om P 
were reflected according to the laws of geometrical optics, the 
total disturbance would be discontinuous in crossing two semi- 
infinite planes, each of which is bounded by the edge of the 
screen. The first of these planes is a boundary of the geometrical 
shadow, when continued across the edge of the screen it passes 
through P. The second plane is the boundary of the geometrical 
shadow for the optical image of P, viz. Pi. 

To obtain the correct solution of the diffraction problem we 
must add to the disturbance just described a second one having 
discontinuities which will annul the above-mentioned dis- 
continuities, the new disturbance must also be chosen so that 
the boundary condition is satisfied at the two faces of the 
screen. We shall now show that the required disturbance can 
be built up by superposition from elementary disturbances 
with singularities along lines such as PQ and P^Q produced. 

Let P be the distance of an arbitrary point {x, y, z, t) 
from Q, then if (0, 0, f) are the coordinates of Q and c is the 
velocity of sound, we know that a function of type 
1 ./ R x-\-iy \ 

rJ c ’ z^^^R) 

satisfies the wave-equation. Let us choose the arbitrary function 
f in such a way that the expression becomes infinite along the 
line PQ produced and returns to its initial value when the 
point X, y, z is rotated twice round the edge of the screen. 
This last condition is added so as to enable us to satisfy the 
boundary condition. 

* This problem has been solved by H. S. Carslaw, Proc. London Math. Soc. 
Vol. 30 (1898), p. 121. A transformation of bis solution suggested the method 
described here. 

t This assumption is usually justifiable. In Prof. A. G. Webster’s experi- 
ments on the reflection of sound from the ground {Phys. Tteview, Vol. 28 (1909) » 
p. 65) it was found that the reflection is more than 90 7o* 
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Now if we write 


X + iy 


— e 


-u+i<f) 


2! Jt 

where it and ^ are real quantities, and use Rq, Uq, (j^o to denote 
the values of R, u, cj) for the point P, it is easy to see that the 
function 


^ (if - ^) sec ^ + i {u + ?(„)] 

satisfies the requirements, for it is periodic in (f> with period 
47r and is infinite along PQ produced, where 


u==—Uoj ^==^o+7r. 

We now imagine sources corresponding to wave-functions of 
this type to be associated with each element of the edge, 
and suppose the strength and phase of the source at Q to 
depend on its position relative to P in such a way that 


(t ^ /-I 

( P -}- Rq\ 


r c j 


where /(i) is the strength of the source at P at time t 
In this way we obtain an integral 

° i.L 

which will be shown to be discontinuous in the way required 
as the point (x, y, z) crosses the boundary of the shadow for P. 
In a similar way we can construct an integral 

i ^ i [-^ + -^0 - i (« + «o)], 

which can be shown to be discontinuous in the way required as 
the point {x, y, z) crosses the boundary of the geometrical 
shadow for Pi. 

Now let r, be the distances of the point x^ y, z respectively 
from P and Pi, then the velocity potential V of the total 
disturbance is given by the following expressions in the different 
regions of space 

+ inS., 

F=-Fi-F, 
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The space Si is bounded by the screen and the limiting 
plane of the geometrical shadow for Pi, So is bounded by the 
limiting planes of the shadows for P and Pi, Ss is bounded 
by the screen and the limiting plane of the shadow for P. 

dV 

The boundary condition is that ^ should be zero over the 

two faces of the screen and it is easy to verify that this 
condition is satisfied. To show that V is continuous for the 
whole of the space outside the screen and vanishes at infinity 
when the function f is finite, we shall transform the integrals 
Vi, Vo to the forms given by Prof. Carslaw. To do this we put 
u + Uo = b, then if p is the distance of a point from the axis 
of we have 

^ f + P, Roy 

R^ = p^ + (s-^ R^ = + {zo ~ 

= P poe-^<^ = Ro-Zo+ 

ppo cosh {u + ^o) = 0 (^0 - ?) + RRo, 

+ po^ + (^ - ZoY + 2ppo cosh (w 4- -Uo) == + RoYy 

d-w + duo = - cZf . 

Hence it follows that 


r^•- 


-1 

4<7r j -.a 


db 


On substituting the expression for R -I- Ro in terms of b we 
obtain an integral which is equivalent to the one given by 
Prof. Carslaw. To see that it is discontinuous* we use and FI 
to denote the values of the integral for <^ = tt + <^o + e and 
<^ = tt + ^0 ~ ^ respectively, where e is a small quantity. The 
difference between these quantities may be regarded as a 
contour integral and can be evaluated by Cauchy's theorem 
We may write 

Taking tke residue for f = 0, i.e. ^ = tt + 6 = 0, we get 


V^-V_ = - 


R + R, 




* A more careful proof is given in Prof. Carslaw’s paper. 
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for the conditions (f) = 7r -{- <f>Q, 6 = 0 imply that the radii 
R and Ro are in one straight line and so give r when added 
together. 

It is now clear that the integral Fj possesses the right type 
of discontinuity and a similar remark holds for the integral Fg. 
The method can no doubt be modified so as to give solutions of 
other types of diffraction problems, the chief difficulty arises 
in the choice of a function which will satisfy the boundary 
conditions. At any rate the method suggests an interesting 
type of boundary problem in which the desired wave-functions 
have specified discontinuities instead of being continuous 
everywhere. This type of problem ought to be studied more 
completely. 

In the general problem of the .diffraction round a moving 
object of the waves issuing from a moving source, the wave- 
functions that are derived by the methods of geometrical optics 
have discontinuities at a certain boundary which is the locus of 
points travelling along straight lines with the velocity of light. 

The points in question start from certain points of the 
moving object and move along tangents to the surface of the 
object, their paths being in fact continuations of the paths of 
particles that may be considered to have been emitted from the 
source. Indeed, if we imagine the source to emit particles in 
all directions as it moves about, the particles which just graze 
the moving object will, when they continue their rectilinear 
motion with the velocity of light, form the boundary at which 
the discontinuities arise. 

In Chapter viii we shall obtain a class of wave-functions 
with singularities moving along straight lines with the velocity 
of light. These functions seem to be just the ones that are 
required for the building up of wave-functions with dis- 
continuities of the type just described. The problem of forming 
in this way the functions which will enable us to complete the 
solution of the diffraction problem is one which awaits solution. 


CHAPTER VI 


TRANSFORMATIONS OF COORDINATES APPROPRIATE FOR 
THE TREATMENT OF PROBLEMS CONNECTED WITH A 
SURFACE OF REVOLUTION 


§ 33. Spheroidal coordinates. 

Problems in which there is symmetry round the axis of z 
can often be treated with the aid of a substitution of the form* 

p + iz =^f{oL + i^) (184). 

Taking a, /S, as orthogonal coordinates, we have 

+ dy^ + dz^ — (da^ + d^^) . .(185), 

and equations (18) of § 8 become 

7 [a^O-) - • 

S 53 5? (t’’***) “ 


where 


p. 9 («. 

d(p,^)- 


* This substitution has been used in other branches of mathematical 
physics by C. Nehmann, Theorie der ElektricitdU- und Wdrme-Vertheilung in 
eineni Binge (1864); E. Mathieu, Cours de ‘physique matheniatique (1873); 
A. Wangerin, Berliner Monatsherichte (1878) ; Hantzschel, Reductimi der 
Potentialgleiclmng ; Michell, Mess, of Math. (1890) ; Basset, By dr ody nannies, 
*Yol. 2, p. 8; F. H. Safford, Amer. Journ. Vol. 21; Archiv der Math. Bd. 13 
(1908), p. 22. The important developments on which the following analysis is 
founded are contained in papers to which we shall refer presently. 
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These equations may be satisfied by putting 


,, JdD, 

p dp’ 


Jda 

'p da’ 


M^=±-n ...(187), 
p 


/dp dfl dp dfl' 

r nTR /!!i 


+j.n. 


: 0 ...( 188 ). 


where Sl—U±iV is a solution of the partial differential 
equation 

0“O 1 fdp 0O dp 0fl\ , ^ o 0 ss^ 

The problem of finding the periods of free electrical oscilla- 
tions on a conducting spheroid is of considerable interest 
because a straight rod of circular cross-section can be regarded 
as approximately equivalent to a prolate spheroid whose major 
axis is relatively much longer than the minor axis. This 
problem has been treated very fully by M. Abraham* * * § , 
R. C. Maclaurinj", M. Brillouinf, F. Ehrenhaft§ and J. W. 
Nicholson ||. The effect of a spheroidal obstacle on a train of 
waves has been studied by K. F. HerzfeldIF. 

For prolate spheroids the appropriate substitution is** 

z + ip = a cosh (a + i/S) (189), 

or ■ p== a sinh a sin /3, z = a cosh a cos yS, 

o / (cosh® a - cos® /9) (190). 

a {a, fiS) 

The partial differential equation is now 


0®fl 


■ coth a - 


a® A;® (cosh® a — cos® /3) 

(191), 


and there are elementary solutions of the form fl ^ 
where A and B satisfy the differential equations 


:A{a)B{0) 


coth a ■ 


(/b®a® cosh® a + X)A 


^ - cot /? ~ + (\ + A:® a“ cos^ 5 = 0 


...(192). 


* JOissertation, Berlin (1897); Ann. d, Phys. Bd. 66 (1898), p. 435 ; Math. 
Ann. Bd. 52 (1899), p. 81. 

t Camhr. Phil. Tram. Yol. 17 (1898-9), pp. 41 — 108. 

t Propagation de PSlectricitS (1904), Ch. vi. 

§ Wiener Berichte (190^, p. 27^. 

\\ Phil. Mag. {190Q). ^ Wiener Berichte (lQll) , p. 15S1. 

** Of. Heine, CreWe, Bd. 26 (1843), p. 1S5] Kugelfunktionen, Bd. 2, § 38 ; 
liomWB Hydrodynamics tp, 
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These equations are discussed in some detail in the papers 
to which we have just referred, they may be reduced to a special 
form of an equation obtained by Prof 0. Niven* in a study of 
the conduction of heat in ellipsoids of revolution. 

For oblate spheroids the appropriate substitution is 

a cosh {a + i/3) (193), 

giving p = a cosh a cos / 3 , 2 := a sinh a sin yQ, 

(cOSh^ OL — cos^ ^). 

9 (a, 

The surfaces /3 = const, are now hjrperboloids of one sheet, 
the surface yS = 0 can be regarded as the surface of a screen 
which is pierced by a circular hole of radius a. 

The partial differential equation for X 2 is now 

a^ + ““ ^ ^ (smh^ a + sin^ /?) O = 0, 

and there are elementary solutions of the form 12 = (a) J 5 (/ 3 ) 

where 

rl-A A 

-j-j — tanh a . — h (X + sinh“ a) ^ = 0 

^ L(IM). 

^ +ta,n/3 .j^+ia^k^shoP ^-X) B = 0^ 

When 12 is independent of t, i = 0 and the elementary 
solutions are of the form 

^=IPnQ)d^iPniv)dvy |^ = cosha, i; = cos/3 ...(195) 
for prolate spheroids, and of the form 

^ == fPn d^JPn (v) dvy sinh a, 97 = sin yS. . .(196) 

for oblate spheroids. In either of these solutions a function 
can be replaced by Qn- The corresponding solutions of Laplace's 
equation are of the type 

F= [APn (f) + BQn m [GPn iv) + DQn iv)l • .(197), 
•where J., P, 0, D are arbitrary constants. 


B. 


Phil. Trans. (1880), p. 138. 


7 
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Examples. 1. Prove that 

[(cosh a cos /3 — cos yf + sinh^ a sin^ jS] ~ ^ 

i (2«.+l) (cosh a) P,i (cos Pa (cos y). 

0 

(C. Neumann.) 

2. Prove that a function B O) which satisfies the differential equation 
(192)Jand is zero for /3=0, ^—ir is a solution of the homogeneous integral 
equation 

E(e)=pam^dJ^ g-il»oosecoss R{ik)^Q, 

where’fi is determined by the condition that the integral equation should 
possess a continuous solution which is not identically zero. 

(M. Abraham.) 

3. If A (a) be defined by the equation 

A (a)=sinh2a 

its satisfies the differential equation (192). A second solution of this 
equation is given by 

I(a)=smh2a R(ik)>0, 
and is suitable for the representation of divergent waves. 

(M. Abraham.) 




§ 34. Paraboloidal coordinates. 

If we write 

z+ip = {a^ + i^^'f, /0o“ = /3 

so that the transformation is 

z=-a-y3, p = 2V^/y (198), 

the differential equation (143) becomes* 


«3^ + (m + l)^-,3 


/ , 18“IV , 


and is satisfied byf 

W==A(ol)B{J3) 


* Of. H. J. Sharpe, Quarterly Journal, Vol. 15 (1878) ; Proc. Camb. Phil. 
'Soe. Vol. 10 (1899), p. 101; Vol. 13 (1905), p. 133; Vol. 15 (1909), p. 190; 
H. Lamb, Proc. London Math. Soc. Ser. 2, Vol. 4 (1907), p. 190. 

+ The existence of elementary solutions for the paraboloid and certain • 
other surfaces is established in Bdcher^s Die Reihenentwickelungen der Potential- 
theorie. See the table on pp. 256-7. 
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if 




(199), 


where h is arbitrary. 

Putting 2ikn = ik (m -f 1) ~ /i we find that the differential 
equations are satisfied by putting 

A = {2ika), B = {^ikj3), 

where Fm^^{s) satisfies the differential equation* 

s + (m + 1 — 5 ) ^ + nF=0 (200). 

When 7^ is a positive integer, one solution of this equation 
is furnished by Sonin's polynomial f Tm^{s), which may be 
defined with the aid of the expansion 


si 

(1 + = i r (m + w + 1) (s). . .(201). 

n—0 


A few properties of this function are given here for the sake 
of reference. 




(m + n + 1) r (m H- n) I n — 1 1 1 

^11-2 


+ 


r (m -h n — 1) Iti — 2 ! 2 


f (s) {s)ds=^0 n 

JO 




r(n-M) r(m+7n-i) 


p=sn 


d^ 


g,T„,\s)==T::z(s) 


...( 202 ), 


..(203), 


.(204), 


^ (s)] = (s) (20.5), 


This is a slight modification of Weiler’s canonical form for an equation of 
Laplace’s type, Grelle^s Journal, Bd. 51 (1856), p. 105. The equation is discussed 
for real values of m and n by 0. Schl6milch,Brf;7tcnm.dnrtZy«i.s-, Bd. 2 (1874), p. 517. 

i- Math. Ann. Bd. 16. Further properties of the function are given by 
L. Oegenbauer, Wien. Ber. (1887), p. 274, who proves that the roots of the equation 
T^'"'(s)=0, considered as an equation of the nth degree in a, are all real, 
•positive and unequal. This is a generalisation of the result obtained by 
Laguerre for the case m=0. A geometrical proof has been given by Bocher, 
Proc. of the Amer. Acad, of Arts and Sciences, Vol. 40 (1904). 


7—2 
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(- l)”e* 




r (m + n + 1 ) r (n + !)&’ 


[e“«s™'+’‘]...(206), 


(«) =T>r+i) 

Equations (201), ( 202 ) and (203) were given by Abel* * * § and 
Murphyf for the case m = 0: the polynomial is then equivalent 
to the polynomial of Tchebycheff| and Laguerre§ which occurs 
in the theory of interpolation and also in the theory of continued 
jfractions. When m = ± the polynomial can be expressed in 
terms of the polynomial ?7^ (^) discussed by Tchebycheff || and 
HermitelT, or in terms of the function of the parabolic cylinder, 
discussed by Weber**, Whittaker ff and others 

The above analysis indicates the existence of a wave-function 
of the form 

n = (2ika) ( 2 */ 3 ) p'^ ... (208). 


This function can be expressed as an integral of the form 
used in § 5, we have in fact the equation 

(kprTrr,^(2iJm)Trr:^{2ik/3) 

(- ir 


/.Btt 


2ttT (m + + 1 ) 


[ To^ [— 2k (p cos 7 -f ^>)] e~‘'^'^ydy 
J 0 

(209), 


from which the required representation can be immediately 
derived. In this formula m is either zero or a positive integer. 
The convergence of a series of terms of type (208) in which 


* Mimoires de matMmatique par N. H. Abel, Paris (1826) ; Oeuvres ^ Sylow 
and Lie, t. 2. 

t Gambr. Phil. Trans. (1833). 

X Mem. de VAcad. de St PUersbourg (1860). 

§ Bull, de la Soc. math, de France, t. 7 (1879) ; Oeuvres de Laguerre, t. 1, 
p. 428. 

11 Loc. cit. See also Sturm, LiouvilWs Journal, Vol. 1. 

IT Gom^ptes Bendus, t. 58 (1864), p. 98. The Hermite functions have been 
generalised by Curzon, Proc. London Math. Soc. Vol. 13 (1914), p. 417. The 
generalised functions are intimately connected with the functions considered here. 

Math. Ann. Bd. 1 (1869), p. 1. 

tt Proc. London Math. Soc. Ser. 1, Vol. 35 (1903), p. 417. 

XX Baer, JDiss. Gilstrin (1883) ; Hantzschel, Zdtschr. fUr Math. Bd. 33 (1888) f 
Adamoff, Annales de St Pitersbourg, t. 5 (1906) ; G. N. Watson, Proc. London 
Math. Soc. Ser. 2, Vol. 8 (1910), p. 393. 
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n takes different integral values can be partially discussed ■with 
the aid of the equation 

{ix) Tm”‘ (- isc) 

1 T Fl I + n + j 

r(K + l)r(m+l)J L ■^l(TO + l)^(m+2)“' 

n {n — 1) (m + + 1) (m + ?i + 2) ^ 

1 . 2 (m + 1)^ (m + 2)2 (m + 3) (m + 4) ^ 

which shows that the modulus of T^(ix) increases with x. 
Hence if a series of terms of type (208) converges absolutely 
for any given value of a, it converges absolutely for all smaller 
values of a. 

For a fuller discussion of the convergence it would be useful 
to have an asymptotic expression for T^{s) when n is large. 
Suitable asymptotic expressions have already been found for 
the case m = ± ^ by Adamoff and Watson. 

The differential equation (200) has been studied for general 
values of m and n by Pochhammer*, Jacobstahlf, Whittakert 
and Barnes§. It usually possesses two distinct solutions which 
can be expanded in power series converging for all finite values 
of 5. If, however, m and n are positive integers, there is only 
one solution which can be represented by a convergent power 
series in s, the other may be defined by the equation 

= r(^l)^‘/o (s- 

it contains a logarithmic term. For negative integral values of 
n we may adopt the definition 

is) = ( 212 ). 

It should be noticed that when | a | is large, TJn!^ {^ika) has 
an asymptotic expansion of which the first term is 

The solutions of type 

iz+ct) ±i7n4 Jj^n (^ika) (2ik/3) 

* Math, Ann. VoL 36, p. 84 ; VoL 46, p. 584* t Ihid, VoL 56, p. 129. 

X Bull. Amer. Math. Soc. (1904). 

§ Cambr. Phil. Trans. Yol. 20 (1906), p. 253. 


•••] ( 210 ), 
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are consequently suitable when a is large for the representation 
of waves diverging to infinity in the positive direction of the 
axis of z. 

It may be worth while to mention here that the functions 
2m” («), (s) both satisfy Gegenbauer’s difference equations 

K~' (s) = (m + n) F:iI (s) + FZl (s) 
n (m + n) (s) - {5 - (m + 2n - 1)} F^"^ (s) + F^~^ (s) = 0 
(n - 1) C^) = {s-^(m + n^l)} F:iI (s) ^ Ft^, {s) 

(n ^ 1) FV ( 5 ) = + 1)} KXi (s) ^ s F:11\ {s) 

(213). 

The function TJrr!^ (s) also satisfies an equation analogous to 
(204). 

I 35. Eelations between different solutions. 

Many useful formulae may be obtained by expanding known 
wave-functions in series of elementary wave-functions of type 
(208) and by identifying our elementary wave-functions with 
certain definite integrals which are known to represent wave- 
functions. For instance, we have the expansion 1 tan | o) | < 1 

^ijcz COB o) Qqp gijj 0)) 

/ . to 

\kp e^'*^sec®Y 2 + 

X tan — Tm^ {2ika) {2ik/3) (214) 


which enables us to represent a plane wave with the aid of 
a double series of solutions of the form (208). 

Further, identities may be obtained by deriving wave- 
functions from Cunningham’s solutions* of the equations 


dr~''dx^^ 

The first equation possesses the polynomial solutions 




71-1 


XT 


<(-s) 


* Proc. Boy, Soc. Ser. A, Vol. 81 (1908), p. 310, See also Wera Lebedelf, 
Dim. Gottingen (1906) ; Math. Ann. (1907). The first result is given by Appell, 
LiouviUe^t Journal, Ser. 4, t. 8(1892), p. 187. 
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and also the solutions 


72. -1-1 

T 2~ e~^T'^ 


(S- 


XT~ 2 e 


The second equation possesses the polynomial solutions of 
type 

T" sin m (</) - <^o) ( 218 ) 

and also the solutions of type 


I (^) m((f>- <j>o) ...( 219 ). 


Wave-functions may be derived from these solutions by the 
method of § 13 . 

We add here a few relations which are obtained by ex- 
pressing the solutions thus formed in terms of old solutions. 


4. ^ _ 2 cos ®] 


00 Vl 

= 2 (- 1)" r (m + n + 1) f 2 . _(220), 

^-^——6^= 2 s”*\”‘+»‘Tm’‘(s) (221), 

(— 1 9 «i+i r°° 

- r(L.+i)»i j,‘~* 

2m+iX»+m^^(Xp) 

= r (?Z 4- 1) r (?Z + m -+• 1) p"*- 

( 223 ). 

The proofs of these are left to the reader. 

Prof. G. D. Birkhoff has remarked to me that the differential 
equation (200) can be regarded as a limiting case of the h}rper- 
geometric equation when two of the singularities coincide at 
infinity*, consequently many properties of the solutions can 
be derived from known properties of hypergeometric functions f. 
It should be noticed that when W is independent of t there 

* Cf. Bdclier, Die ReihenentwicTzelungen der Potentialtheorie, p. 137. 
t This method was used in a particular case by Kummer, Crane's Journal, 
Bd. 15 (1836), p. 138. 
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are elementary solutions of equation (143) of the form W = AB 
where 

(2i V Aa)^ ’ (2^ 

We thus obtain elementary solutions of Laplace’s equation 
of the form 

Jm Va) Jm Vy3) cos m ((;& - ^o) (224), 

where m, c^o are arbitrary parameters. 


§ 36. Toroidal coordinates. 

If we put 

^ = pcos^, y = psmc[), ^=fcoshft), fsinh ft)...( 226 ), 
p + acoth^— ~ , 


«»mb^ asint ^ 

cosn <r — cos 'yjr ^ cosh cr — cos i/r ^ 

the wave-equation becomes 

_9 I sinh <7 sin du} 9 f sinh<rsini/r 9't^) 
d<r [(cosh cr — cos i/r)2 9<r] [(cosh <T — cos 9'i/r) 

sin ^fr ^ ^ sinh cr d^u 

sinh cr (cosh a — cos d<p^ sin ^fr (cosh cr — cos i/r)^ ” 

(227). 

This is satisfied by 

U = F((r) G(yfr) (cosh <r — COS 'yjr) e**" COS m (<^ — <j^)o)...(228) 


if 


cosech „ ^ (si* , g) -{„(„ + 1 ) + 

cosecV- A (3i„ 3. 1) 


Hence "we obtain wave-functions of the form 


^’==0 

F=0 

i 

(229). 


(cosh (r — cos^lr) (cosh a-) (cos i/o) e**" cos m (<^ — (^o) 

(230). 

Other solutions of the wave-equation may be obtained by 
replacing the functions PJ‘ by Q/. 

Many useful formulae may be obtained by expanding 
particular wave-functions in series of wave-functions of type 
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(230). The expansion of unity, for instance, gives rise to 
Neumann’s expansion 

co ih <r -cos^ = ^ 

It should be noticed that when we make the substitution 
(225) the wave-equation becomes 


d^n Idu 1 d^u 
dp^ pdp^ d(p^ 


d^u 1 du 


1 ^ 
■ f ^ dco^ 


= 0...(232): 


it thus possesses elementary solutions of the forms 

Kp (X^) J,Y, (Xp) cos m ((^ — <^o) (233), 

Jp {iX^) (Xp) cos m ((^ - ^o) (234). 

The expression of solutions of type (230) in terms of the 
solutions just found leads to some interesting identities. Thus 
we have the equation 

r Kp (XO (V) Jn (xa) Xf +“-«+' dX 


Jo 


_ o«+«.-n-i r (p + m + 1) r (p + 1) r (m + 1) 

“ r(n+i) 


X(C0sh or - COS ■f)P~pl„.^ (cos ■f) Pp+m-n 

p> — lj m> — l, p4'm>— 1 (235). 

Many important formulae connected with Bessel functions 
are simply particular cases of this one*. It should be remem- 
bered that 



The corresponding integral in which Kp (X^) is replaced by 
tf^p (X^) can also be evaluated in terras of Legendre functions, 
but the formulae are more complicated. The case p = — m is 
discussed by Macdonald i*. 

It should be noticed that if we write 


cosh (« — ©) = i cot yjpj cos (j3 — </>) = coth cr "j / 90 ' 7 \ 

sinh (a ~ co) = i cosec \J/, sin (0 — ^) = ±i cosech <t J ' 
three relations of type 


9 («. /3) ^ f S (a, 

9 (*= 2/) “ c 9 (a, t) 


.(238) 


* See, for instance, the formulae given by H. M. Macdonald, Proc. London 
Math. Soc. Ser. 2, Vol. 7, p. 147, and by the author, ibid. Vol. 12, Abstracts. 
f Log. cit. p. 142. 
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are satisfied and so the functions a, /3 can be used to obtain 
an electromagnetic field by the method of § 5. It is easy 
to verify that the function 

u = (cosh <7 — cos /3) (239) 

satisfies the wave-equation, /being an arbitrary function. 

We add here a few formulae for PrU^ (cosh <j), Qn^ (cosh cr) ; 
these and other formulae will be found in the memoirs of 
Dr Hobson and Dr Barnes to which we have already referred. 

(cosh cr) = \ coth^ ^ 

r (1 — m) 2 


X JP |— n, 7^ -I- 1 ; 1 — m ; — sinh^ 
92m 


, ^2. ^ — 2<ry—m g— (?i+l)<r 


r (1 — m) 

X F{^ — m,l+n — m; 1 — 2m; 1 — 

cr>0, 


(cosh cr) = (- 1)^ Vtt (1 - 

r(n-l-f) 

XjP{m-f-|, n + m + l; + cr > 0. 

Various asymptotic expansions for these functions are given 
by the authors just named and by Dr Nicholson*. 

It should be mentioned that the solutions of the wave- 
equation that have just been obtained are not directly useful 
for the treatment of the boundary problems of mathematical 
physics. They may, however, be used to construct useful 
solutions of the equation Au + 1(^—0 by means of various 
artifices. If, for instance, we multiply one of our wave-func- 
tions by and integrate with regard to t between z and 00 , 
the resulting function will often be a solution of Av -f- k^v — 0. 
This may be illustrated by taking the wave-function 
u = Jq (i\^) (\p) cos m ((f)- <po) 

and using the formula 






x= 




0 

X=> 

....(240). 


British Association Reports, Winnipeg (1909), p. 391. 
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The integration can be taken between other limits in 
certain cases ; for instance, the integral 
r'/zHlp+a)'" dt 

^ J VS^+(7-W 7{(P + -z^-{p- a)=l 

......(241) 

represents the solution of Ai; H- = 0 corresponding to a 
circular ring of sources. In this case our wave-function u 
is a constant multiple of cosh a* — cos 

The theory of electrical oscillations on a conducting anchor 
ring has been treated by H. C. Pocklington*, W. McF. Orrf 
and Lord EayleighJ, without the use of toroidal coordinates ; 
the results are of course only approximate. 


§ 37. Solutions of Laplace's equation. 
If we put 




a sinh < 


cosh a- — cos ^ 


__ a sin 
cosh <r — cos 


...(242), 


the angle yjr may be interpreted as the angle which two fixed 
points Aj B whose coordinates are 0 = 0, p=^ ± a, subtend at 

PA 

a point P (p, 0 ); the quantity a may be interpreted as log . 

The surfaces yjr = const, are spheres having a real circle (p = a, 
0 = 0) in common, the surfaces a* = const, are anchor rings. 

If we use the toroidal coordinates cr, yjr, <^, Laplace's equation 
becomes§ 

d I sinh a du 
dyjr (cosh cr — cos yfr dyfr 

1 = 0 

^ sinh cr (cosh cr — cos yfr) d(jy 


Au = 


__ 9 f sinh cr 


du 


dcr (cosh a — cos yjr dcr 


* JProc. Caml). Phil. Soc. Vol. 9 (1897), p. 324. 

+ Phil Mag. Vol. 6 (1903), p. 667. 

X Proc. Boy. Soc. Ser. A, Vol. 87 (1912), p. 93. See also C. W. Oseen, 
Phys. Zeitschr.i Dec. 1st (1913) ; Arkivfdr Mat. Ast. och Fysikf Bd. 9 (1913). 

§ B. Biemann, PartielleDiferentialgleichungen^BisittendoTcfs edition (1861); 
C. Neumann, Theorie der Elektricitdts- und Wdrme-Vertheilung in eimm Binge, 
Halle (1864) ; W. M. Hicks, Phil. Trans. (1881), p. 609; A. B. Basset, Amer. 
Journ. Vol. 15, Hydrodynamics, Vol. 2. For an alternative method see F. H. 
Saf ord, Annals of Mathematics, Vol. 12 (1898), p. 27. 
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and possesses solutions of the form 

P^_j(cosh a) 

IC == (cosh cr — cos COS n(‘yfr -- 'xItq) cos — <^o) 

Q;;^_,^(cosh^) 

(243), 

which are suitable for the treatment of problems connected 
with the anchor ring, circular disc and spherical bowl*. 

For problems connected with two spheres bipolar coordinates 
may be used ; the appropriate substitution isf 

— <^sin-\/r _ a sinh cr /9A/L\ 

^ cosh cr — cos \{r ’ cosh <r — cos 
The surfaces o- = const, are now coaxal spheres with imaginary 
common circle. The radius of the sphere o- = ctq is a | cosech ctq | 
and the distance of its centre from the origin is a | coth o-q |. 

The ratio of the distances of a point from the limiting points 
of the system of coaxal spheres is and the angle between the 
radii from these points is 

The appropriate solutions of Laplace's equation are now of 
the typej 

u = (cosh (T — cos [A cosh (nA^)<r + B sinh {n + o-] 

X cos m(</) - <^o) [fPn^ (cos (C0S'\/r)]...(245). 

It should be noticed that when we are using toroidal 
coordinates the function 


u = (cosh <r — cos (j) ± % log tanh ~ cos — ^frQ) 

(246) 

satisfies Laplace's equation and that when we use bipolar 
coordinates the corresponding solution is 

(cosh a -- cos yfr)^ f\ (j) ± i log tan |'lcoslii(a-a-,)...(247). 


^ See for instance E. W. Hobson, Camhr, Fhil. Tram. Vol. 18 (1899) ; C. 
W. Oseen, Arhivf&r Matematik, Bd. 2, No. 5; H. C. Pocklington, Phil. Trans. 
A, Vol- 186 (1895), p. 603. 

t W. Thomson (Lord Kelvin), LiouviUe's Journal (1847). 

J Gr. B. Jeffery, Proc. Boy. Soc. Ser, A, Vol. 87 (1912) ; G-. B. Dean, Phys. 
Peview (1912) ; G-, Darboux, Bull, des Sciences math. t. 31 (1907), p. 17. 
Another method of dealing with problems connected with two spheres is 
described by A. Guillet and M. Aubert, Journal de Physique ^ t. 3 (1913). 
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It should be mentioned here that other simple [solutions 
may be obtained by using the formulae 

/ 2 

X (cosh 0-) = — . / — ^ — cosh ncr. 

EXAMPLES. 

1. If with the notation of § 36 we write 

f 4- 2 ;p = a cosh (a + ^ ^) 
the wave-equation becomes 

^ + 2 coth 20 ^ 4-2 cot 20 (cosech^a+cosec^jS)^ 

7 ^ 2/1 

+ (sech^ a ~ sec^^) = 0. 

Hence show that there are wave-functions of the form 
u^A(a)B O) e**" cos <^o), 

where a, m and < 3 E)o are arbitrary constants. 

2. Prove that if p+i2=/(a4-^j3) the wave-equation becomes 

0a2 3/32 p 0a 0/3 3)37 \p^ 302 c2 3^2 j 

and obtain elementary solutions of type A (a) -5 (/S) when 

0 + ip=a cosh (a-i-^/S). Notice that the solutions of equation (188) are not 
wave-functions, they are analogous to the stream-line functions of hydro- 
dynamics. 
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§ 38. The method of Stieltjes*. 

Wave-functions which are homogeneous functions of x,y,z,t 
may be studied with the aid of the substitution 


x — s cos 6 cos (j), s Bin 6 cos % 

j/ = 5 cos d sin cj>y ict = 5 sin 0 sin ^ 

The wave-equation in these coordinates has the form 



d“U 3 1 d^u 1 d^u 

ds^ s 05 ^ cos^ 0 d<j>^ 


1 d‘^u cot 6 — tan O du ^ 
s^sin^Odx^’^ dd 

Putting cos 2$ == yet, we find that there are elementary solutions 
of degree 2n of the form 

(249), 


if 


A 

dfjb 





+0 



2!_l = o 

2 ( 1 +/.) 2 ( 1 -/.)! 

(250). 


This equation is satisfied by 

/n . nT/t S [, 1 , 'm+p m+p I - fJ\ 

(1 + /.)- (1 — /.)- ^w + lH 2 — > ~2 

(251), 

with the usual notation of the hypergeometric function t. 


* Comj^teB Rendus, t. 95 (1882), p. 901 ; Liouville^s Journal, Ser. 4, t. 5 
(1889), p. 55- See also Tisserand, Traite de micanique eileste, Paris (1889) ; 
H. Bateman, Proc. London Math. Soc. Ser. 2, Vol. 3 (1905), p. 111. 

t This gives a polynomial if» either w + l+ or — J"- — - n is zero or a 

ji It 

negative integer. 
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It should be noticed that if we write 
s 0 = sec a, sin 0 = i tan a, sin % == cosh w, cos ^ = — i sinh w, 
fL = sec- a 4- tan^ a, 

^ obtain real wave-functions of the form 


u = 5-^ ©(//.) cos — (f)o) (252), 

variable in the hypergeometric function is now — tan^ «. 
'^hen m—p the equation for @ is the equation satisfied by the 
sociated Legendre functions. We thus obtain wave-functions 
■ the form 


u = ^ (cos 2^) +x) (253). 


omparing this with the elementary solution of Laplace’s equa- 
on in polar coordinates (r, 6, <^), we see that if /(r, 0, <}>) is a 
•lution of Laplace’s equation f (s% 20, <^ + %) is a wave-function, 
may thus derive wave-functions from harmonic functions ; 

L particular, the fundamental harmonic function ^ gives us the 


indamental wave-function \ = -v-- — — r . We have 

ready remarked in § 13 that Lord Kelvin’s method of inver- 
on may be extended to wave-functions, it is easy to see that 
le result is an immediate consequence of the fact that the 
ifferential equation (250) is unaltered when — (?^^-l) is 
ritten in place of n. 

It is easy to see that there are (n 4- 1)® linearly independent 
olynomial solutions of degree n, for a general polynomial of 

egree n contains - (^^ + 1) (?2 4- 2) (?^ 4- 3) coefficients and when 


\is is operated on with H the vanishing of the resulting 

olynomial of degree n — 2 gives ~ (n — 1) (n 4- 1) conditions. 

'he difference between these two numbers is (7^4- 1)^ 

A polynomial solution of degree n is given by the integral 

/ ’2ir 

{x cos a 4- 2/ sin a 4- iz)^ (x sin a — y cos a — da, 

0 

k, set of (n + iy linearly independent polynomials is obtained 
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by allowing m and p to take the values (0, 1 . . . n). A better 
set of solutions is obtained by using the integrals of type 

u = [ [(iT — iy) + ci)]^ 

Jo 

X [(x 4“ iy) -f i (-s' + ct)Y~^ 

Polynomial solutions may also be obtained by differentiating 

the fundamental wave-function ^ and using generalised inver- 

s 

sion*. The polynomial solutions were first discussed by 
Cayleyf. WaelschJ has recently studied them from a new 
point of view- 

Example. Prove that when n is a positive integer 

Fn (cos 2^) = ~ (cos ^ cos ^ + i sin B cos dx- 

§ 39. The method of Green§. 

Homogeneous solutions may also be investigated with the 
aid of Green’s substitution 

^ = 5 sin a sin /3 cos y — s sin a sin /3 sin < 5 ^> 

2 ? = 5 sin a cos /3, ict = s cos a 

The wave-equation now becomes 

3 01 ^ 2 ,016 1 0% 

05^ s ds 5® da 0a^ 

1 0 f . , 1 /e>KK\ 

suF a sin /3 0j8 d^J sP- sin^ a sin^ /3 d<fF ~ 

and possesses elementary solutions of the form 

u^s'^ A{ol) B (y8) cos m (^ — 

* CL F. Didon, Annales de vAcole Normale (1), t. 5 (1868), p. 229; t. 6 
(1869), p. 7; t. 7 (1870), pp. 89, 247; P. AppeU, Rend. Palermo, t. 36 (1913) ; 
K. de F4riet, Comptes Rendus, Nov. 17th (1913). 

t Liouville^s Journal, t- 13 (1848) ; Phil. Tram. Yol. (165) ii. (1875), p. 675, 
See also Hermite, Oeuvres, t. 2. 

$ Deutsche Math. Verein, Bd. 19, p. 90. 

§ Camhr. Phil. Tram. Vol. 5 (1835), p. 395 ; Collected Papers, p. 187 ; 
Cayley, loc. cit. See also Heine, Handbuch der Kugelfunktionen, Bd. 1, p. 449 ; 
Crelle, Bd, 60, 61, 62 (1862 — 1863) ; E. W. Hobson, Proc. London Math. Soc. 
Ser. 1, Vol. 24, p. 67, Vol. 25 ; F. G. Mehler, Progr. Danzig (1864) ; Crelle, ' 
Bd. 66 (1866), p. 161; C. Neumann, Zeitschr. Math. Phys. Bd. 12 (1867), 
p. 116 ; V. Giulotto, Gior. d. Mat. 39 (1901), p. 162. 
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^ + 2 cota-^+ n(n+2) — ^ ^ J. = 0 . 
•*' aa L ^ sin^ a J 


sm/3 d^ \ 
d-A 


dd 

We may thus take 

-S = Cl P,; (cos ^)-(kQm (cos /9) 


bi Pill (cosa)+62 cl + J(cosa) r 


-4 = V cosec 
where Ci, 62 arbitrary constants. 

Wave-functions of degree 25 ero. 


§40. 

If XI is a wave-function of degree — 2, the formulae 


Hz-x 


9Xi 

an 


X 3X1 

+ ct 

an ^ 

dz ■ 

■^0y ’ 

Px-- 

c dt 

dx 

an 

dx 

an 

Ey = 

II 

+ ct 

an 

dy 


an 


zdfl 

“h ct 

an j 

dy ‘ 


Ez = 

c dt 

dz J 


...(257) 


give a solution of Maxwell’s equations. 

A homogeneous wave-function of degree — 2 can, of course, 
he derived from a homogeneous wave-function of degree zero. 
If F is an arbitrary function of two variables subject to suitable 
restrictions, the integral 


Xl= F 


X sin a — y cos a — ct 


da (258) 


X cos a + y sin a -I- J 
represents a wave-function of degree zero, and when this is 

multiplied by ~ ^ wave-function of degree — 2 is obtained. We 

S“ 

add here a few particular wave-functions of degree zero : 




±W 

ct. 


tan' 


log?^ (259). 


-1^ 


X ' z ' X'-^y^ 

Electromagnetic fields which are derived from this type of wave- 
function of degree -- 2 may be generalised by writing a? — f (r), 
y — ^ (t), z—^ (t), i — t instead of x^ y, t respectively and 
integrating round a closed contour in the complex t plane. 
The integrals thus obtained can generally be evaluated by 

B. B 
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means of Cauchy's theorem. Many of the results given in the 
next chapter are suggested at once by this method and may 
be thoroughly established by a method of direct verification. 

It is worthy of note that if we write r in place of ct in a 
wave-function of degree zero the resulting function is a solution 
of Laplace's equation Aw = 0. A general solution of Laplace's 
equation of degi'ee zero can be derived at once in this way from 
the first of the solutions (259). We thus obtain Donkin’s 
formula^ 



A similar result is that if 


a^F{x,y,z,w, t) 

is a homogeneous function of degree — ^ satisfying the equation 

^ ^ ^ 1 

doc^ dy^ dz^ c- ' 

and s be written in place of w, the resulting function is a 
wave-function. Now if f{x,y,z) is a solution of Laplace's 
equation, the function 

^w — ct — w — ct' w-ctj 

satisfies the requirements, consequently we may conclude that 
the function 

Ai' 

is a wave-function*]'. Other wave-functions may be derived 
from this by generalised inversion or by interchanging the 
variables y, z, ict 


* Fhil, Tram. (1857). This solution may be obtained at once from Jacobi's 
theorem that if p, q, r are three functions of u which satisfy the equation 
JJ® + 3^ + r2 z= 0 and u is defined by the equation au = xp (w) + yq (u) + zr (u), 
then an arbitrary function of w is a solution of Laplace’s equation, Werke, 
Bd. 2, p. 208. See also Forsyth, Mess, of Math. (1898). 

t This result is obtained in another way by Pockels, VUr die partielle 
JOiferentialgleichung Aw + hhi = 0. Teubner, Leipzig (1891). 
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ELECTROMAGNETIC FIELDS WITH MOVING 
SINGULARITIES 


§ 41. An electromagnetic field with a simple singularity 
or electron, first model of a corpuscle*. 

We shall now derive a family of wave-functions from the 
fundamental wave-function 1/s^ where 


= [^ - f (t)]^ +[y-v (t)]^ + [s- ^(r)f tJ. . .(262) 

and T is a variable parameter, which is at first independent of 
(T, 2 /, 0 , t. Using a method invented by Prof. A. W. Conway f 
we consider the integral 


o = ~ 


_L 

27rz J ’ 


taken round a closed contour in the plane of the complex vari- 
able T. If this contour contains only one root r of the equation 
= 0, the value of the integral is 

./(i) 

2v ’ 




where 


Z, = r (t) ~ f ( t) r (t) (^ - 0 - - 0 

(263) 


and T is the root in question. 


* I have ventured to use Johnstone Stoney’s term “ electron ” to denote the 
simple point singularity and Sir Joseph Thomson’s term “ corpuscle ” to denote 
the elementary charged particle which has been discovered by experimental 
work. 

t Proc. London Math. Soc. Ser. 2, Yol. 1 (1903). Integrals over complex 
paths had been employed previously in electromagnetic theory by Sommerfeld 
and other writers. 


8—2 
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This function v vanishes when x = ^,y = 7j, = T and 

so the wave-function XI has a singularity which moves along the 
curve r represented by 

^ = y = "n(j), (264). 

If, moreover, the velocity of this singularity E is always less 
than the velocity of light, it is easy to see that v does not vanish 
for any real values of {x, y, z, t) other than those just mentioned. 

When the velocity of the singularity E is always less than c, 
there is only one value of t less than t for which the equation 
^2 — 0 is satisfied : Xy jr, Zy t being supposed to be given. 

To prove this we surround each point ^ on F by a sphere of 
radius c (t — t) having E as centre ; then it is clear that each 
sphere lies entirely within the neighbouring one corresponding 
to a smaller value of r, provided T<t and cZp 4- drf + 

This shows that one and only one of these spheres passes 
through a given point of space and so there is only one value of 
T< t for which the equation 

f (t)P + [y^7j (t)P + ^ r (t)? r )\ . .(265) 

is satisfied^. 

Now let a point Q (xr, y, z, t) move with a velocity less than c 
along a curve 0 and let us consider the variation of r with t. 
As t increases from t tot dt the radius of the sphere associated 
with each point E will increase by cdt and since Q moves a 
distance less than cdt in the interval dty its new position will 
lie within the new sphere associated with the time t. Conse- 
quently the new position of Q lies on a sphere associated with a 
greater time r. 

Hence if Q moves in any manner with a velocity less than the 
velocity of light, r increases with t 

Things are quite diiferent when the velocity of E is greater 
than c. The spheres then have a real envelope and there may 
be more than one sphere through a given point in space, also t 
may sometimes decrease when t increases. 

In this case, however, v vanishes for values of x, y, z, t other 
than x = ^y y = 7), z—^y ^ == r, and so II has oo ^ singular lines 

* Cf- A. W. Conway, loc, dt . ; H. Bateman, Manchester Memoirs (1910) ; 
G. A. Scb-ott, Electromagnetic Radiation (1912). The theorem is due to Li4nard. 
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through each point E. These singular lines form a right 
circular cone whose axis is along Es direction of motion : each 
singular line is described by a singular point that travels with 
the velocity of light. 

When the velocity of the singularity is less than o we can 
obtain a solution of Maxwell’s equations having the moving 
singularity by using the potentials* 

A A a =??Kj) 

^ 47 ri' ’ 4i7rv ’ ^ irrv ’ 4irrv 

(266). 

It is easy to verify that they satisfy the relation 

div^+I^=0 (267). 

0 (Jv 

When the electric and magnetic forces are calculated from 
these potentials with the aid of the formulae 

H = votA, ^=-i^-grad«I> (268) 

it is found thatf 

_ e d(T, a) p _ e d (r, <r) 

“ 47r a (y, z)’ ” 47rc d (x, t) 


where 

av=:r(^-^) + v"(y-v) + r{^-0-(^'^+v'^ + n + 

(270). 

It is clear from these equations that the magnetic force is 

* A. Li4nard, U^clairagc ^leetrique, Vol. 16 (1898), pp. 5, 53, 106. See also 
E. WiecLert, Arch, n^erlandaises (2), Vol. 5 (1900), p. 54 ; K. ScLwarzscbild, 
Gott. Nadir . (1903). The potentials are usually written in the form 

(|> f — !■■■ A = — u , where the square bracket indicates 

that the quantity enclosed is to be calculated at time T=i--. Of. H. A. 

Lorentz, The Theory of Electrons, p. 50. To obtain a model of a corpuscle we 
must write de instead of e and integrate over a small region. 

t These expressions for the components of E and N were communicated to 
me by Mr E. Hargreaves in 1909 ; they should be of some historical interest 
in connection with the general theory of § 5. This was, however, the outcome 
of some independent work. Cf. Proc. London Math. Soe. (2), Vol. 10 (1911), 
p. 96. 
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perpendicular to the electric force and also to the radius 
from the effective position of E', for we have the relations 


3t ^ 


dr 

' ay' 

dr 


‘-y-v 


: — (35 — t ) 


.(271). 


It also follows jBrom these relations that r satisfies the char- 
acteristic equation 



This is to be expected because, as Jacobi has remarked^ for the 
case of Laplace’s equation, the argument t of an arbitrary func- 
tion occurring in the solution of a partial differential equation 
must satisfy the partial differential equation of the character- 
isticsf. 

To prove that there is a constant charge e associated with 
the singularity of our electromagnetic field we shall calculate 
the integral of the radial component of E over a sphere having 
the singularity as centre. We have to evaluate the integral 


ce ff/dcrdr dcrdr dcr dr 1 do- dT\ 

4i7r JJ [da; dx ^ dydy’^ dz dz dt dt) ’ 
which is easily transformed into 


Transforming the axes so that the axis of z is in the direction 
of motion of the singularity, we may put 

j/ = r (t; cos ^ - c), dS == sin ddddcj), + ^ = 

and our integral becomes 

ce P p”’ (c^ — v^) sin 0 ddd<^ __ 

^rrJoJo {vcosd — cy 


* Werke, Bd. 2, p. 208 ; Crelle^s Journal, Bd. 36 (1848). 
t For the general theory of characteristics see Hadamard, Propagation de$ 
Ondest (1903), Chapters vn. and vin. ; J. Coulon, Comptes Rendm, t. 128 
(1899), p. 1386 ; A. V. Backlund, MatK Ann. Bd. 13 (1878), p. 411 ; J. Bendon, 
Comptes Rendus, t. 124 (1897), p. 124. 
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Prof. E. T. Whittaker* has calculated potentials P, 11 from 
which A and ^ can be derived by using (7). He finds that 

n=(o,o,5f), r=(o,o,W), 

where 


= e sinh~‘ ^ , 4'rrJV = — e tan"' 

^ttK = — e log V {x — + (y *- ■)?)“ 


y-v 


.(273). 






It may be verified without difficulty that the functions are 
wave-functions. This result is a particular case of the 
following general theorem. 

If f{x, y, z) is a homogeneous f unction of degree zero satisfy- 
ing Laplaces equation Au = 0, the function 


(r), (274) 

is a wave-function. 


§ 42. The electromagnetic field due to a moving doublet. 
Let us now derive an electromagnetic field by superposing 
two electromagnetic fields of the type just described wherein 
the singularities move along the two neighbouring curves 

y=v(ry ^=?(7'), 

(ti) + €« (ti), y = V (t"i) + eyS (ti), 2 =^ (ti) + 6y (ti), 

€ being a quantity whose square may be neglected. 

If Ti is defined in tenns of x, y, z, t by the equation 

~ ^ (n) - ea (ti)] 2 + [y - 1? (ti) ~ €y8 (ti)]^ 

+ ? (ti) - 67 (ti)]2 = c" (it - Tff 

and Ti = T -f eO, we easily find that 

1/0 -f a (iJ? — f ) + /3 (y — 97 ) 4* 7 (. 2 ^ — 0 = 0. 

Also if Vi is the quantity corresponding to v, we have 

T/i = v + e [dvo- 'h a' (x — ^) + 0'{y —v) 4 - 7 ' (^r — 

— ^7)' - — v e[6v(r^ p\ say. 


Now if 


~ 47rv, ’ “ 4'rrvy ’ 


Proc. London Math. Soc. Ser. 2, V'ol. 1. 
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we find that 
(Xa; = ~ J “ ^ 

^ = 1 [$'-4)] = - 


[x-a' (T)-pr (t) + ve^''(j)-vea^\r)l 
[p + vecrl. 


vo! ~ + vd^^' — vda^* = n' (y — 7}) — 

H- c^a — 77 ^'m — <r {a(i — t) — n(y — 77) 4 - w (^ -• S')}, 
where i = /3S' — 7 ^^ m^y^' — n — ar}' — 

Hence we may write 


r d /n\ ^ /m' 

|_3y \v/ 3.2: \ I/, 


4i7r doc \v 


-...(275). 


The electromagnetic field derived from these potentials is due 
to a moving electric doublet. It should be noticed that we 
have the relations 


la-hm^ + ny = 0] 

l^' + mv' + n^' = o\ ^ 

We can write down by analogy the potentials for an electro- 
magnetic field 'due to a moving magnetic doublet. They are 


47r [_dy \ V 


9 f 0Q 
dz [v , 


= - 


?oa„ + ma/3o + no7o = 0 1 


all 


■where Zo^o + mo^o H- ^o 7 o = 0 \ 

Zof + mo77'4-^,or = Oj ^ 

and ofo? A> 7 o> ^o, no are functions of t. When a, y 8 , 7 , m, 
n are functions of r which are not connected by the relations 
(276) the potentials (275) can be used to construct an electro- 
magnetic field which must be regarded as that due to an 
electric doublet and magnetic doublet which move together. 


§ 43. Electromagnetic fields in which singularities are 
projected from a moving point or curve and travel with the 
velocity of light. 

We shall now develop some mathematical analysis of con- 
siderable interest whose physical significance is not yet fully 
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understood. At first sight it seems appropriate for a discussion 
of an emission theory of light in which waves in the aether are 
either produced or guided by small particles which move in 
straight lines with the velocity c. After further study I have 
thought it may be useful for a discussion of the question Has 
the aether a structure ? 

This question has already been raised by Sir Joseph 
Larmor^ and Sir Joseph Thomson f. The latter has, indeed, 
developed a theory in which the aether has a kind of atomic 
structure of which the elements are Faraday tubesj. In the 
most recent form of the theory it is assumed that the electric 
and magnetic forces are zero outside the tubes and that a 
certain amount of work is performed when one corpuscle 
crosses a tube of force attached to another. In an application 
of the present analysis to Sir Joseph Thomson’s theory the 
aim would be to build up his discontinuous electromagnetic 
fields from electromagnetic fields with certain types of sin- 
gularities, making use of discontinuous definite integrals. To 
illustrate the possibility of doing this it will be sufficient to 
mention the definite integral 

~ da ^ ^ 

Jo + ios cos a + ^2/ sin a r 

= 0 z<0. 

The electrostatic field derived from the function V is zero on 
one side of the plane ^ = 0 and has the character of the field 
due to a point charge on the other side. It should be noticed 
that the integrand is a potential function which becomes infinite 
along the line z = 0, x cos a + 3/ sin a = 0, and as a varies this 
line sweeps out the plane of discontinuity of our electrostatic 
field. 

To generalise this result we must endeavour to solve the 

* Aether and Matter (1900), p. 188. The question as to whether the aether 
is continuous or discontinuous is discussed by H. Witte, Ann. d. Phys. (4), 
Bd. 26 (1908). 

t Presidential Address, British Association JRejportSj Winnipeg (1909). 

X Becent Researches on Electricity and Magnetism', Electricity and Matter; 
Phil. Mag. Vol. 19 (1910), p. 301, Oct. — Dec. (1912). See alsoN. E. Campbell, 
The New Quarterly (1909) ; Phil. Mag. Vol. 19 (1910), p. 181. 
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general problem of finding electromagnetic fields whose 
singularities lie on moving curves*. 

A partial solution of this problem may be obtained by con- 
sidering first of all the field represented by equations (10) 
of § 5. We may obtain a suitable pair of functions a, yS by 
solving the equations 


\dxl- 

CD’ 

da d/3 


/day 

[dyj 
(djY 
\dy) 
da 9/3 


.(279); 


da d/3 da 9/3 da 8/3 
dx dx dy dy dz dz 


for clearly 

/da '^_da 9/3' 
VBy dz dz dy 


J. 

dt dt 


.9y dy dzdz ) 


/9aY 

/0ay 

'1/d^V 

fMY' 

/IdadjS 0a 0/3 Y 

\dtj 

WJ 

y[dtj 

[da;/ _ 

St dt dx dxj 


_ 1 rda 3/3 da d^V 

[dx dt dt dx) * 

Two other equations can be obtained in a similar way and so it 
follows that if we make a suitable choice of an ambiguous sign 
which is involved in the definitions of the functions a and 
/3, the equations (10) will be a consequence of equations (279). 

A more general electromagnetic field is obtained by 
multiplying the components of M in equations (10) by an 
arbitrary function f(a,^). Since the components of M are 
necessarily solutions of the wave-equation it follows that, if 
^ = 1//, an expression of the type 

is a solution of the wave-equation (8). 


* The aim must be not only to obtain a coxnplete generalisation of Green's 
equivalent layer which will be applicable to the case of a moving surface, but to* 
obtain, if possible, an analysis of the electric and magnetic current sheets which 
are required. Of. p. 29. 
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The electromagnetic field which has just been obtained 
generally has singularities at space-time points for which 
/(a, /3) is zero. Let us write f—u{x, y, z, t) -h iv (a?, y, z, t) where 
It and V are real when x, y, t are real ; then the points for 
which /= 0 lie on the moving curve, defined by the equations 
u =0, 'y = 0. Now it follows from (279) that / is a solution of 
the equation 



and consequently 



dll dv dit dv dll dv du dv 

dx dx dydy'^ dz dz dt dt ' 


Now let F(u,v) = 0 be the equation of a moving surface which 
always contains the moving curve, then if 


dFdv 

dFdu 

dF dv 

dF du 

du dx 

dv dx' 

^ du dy 

dv dy ’ 

dF dv 

dF du 

dF dv 

dFdu 

du dz 

dv dz ' 

^ du dt 

dv dt ’ 


we have 


mo 


,= (dFY 

1 



Ins) 




VfdFY i 

fdFY 

rdFV 5 =“! 

.fe) ' 

Kdy) + 

\dz 



= 2 


dF dFY f ^F dF\- , / dF dFY 

r* I ( r n — i -I- I r) a - — 1 : 


9^ 


dyl 


( dF 
V 


dz 




dy ^ dx) 


consequently 



This means that the component velocity of the surface F=0 in 
the direction of the normal at (x, y, z, t) is less than the velocity 
, of light, it is equal to the velocity of light only in an excep- 
tional case. Since this is true for any surface that always 
contains the moving curve it follows that the curve can be 
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regarded as moving with a velocity less than that of light. It 
should be understood that the curve generally changes in shape 
as it moves but this is not necessarily the case. 

It often happens that the moving curve /===0 can be 
regarded as made up of the instantaneous positions of a series 
of points which move in straight lines with the velocity of 
light. To see this let us suppose that a function (j> (a, 13) can 
be found such that is a real function of x, y, z, t It is evident 
from (279) that 


dfdcj) 


dfd^^ldfdj^ 
^ 2 - - • 




dx dx cy dy dz dz di dt 

and this means that if a point starts from {x, y, z, t) and moves 

with the velocity of light along a straight line whose direction 

cosines Z, m, n are proportional function 

/ will remain constant along its path, and consequently if the 
point once lies on the moving, curve /== 0 it will always lie on 
this curve. It should be noticed that the function (j> and its 
first derivatives with regard to x, y, z^ t all remain constant along 
the path of the moving point. 

The case in which no such function (j> exists may be oi 
importance in future developments of the subject; this case 
has not yet been discussed in any detail. 

Two methods of solving cejuations (279) are known, bul 
they are not really distinct. In the first method the function 
Of, j3 are defined by ccpiations 


f m + [// ~ ^ (a, m ? («* /9)P - - T (a, m 

I (a, - f (a, ^)] + m (a, /3) [y - r) (a, ^)] 

+ n (a, /3) [z - f (a, fi)] == (a, /3) [Z - r (a, ^)]J 

(280), 


where f P arbitrary functions satisfying th 
relation 

Z® + =» 

The functions ot, 0 may evidently be replaced by two othc 
functions a', defined by equations such as 

consequently we may without loss of generality introduce 
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further relation connecting the functions t], t, r, I, m, n,p. The 
relation we shall choose is 


j 9^ 9i? 9f ^ dr 

9/9 9 / 3 ■^”' 9/3 ~ 


.(281). 


When a and /S are defined in this way we can obtain the simple 
specifications of two types of electromagnetic fields ; one type 
is obtained by writing 



and two similar equations. A second type of electromagnetic 
field may be derived from the potentials* 






_ («, /3) 

FS-QR ’ 

nSg(<x,F) ^_cpSg{a,j3) 

rja ^ rv T-k 


^ _ m8g{a,R) \ 

" F8 - QR 


where 

F 


FS-QR' 


FS-QR J 


..(283), 


^ = H + S <^2/ -^) + H (^- (t-r), 

T> , S’? , 9? „ dr dl , 

-^(2/-’?)-|^(^-r)+o^^ (t-T), 




■T)-|^(aj-a 


dm 


iy- 


doc 

•.)-3^(-?)- 


' dF d/3 

It is easy to see that these potentials satisfy the relation (267). 
To prove that they are wave-functions we remark that 


A.X — 


9{o^, / 9 ) 


FS-QRldF 

dn 


dm . 




+ 


9y8 


n(»-|)-Z(A-01 -c=^{p(«-f)-;(t- 


r8m 9 (a, /3) ^ 9 (a, ff) ^ 9 (a, /3) 

^ _d^'d{x,y) 9/3’9(a;, a) 9/3 ’ 9 («, t) J ’ 


^)}] 


Now it has already been proved that an arbitrary function 
^of a and /3 multiplied by one of the Jacobians represents a 

* We can also call these potentials Xy, X^, A and derive an electro- 
magnetic field from them by the method of § 4. 
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where h, etc. are functions of a which satisfy the equations 

lo^ + -h Uo^ ) 

4- > (288). 

Zo^i + +• 7^o% = d^jpQpi j 

It is easy to verify that the wave-equation is satisfied by 
a function of type 




thus we have a generalisation of the theorem of § 41. 

We shall assume that Iq, rrio, no, po are real and that some 
or all of the quantities ki pi are complex. It is easy 

to see that if we assign a real value to a and a complex value 
to /9 the corresponding space-time points Qc, y, z, t) can be 
regarded as the successive positions of a point which starts 
from the point 

^ = y = i = ^ (289) 

and moves with the velocity of light along a straight line 
through this point. There is clearly just one line through 
this point for each complex value of ^ and vice versa. If we 
consider all the points in space at a particular time t we can 
specify each point uniquely by a real parameter a and a real 
or complex parameter 

Let us now consider the electromagnetic field which is 
specified by the potentials 

A J. -d _”/ dS_2^ 

jp y 


A 


= ^ (290), 


where Z = /3Zi - Zo, m = /9mi - mo, n — /3ni -^no,p^ /Sp^ -po and 
/ is an arbitrary function of a and /3. These potentials are 
derived from (283) by putting Q = 0. 

After a long calculation we find that the component of the 
electric force along the radius from (^, r), a) to {x, y, z, t) is 


pi _ 


To obtain an electromagnetic field in which there is a 
constant electric charge associated with the singularity 
(^, 7), a), we assume that jPo==Pi=/=I 
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da 




drj 

Ta' 


dfx 


4p + mo|5+^„| = (|l^ 


+ 


(S) 


+ 


da) 


...(291). 


da ' "(ja ’ ''"0a vW 
The expression for the radial electric force then becomes 
1 




0 “- 




'V' 


y/2 

5 


Comparing this with the expression for the radial electric 
force in the case of an electromagnetic field with a simple 
sing-ularity (|=, t), r), we see that there is a constant electric 

charge iv-Jc associated with the singularity (^, tj, a). 

I-fc should be mentioned that the second of equations (291) 
is a consequence of the other equations satisfied by 4, Wc Wj. 
To "prove this we take the axis of a; in a direction parallel to 
the velocity of the singularity (^, r), at time a We then 

ha VO for this instant = <5- = 0. If, moreover, we choose the 

oa ca 

axis of y in such a way that = 0, we may satisfy the first 
of equations (291) by writing 

= c cos 6, = c sec 0, nii = ic tan 6, lo = c cos 0, 

mo == 0, Uq = ± c sin 0, 

and then it is clear that the second of equations (291) is 
satisfied. 


Let us 

now 

write 

Zq — c 

cos 6, 

= 

= 0, n^ = 

= + c sin 0, 

li — c sec 6, 

Wi 

= — ic tan 0, 

ni = 0, 

X — 

r=x 

y-v = T, 


t—T 

= T ; then it is easy to see ■ 

that if 


c^T- 

Z,X- 

— miY— 

■ th z, 

S = 

&T~ 

■liX — mi 

Y-riiZ, 

Z7= c^T- 

ZoX- 

-m, F- 

71q Z, 

u= 

&T- 

-loX — Wi 

1 F— rid Z, 

we have 

cos^ 0SS = 

: UU, 

ff+17= 

-2P. 



since yS/S= U it follows that the potentials (290) become 
infinite when U=0, i.e. when 

]/Q tiq 1 

""When a is given these equations are satisfied by a point which 
starts at (^, rj, a) and moves with the velocity of light along 

9 
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a straight line whose direction cosines are proportional to 
hi Wo, rio. This line makes an angle 6 with the direction of 
motion of the point 77, a). 

The electric and magnetic forces in the electromagnetic 
field derived from the potentials (290) are generally infinite 
for ?7 = 0 and so our field possesses a number of singular points 
which are projected from the moving point rj, a) and travel 
along straight lines with the velocity of light. It should be 
remarked, however, that if we retain only the real parts of the 
potentials (290), the projected singularities disappear as soon 
as the singularity 77, a moves in a straight line with uniform 
velocity and I, m, n are independent of a. The field then 
becomes identical with that derived from Li4nard’s potentials. 
To prove this we shall show that, on the above assumptions, 
the field derived from the potentials 




bi: 

P ’ 




m — 77' 


A<> = R , <!>» = R (292), 


is everywhere null, R being used to denote the real part of a 
quantity following it. In the first place we remark that 
we now have 


R = l^' + mij' + c^p == — 


Hence 2R 




o9/l DN 


U 3 , 


SS „ a 




V 0 , V „ 3 


Now since f/'+ P+ 2P = 0, it follows that P/P is a function of 

— 1 . * 

JJjJJ and so Ay^, A^, — - <I)0 are the derivatives of a single 
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function, consequently the electromagnetic field derived from 
these potentials is everywhere null. 

Summing up our results we can say that when the 
conditions (291) are satisfied, the electromagnetic field derived 
from the potentials (290) contains a point charge which moves 
with a velocity less than that of light ; attached to this point 
charge there is a certain curve which becomes the locus of 
a series of moving point singularities whenever its form differs 
in any portion from a straight line or its direction changes. 
The form of the curve at any instant is subject to the condition 
that the points of the curve can be regarded as having been 
projected from the moving charge at different instants, the 
direction of projection being partially determined by the law 
cos 6 = vjc where v is the velocity of the point charge, and 9 is 
the angle between the direction of projection and the direc- 
tion of motion of the point charge. 

We may now obtain a new model of a corpuscle by con- 
sidering an aggregate of elementary fields of the type just 
described, the point charges and exceptional curves being 
nearly coincident. If we write de for the charge associated 
with one of the elementary fields we may obtain a field in which 
the electric and magnetic forces are finite by a suitable process 
of integration. According to this idea a corpuscle has a kind 
of tube or thread attached to it. When the motion of the 
corpuscle changes a wave or kink runs along the thread; the 
energy radiated from the corpuscle spreads out in all directions 
but is concentrated round the thread so that the thread acts as 
a guiding wire. This theory of radiation is in some respects 
similar to that given by Sir Joseph Thomson in his theory 
of the Rontgen rays*. It is in accordance with his idea that 
the energy may he concentrated round certain points of the 
wave-front. 

The following figure indicates roughly the changes in 
the form of a tube which always lies in one plane and is 
attached to a corpuscle performing a simple harmonic motion ; 
^it is seen that a type of progressive wave travels along the tube. 

* Electricity and Matter, London (1904) ; Phil. Mag. Vol. 19 (1910). 

9—2 
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In Sir Joseph Thomson's theory of the Rontgen rays the 
kink in the tube of force becomes longer and longer as it 
recedes from the charge. A similar remark applies to the 



oscillations of the thread attached to our point charge. This 
phenomenon may be due entirely to the fact that the tube of 
force and thread extend to infinity. If we suppose that the 
tube or thread does not extend to infinity but ends at some 
other point charge, the circumstances of the motion will be 
different4 If in this case we treat the thread as a singular 
line of an electromagnetic field and suppose that it is given by 
an equation of the form 

/(«,/3) = 0 

where a and yS are functions which satisfy (279), we must 
conclude that there is no function of type F {a, /3) which is 
a real function of for if this were the case the moving 

thread would be the locus of points travelling in straight lines 
with the velocity of light and would consequently extend to 
infinity. 

Electromagnetic fields with moving point charges joined by 
singular curves which do not extend to infinity have not yet 
been obtained, but I think there is some hope of deriving them 
by the general methods of § 43, when the quantities a, /3 are 
both complex. 

EccamjpUs. 1. Discuss tke properties of the electromagnetic fields 
that can be derived from the potentials 


'*•-7+?' 

mi mi 

s ’ 



s s ' 


* is is ^ 

A r 

^ “* is is ^ 

a; 

7^1 ni 

is is ’ 

is is ' 



respectively and determine the lines of electric and magnetic force. 
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2. Prove that if a and 0 are defined as in § 44, the electromagnetic 
field specified by the equations (282) is conjugate to the field specified by 
the potentials (290), or to the field specified by Lienard’s potentials. 

§ 45. Electromagnetic fields with singularities moving 
with velocities greater than that of light. 

Some of the preceding analysis holds and provides ns with 
solutions of Maxwell’s equations when the velocity of the 
primary singularity (^, rj, r) is greater than that of light, but 
in the case of a field specified by potentials of type (266) 
a transition from a velocity less than that of light to a velocity 
greater than that of light does not seem, to be physically pos- 
sible on account of the occurrence of infinite values of the 
electric and magnetic forces in the critical case. Moreover, it 
is difficult in the general case to give a rule which will enable 
us to pick out just one root of the equation (265). An interest- 
ing type of field may, however, be obtained by a process of 
summation over some of the roots of the equation*. 

The case of infinite velocity is of some interest, for then 
we obtain electromagnetic fields with singularities along a 
fixed curve at a given instant of time. The following example 
indicates that the case in which the primary singularities are 
imaginary may be associated with another case in which they 
are real. 

Consider the two equations 

{x—a cos aY + (y — a sin oCf 
^ q- _ ai cosh /3)2 = {ot — id sinh /9)^ 
and write in analogy with (263) 

j/ = a cos a{y — a sin a)— a sin a{x — a cos a)= a {y cos sin a), 

j/j = ia sinh /3 (^ — ia cosh — mcosh /3 {ct — ia sinh /3) 

= m {z sinh ^ — cosh /S). 

We evidently have 

4i^^ == 4a2 (aP' + y^) — (.-r^ y^ + z^ — 

W = (z^ - cH^) + (^’" + y" + - cH^ - a% 

hence + Vi = 0. 

* A more complete discussion is given in G. A. Schott^s Electromagnetic 
Radiation. 
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It is easy to verify that the function 

“ = ;/(«> -5) 

is a wave- function and that the functions a, ^ are suitable for 
obtaining electromagnetic fields by the method of § 6 ; they are, 
in fact, the functions considered in § 36. 


§ 46. Second solution of the fundamental equations. 

The fundamental equations (279) are also satisfied when 
the functions (a, are defined by the relations 

(a, /9) y -f- n (a, {a, /3) t + g(a,^)=: 0| 
L(a,/3)x’^M (a, j3) y (a, ^)z ~ d^P {a, ^)t+Q {a, IB) = oi 

(293) 


where I, i, etc. are arbitrary functions satisfying the relations 

D + ^2 

IL +mM + nN = (294). 

If now 

^ dl dm dn , dg 

dl 9m , dn . , 9p dg 


dL dM dN 


SP dO 


V=a;3- + W — + > 

9a 9a 9a 9o 9a 


dL 

'9/3' 


, dM dN dO 

‘*9/3’^^ 9/3 ■^^9/3 ®*9/9‘^ 9y8’ 

we can derive an electromagnetic field from the potentials 


A-Jiff- A A--^^ 

XL jjj — j XL y ^ f Z ”” -V ? — "x 

AtiT — /MV Aor — /IV Act — /MV Act — /mv 


<E> = : ■' 


where f is an arbitrary function of a and /3, and 
, 9P , 9ilf , dM , 9P 


.(295) 


( 296 ). 
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To verify tliat these potentials satisfy the wave-equation and 
the relation (267), it is sufl&cient to remark that 




3/3 3 {x, y) ■ 


dN d (a, j3) 

3/0 d (Xj z) 


aP 3(a,/3) 

3/3 3 {x, t) 


_ idM d(a,$) idN dja,^) . 3P 3(a,y0) 
c d/3 d (z, t) ”<3 3/3 3 (y, t) ”^^3/3 3 (y, -s:) 


L..(297). 


When similar expressions are obtained for Ay, Az, ^ it is clear 
that the relation (267) is satisfied. Other types of electro- 
magnetic fields are obtained by writing a instead of /3 in 
equation (296) or by writing I, m, n, p in place of L, M, N, P 
in (295) and (296). If AJ, Ay, Az, <!>' are the potentials 
obtained in either of these ways we have clearly 

A^A^ -f- Ay Ay' ^AzAz'-<i>^'^0 (298). 

I have noticed that this relation is often satisfied by the 
potentials of two conjugate fields. 

Another type of electromagnetic field may be obtained 
by the method of § 5. 

Some particular cases of the preceding theorems may be 
deduced by contour integration. To illustrate the method let 
us suppose that l,m , ... X, AT,... are functions of a parameter a 
which satisfy the relations (294). If we regard these quantities 
as independent of x, y, z, t, the contour integral 

1 ^F{Lx’\-My-\-Nz — &Pt-\-G,a)^^ 

^Tvi J lx + my -\-nZ'- &pt -H g 
will represent a wave-function. Now let us suppose that the 
contour encloses only one root of the equation 

xl (a) + ym (a) -f zn (a) — cHp (a) -hg(a) = 0 , ..(299), 
and that the numerator is finite and single-valued within the 
contour and on its boundary, then by Cauchy’s theorem the 
value of the integral is generally 




(300) 


where /3 = (a) + yM (a) -h zJV (a) — (a) -f G (ce). . .(301) 

and X is defined in the same way as before. 

We have then the result that the function (300) satisfies 
the wave-equation. This is a particular case of the general 
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theorem of § 43 and is a generalisation of a theorem due to 
Forsyth*. In Forsyth’s work the functions L, My Ny P are 
assumed to be the derivatives of ly m, iiy p with regard to a. 
The functions a, jS are evidently particular cases of the func- 
tions that have already been defined and so may be used to 
construct an electromagnetic field by the method of § 5. An 
interesting electromagnetic field may also be derived from the 
potentials 




X 


Ay — 


m' {a) 


. («) 


X 


(302); 


it is easy to verify that the relation (267) is satisfied. The 
case in which Z, m, riy p, g are all real functions of a is un- 
interesting because then our potentials become infinite for 
00 ^ space-time points which lie in oo ^ planes. When, however, 
ly m, Thy p, g are complex functions of the type ^ (a) -h (a), 
the singularities of the electromagnetic field generally lie on 
a moving curve. 

It should be remarked that when a and /? are defined by 
the equations (299) and (301) a function of the type 

ft = F (a, /3) 


satisfies the wave-equation, and at the same time satisfies the 
differential equation 


/aft Y /aft Y /^y = 1 

\dx) Vay/ \dz) c^\dt) ' 


This is a generalisation of a theorem due to Forsyth f and 
JacobiJ. A solution of MaxwelFs equations may be derived 
from the potentials 

A^^^k{a)y A/ = mo(a), -4/ =%(«), <I>® = cpo(a) ...(303) 


if 


^aa' 


drriQ , a^o 
m + n 
da da 




■ 0 . 


k 


di 


dm 

da 


dn 

^Ta 


dp 

,0 f = 0. 


’aa' 


* Meu&nger of Mathematics y Vol. 27 (1898), p. 188. The theorem is 
obtained in another manner by Forsyth, 
t Loc. cit. 

fFerkCy Bd. 2, p. 208. 
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One way of satisfying these equations is to put 
Zo = I, 77io = m, 7io = n, po = p. 

An interesting type of wave-function may be obtained by 
a generalisation of a method due to Schottky 
Let 

Qq (a)=jl (a) da, bo (a) =/ m (a) da, Cq ( a) —Jn (a) da, 
do (a) =Sp (a) da, eo (a) = Jg (a) da, 
then il = xao (a) + ybo (a) + zCq (a) — cHdo (a) -h Bq (a) 
is a wave-function whose derivatives with regard to x, y, z, t 
are all functions of the single variable a. The function a is 
supposed to be defined by equation (299). 


Example, If ^ = 1 — w = s + a, n-=^i(s — a), cp=^l-\-sa, g= a^, where s 
is a constant, the function X vanishes when {x, y, z, t) lies on the moving 
curve 


l-s2 




; o2 _ m2 


5 — ^2) — 2s 


Z — ’^'V, 


A point for which v is constant moves in a straight line with the 
velocity of light. 


§ 47. A wave-function witli a fixed curve of singu- 
larities. 

Let T be defined in terms of x, y, z by the equation 
ar = xp (t) + yq (r) -f- zr (r) 
where p® -h + r® = 0 and a is a constant. 

Let 1 / = a — xp' (t) — yq' (r) — zr (r), 

e^ot± 

{p'^+q'^ + r'^y 
then is ^ wave-function "f. 

A solution of Maxwell’s equations may be derived from the 
potentials 

Ax=^/{d,r), Ay-^f{6,r), Az-'^f{0,T), €> = 0 . 

* Berlm. Sitzungsberichte (1909), p. 1152. 

t This result was derived from a theorem given by Prof. Porsyth, loc. eit. 
If we put r=a the functions a, jS can be used to obtain an electromagnetic 
field by the method of § 5- These functions are, of course, particular cases of 
the functions defined at the beginning of § 46. 
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To obtain a real electromagnetic field we must retain only 
the real parts of these expressions. 

It can be shown by putting 

jp = Pi (zc, v) 4- iPa (% v\ q = Qi (% z)) + iQa 

(u, v) + iPo (Uy v), a— (u, v) 4- iA^ (u, v), 

that the electromagnetic field has generally a fixed curve of 
singularities. In the special case when 

g=2T, r=:i(l + P), Ai=:hy A^ — h 
the fixed curve is the circle ^ 4- J y = |- A. 

§ 48. Cylindrical wave-functions with moving singu- 
larities. 

If we define r in terms of Xy y by the equation 
[x - ^ (t)]2 4 - [y - '7 (t)]2 = t)2 

and define I (r), m (r), p (r), so that 

(t) + mr}' (r) = c^j), 

P 4- 

the function v = Z 4- m (y — 77) — (^ — t) 

is such that ^f(P) 

is a wave-function. 

In particular, if ^ = 77 = 0, Z = 1 , m=i, p = 0 we obtain 
Poisson’s wave-function 



Another interesting result is that if F {Xy y, t) is a homo- 
geneous function of degree ^ which satisfies the wave-equation, 
the function 

n = P(a?-^, y-77, t-r) 
also satisfies the wave-equation. 

If 0 - = f (T)(a;-|) + i?'(T)(y-i;)-c“(«-T), 

the function 


1 
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satisfies the wave-equation only when f '^ + 17'^ = When 

iz is written in place of ct these wave-functions can all he 
regarded as solutions of Laplace’s equation. 


EXAMPLES. 


1 . Let 


z=xf(r) - (s) $ («) ds, 2 = 2 / (r) - (*) C W 

y=yf{-^) - («) V (s) t= </(r) - J^/' W sds ; 

where r is defined in terms of y, z, t by means of equation (265) and the 
inequality Prove that if an electromagnetic field H) is such 

that 


c di \c dt 




with similar equations, where M=^H±iE, then an electromagnetic field in 
the variables i*, y, t can be found such that we have identically 
M^d (y, z)-\rMyd {z^ x) + M^^d {x, y) + icM^d {x, t) 

Hh ic My d (y, t) + ic M^d (2, t) 

s Mxd (y, z) + Myd {z, x) -\-Mgjd {x, y) ^ 

+ icMxd (x, t) + icMyd (y, t) + icMgd (z, t), 
where d(y, z) denotes dyhz-dz^y and dx^ hx^ etc. are two independent 
sets of increments of the variables. 

2. Prove that if /'(r) is always positive the variable t increases with t. 

Show also that if (J, J, f) correspond to (^, 7, r), the point (|, 7, C) 

moves along a curve r with a velocity less than that of light and that the 
velocities at two corresponding points of F and f are the same in magnitude 
and direction. 

3. Prove that the conditions imposed upon the electromagnetic field 
in Example 1 are all satisfied in the case of the field specified by the 
potentials (266). Hence show that the transformation transforms the 
field of an electron moving along the curve F into the field of an electron 
moving along the curve F. 

4. Prove that the conditions of Example 1 are also satisfied for any 
electromagnetic field of type (282) where a and /8 are defined as in § 44. 
Hence show that a field of this type is transformed into another field 
of the same type associated with the curve F. 
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5. Prove that the conditions of Example 1 are equivalent to only two 
conditions which imply that the field {E, H) is conjugate to any electro- 
magnetic field of type (282). 

6. Prove that electromagnetic fields of types (302) and (303) can be 
transformed into fields of the same general types with the aid of the 
transformation 

I (s) f (s) ds, z=z+J^ n (s) f (s) dsj 

m{s)f(s)ds, i=t+ p («)/(«) ds. 

7. If r be defined in terms of x, y, t by the equations 

ft r)]2-|-[y-97(a, ft t)]^ + [2 - (a, r)P = ^ (a, ft r)P, 

it satisfies the partial differential equation 



8, If 
and 


A'—p cos 


y=psin0, 

. , , — ^2 


the function +y ^ + (2 ± a)^] [0 + 

satisfies Laplace’s equation. 


9- Prove that if in the last example we write 

f[^ + ^ 

and integrate with regard to u between 0 and 2ir we can obtain a potential 
fimction which is zero outside the tube \ <l> + iyf/\^L 


10. Particles are projected in certain directions from the different 
positions of a moving electron and travel along straight lines with the 
velocity of light. Prove that if the law, according to which the direction 
of projection varies with the velocity of the electron, be suitably chosen 
the particles will at each instant form a line of electric force in the 
electromagnetic field due to the moving electron. 


CHAPTEE IX 


MISCELLANEOUS THEOEIES 

§ 49. KirdihoflTs formula and its extensions. 

An. important solution of the wave-equation is embodied in 
KirchhofFs formula* which is usually interpreted as the mathe- 
matical expression of the principle of Huygens. This formula 
has been extended by Lovef and MacdonaldJ so as to give 
a representation of .an electromagnetic field outside a surface 
in terms of the electric and magnetic forces tangential to 
the surface. In Macdonald’s formula it is the time derivatives 
of E and H that are so expressed. Tonolo§ has given a 
formula in which E and H are expressed in terms of their 
surface values. The formulae are given in examples 3 — 5 at 
the end of this chapter. 

When the sui'face is a sphere Kirchhoff s formula reduces to 
the formula of Poisson || (Ex. 5) which enables us to find a wave- 
function which satisfies the conditions 

« =/(«> y> ^). y> 

Poisson’s formula may be used to derive the theorem^ that the 
mean value of a wave-function w over a sphere of radius cr at 
time t is equal to the mean value of at the centre of the 

* Berlin. Bet. (1882), p. 641; Wied. Ann.. Bd. 18 (1883) ; Ges. Ahh. t. 2, 
p. 22. Simple proofs of the formula have been given by Beltrami, Rend. Ace, 
Line. Rom. (5), t. 4 (1895) ; Larmor, Proc. London Math. Soc. Ser. 2, Yol. 1, p. 1 ; 
Love, Ihid. p. 37 (1903); Lamb, Hydrodynamics, 2nd edition (1906), p. 477; 
H. A. Lorentz, The Theory of Electrons, p. 233 ; E. Laura, II Nuovo Cimento 
(1913). 

t Phil. Trans. A, Yol. 197 (1901). 

X 'Electric Waves, p- 16; Proc. London Math. Soc. Ser. 2, Yol. 10 (1911), 
p. 91; Phil. Tram. A, Yol. 212 (1912), p. 295. This theorem gives an 
analytical specification of a generalised Green’s equivalent layer. See p. 29. 

§ Annali di Matematica, Ser, 3, t. 17 (1910). 

11 The details of the calculation are given by Love, loc. cit. A simple proof 
of Poisson’s formula is given by Lamb, loc. cit. p. 471. 

H Cf. Bayleigh’s Sound, appendix, and H. Bateman, Amer. Journ. (1912), 
where some other theorems of a similar kind are given. 
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sphere during the interval ^ — r to t + r. The function u is 
subject to the conditions in Kirchhoff’s theorem. 

When the function ii is independent of Poisson’s 
formula reduces to ParsevaFs formula for a cylindrical wave- 
function. Volterra^ has extended ParsevaFs formula so as 
to obtain a two-dimensional analogue of Kirchhoff’s formula. 
His formula indicates that the propagation of cylindrical 
waves is essentially different in character from that of 
spherical waves. In the three-dimensional case the value of 
a wave-function u {x, y, z, t) at a point (x, y, z) at time t 

is completely determined by the values of u and — over a 

concentric sphere of radius or at time t — r. In the two- 
dimensional case, on the other hand, the value of u {x, y, t) at 
a point (x, y) at time t is not determined by the values of u and 

^ over a concentric circle at time t~-r. To find u {x, y, t) we 
ot 

must know the values of u and over a series of such circles 

ot 

in which the radius cr varies from 0 to some other value ctj. 
The essential diflference between the two cases may be attri- 
buted to the fact that in the three-dimensional case the wave- 

function for a source is of type while in the two- 

dimensional case it is of type ^ cosh and a wave 

does not leave the region undisturbed after it has passed, but 
has a tail or residue f. 

When u is a periodic function of t, Kirchhoff’s formula may 
be replaced by the simple formula of HelmholtzJ. In this case 
there is an analogous formula for cylindrical wave-functions, the 
function Kq (ipk) taking the place of 


§ 50. Greenes Functions. 

The solution of a problem in which a periodic wave-function 
is to be determined from a knowledge of its behaviour at 

* Acta Math. 1. 18; Lectures at Clark University (1912), p. 88. 
t See Lamb’s Hydrodynamics j p. 474. 

X See also J. Hadamard, Bull, de la SocUU math, de France, t. 28 (1900), 
p. 69 ; J. Larmor, Proc. London Math. Soc. (2), Vol. 1 (1903), p. 13. 
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certain boundaries can be made to depend on that of an 
auxiliary problem, viz. the determination of the Green’s 
function*. 

Let G (x, y,z\ x^.y^^ be a solution of Aw + = 0 with 

the following properties: It is to be finite and continuous, 
as also its first and second derivatives, in a region bounded 
by a surface /S, except in the neighbourhood of the point 
where it is to be infinite like cos foy47rr, when 
r 0. At the surface >?, G satisfies some boundary condition 

such as (1) ^4 = 0 or (2) ~ = 0. 

Adopting the notation of Plemelj f and KneserJ we shall 
denote the values of a function (^, r), f) at the points (x, y, z\ 
(^’i j 2/i, <2^0 respectively by ^ (0) and (1). The Green’s function 
is then denoted by the symbol (t(0, 1). The importance of the 
Green’s function depends chiefly on the following theorem. 

Let be a solution of 

+ k-<l> +f{x, y, z) = 0, 

which is finite and continuous, together with its first and second 
derivatives, through the interior of the region and satisfies the 
same boundary condition as Q (0, 1), then 


0 (1) = [[[/(0)G (0, 1) dxdydz. 
J JJ 


This theorem is proved by applying Green’s theorem to the 
region between a small sphere 2, whose centre is at (xi,yi,Zi), 
and the surface S. For since 


jJI (<j)AQ - GA<p) dxdydz =lf(^4>^-G dS 




we obtain the required relation by making 2 ->• 0 and using the 
boundary conditions. 

* This function was first used by Green in the solution of a problem of 
electrostatics, Essay on the application of mathematical analysis to the theories 
of electricity and magnetism, Nottingham (1828) ; Math. Papers, p. 31. 

t MonatsJiefte fur Math. u. Phys. (1904) and (1907). 

X Die Integralgleichuyigen und ihre Anwendung iri der mathematischen Physik, 
§ 31, Brunswick (1911). 


144 


MISCELLANEOUS THEORIES 


[CH. 


If ^(2,0) is the Green's function for the same boundary 
condition as G (2, 0) but for k = cr, we must also surround the 
point (2) by a small sphere when we apply Green’s theorem 
with 0(0) = ^(2,0). We then obtain the equation 

g {% 1) = © (2, 1) - - <r=) (2, 0) (0, 1) dcodydz. 

This important relation indicates that 5 ^( 2 ,!) is the solving 
function of an integral equation of which (? (0, 1) is kernel and 
mce versd. The theory of integral equations tells us that when 
(0, 1) is given there may be certain singular values of a-- 
for which g {2, 1) is not finite. These are the values of cr^ for 
which the homogeneous integral equation 

^ (1) = {a- -k^)l jj ^ (0) 0 (0, 1) dxdy dz 

possesses a continuous solution <f> (0) which is different from 
zero. Formula (2) indicates that for such values of cr^ the 
differential equation Ac^) = 0 possesses a solution satisfy- 
ing the boundary condition and the other conditions imposed 
on <j). The solutions of this type are of great importance in the 
theory of sound and have been discussed by many writers*. 

If we put f{0) - — k^) g (0,2) and proceed as before, 

Green’s theorem gives 

g (1, 2) = (? (2. 1) - j g(0, 2) G (0, 1) dxdydz. 

Putting cr = k and comparing this with the previous equation 
we get 

^( 1 , 2 ) = ^( 2 , 1 ). 

Hence the Green’s function is a symmetric function of the 
coordinates of the points 1, 2. When the boundary condition is 

^ = 0 this result is equivalent to Helmholtz’s theorem f. 

Since G is a real symmetric function when A? = 0 it follows 
from the general theory of integral equations that there is 


* See especially Lord Eayleigli , Theory of Sound, Vol. 2 ; Pockels , Die partielle 
Differ entialgleichung Au + k’hi=i0 ; A. Sommerfeld, Encykloyddie der Math^ 
Wiss. Band n. 1, Heft 4, p. 540. 

t Cf. Rayleigh’s Sound, Vol. 2, p. 131. 
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at least one real singular value of ; that all the singular 
values are positive may be deduced at once from the. 
equation 



dxdydz 



The Green’s function is usually obtained in practice by finding 
a suitable expansion in terms of elementary solutions of the 
equation A?,/; + = 0. This method is explained in Heine’s 

Kugelfunktionen and many examples of Green’s functions are 
given for the case k = 0. The general case has been discussed 
at length by A. Sommerfeld^ who also obtains a number of 
definite integrals which represent Green’s functions. These 
expressions lead to interesting generalisations of Fourier’s 
theorem. 

The problem of electrical oscillations in a cavity has been 
discussed by Weylf. With the aid of a generalisation of the 
Green’s function, viz. a Green’s tensor, he obtains a number 
of inequalities satisfied by the periods of vibration. 

The Green’s function for the equation can 

theoretically be found when the corresponding Green’s function 
for the equation == 0 is known. Considerable progress has 
been made in the theory since the appearance of Heine’s work 
and so a few references to recent literature will be useful J. A 


* P%«. Zeitschr, Bd. 11 (1910), p. 1087; Jahresbericht der deutsch. math. 
Verein, Bd. 21 (1913). 

+ Math. Ann. Bd. 71, p. 441 ; Crelle, Bd, 141 (1912). 

X For tlie determination of special Green’s functions see E. W. Hobson, 

Cambr. Phil. Trans. Vol. 18 (1899), p. 277; H. M. Macdonald, Ibid. p. 292, 

Proc. London Math. Soc. Vol. 26 (1895), p. 161 ; A. G. Greenhill, Proc. Cambr. 

Phil. Soc. Vol. 3 (1880) ; J. Bougall, Proc. Edinburgh Math. Soc. (1900) ; 

H. S. Carslaw, Ibid. (1912), Proc. London Math. Soc. (2), Vol. 8, p. 365; C. W. 
Oseen, Arkiv for mat. Bd. 2; C. Neumann, Leipziger Berichte, Bd. 58 (1906), 
Bd. 62 (1910) ; W. Burnside, Proc. London Math. Soc. Vol. 25 (1894), p. 94. 
For the general theory H. Poincard, Palermo^ t. 8 (1894), p. 57; 

S. Zaremba, Ibid. t. 19 (1905) ; E. B. Neumann, Studim vber die Metliode 
von C. Neumami nnd G. Bohin zur Losung der beiden dwertaufgahen der 
Potentialtheorie^ Leipzig (1905) ; D. Hilbert, Gott. Nachr. 1904) ; M. Mason, 
Newhaven Math. Colloquium (1910) ; E. Picard, An7i. de VJ&cole Normale (1906), 
p. 509. 


B. 
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good account of the developments up to 1900 is given in 
Sommerfeld’s article in the Encyhlopddie der Mathematischen 
Wissenschaften, 

§ 51. The transformation of the electromagnetic 
equations. 

The transformations which can be used to transform any 
solution of the wave-equation into another solution or any 
electromagnetic field into another belong to a group which 
is characterised by a relation of the form* 

daf^ -f dy''^ 4- — c^dt'^ (dx^ + dy^ + dz^ — c^dt^). 

The linear transformations belonging to this group are of great 
importance in the modern theory of relativity f ; two of the 
non-linear transformations have been mentioned in § 13. 

In addition to these transformations there are other trans- 
formations, involving arbitrary functions in their specification, 
which can be applied to certain types of wave-functions, and to 
certain types of electromagnetic fields. There are often two 
families of wave-functions to which a given transformation can 
be applied, when the transformation is of a suitable character ; 
each of these families may be defined by a linear relation which 
exists between the wave-function and its derivatives, sometimes 
between the derivatives alone. Some idea of the theory may be 
derived from the examples. It also happens that there is often 
a family of electromagnetic fields to which a given transforma- 
tion can be applied and this family is defined by means of two 
linear relations between E and which can be interpreted to 
mean that the field is conjugate to some definite electromagnetic 
field or family of electromagnetic fields determined by the 
transformation. In some cases these last fields are self-con- 
jugate and the transformation is applicable to them also. 

* H. Bateman, Proc. London Math. Soc. Ser. 2, Vol. 7 (1909), Vol. 8 (1910) ; 
E. Cunningham, Ibid. Vol. 8 (1910). 

t For this see A. Einstein, An7i. d. Phys. Bd. 17 (1905) ; Lane, Das 
Relativitdtsprinzip, Brnnswich (1911); E. Cunningham, British Association 
Reports (1911) ; H. Minkowski, Gott. Nachr. (1908) ; E. B. Wilson and G. N. 
Lewis, Proc.Amer. Acad, of Arts and Sciences, Vol. 48 (1912), p. 389 ; J. Ishiwara, 
‘‘Bericht iiber die Relativitatstheorie,’’ Jahrbuck der Radioaktivitat und • 
ElektronikfBd. 9 (1912), pp. 560 — 648; L. Silberstein, The Theory of Relativity , 
Macmillan and Co. (19i4). 
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The fact that the condition of conjugacy between two 
electromagnetic fields often implies the existence of one or 
more transformations depending on arbitrary functions, may 
be regarded as of some philosophical interest. 


MISCELLANEOUS EXAMPLES. 


1. Show that the most general periodic solution (valid for all 
space outside a given closed surface) of the wave-equation is 

where P and Q are arbitrary functions, r is the distance from the element 
of surface djS to the point where ^ is estimated, and is the angle between 
r and the outward drawn normal. Show further that the necessary and 
sujBficient condition that the value of given by the same analytical 

00 

expression, should vanish for points inside the surface, is that P= ^ . 

(Cambr. Math. Tripos, Part II, 1904.) 

2. Let Q be a function which satisfies the wave-equation and is such 
that its differential coefficients of the first order are continuous functions 
of y, t within a region bounded by a closed surface If either Q or 

— be given for points on the surface S there is only one function Q which 

reduces to a given function /(^, y, z) for t—tQ. 

(A. E. H. Love, Proc. London Math. Soc. Ser. 2, Vol. 1, p. 42 ; 
J. Hadamard, Bull, de la Soc. Math, de France.^ t. 28 (1900).) 


3. If throughout a specified region of space and a specified interval of 
time % and its differential coefficients of the first order are continuous 
functions of a', y, z and of if also the differential coefficients of the 

second order such as ^ are finite and integrable, then a solution of 


the equation 

Qu+cr (. 2 ?, y, Zj t)=0 

which is valid for this region is given by the formula 




r^' 

dn 


1 dr pw'll 
r dn LSUJ 


dS 


where r2=(^-Xo)^ + (y-"yo)^+(^""- 2 ro) 2 , n denotes the normal to dS drawn 
into the specified region and the integration is taken throughout this 

10—2 
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region and over its boundary. The function a- is supposed to be finite and 
integrable and a quantity within square brackets is calculated at tinio 

(G- Kirchhoffi) 

4. If an electromagnetic field is such that the specified region docs 
not contain any charges or convection currents and the value 

of the vector if at (.a;o, 3 / 0 » 2 ^ 0 ) is given by the formula 


47ri/; 


I cos [i/^] cos r:^ — [if J cos rn} ~ 


/ fm cos nx - [if^] cos r/v ~ [i4] cos ni 

(A. Tonolo.) 


C dtn 


[if J cos — t [if y] cos 


6. Prove that in the same circumstances 


47r 


0 

3^0 c 

+t 


klMd) 


000 


0^0 9yo 


dS 


- 0^0000 \rj dt(? \r, 




where a==/jL[ifJ — v [if,,], etc. and (X, jx,v) are the direction cosines of the 
normal drawn into the region bounded by S. 

(H. M. Macdonald.) 

6 . If u is a wave-function independent of 0 and periodic in t like 

M(«o.yo)=^ J 1^0 W ^ Jro(iir)jdt. 

7. A wave-function which satisfies the conditions 

«=/(«’, y. 2) *) 

is given by the formula 

(<7)+<y, 

where /, ^ denote the mean values of /, ^ respectively over the surface of 
a sphere of radius cjf having the point y, 0 as centre. 

(S. I). Poisson.) 

8 . If u satisfies^ ^ ~ 0^2 finite second derivatives within 

a suitable domain 

1 0 r . 

« (^x, yx, 0 = 2 ;; 


where 




p’‘=(x-Xiy+(y-yi)<‘, 
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and a denotes the area within the circle cut out on the plane T= t by the 
cylinder 

(Parseval and Yolterra.) 

9. Prove that if a transformation of variables from {x, y, z, t) to 
(^j y, i, t) is such that 

(& 4- dy^ — (dx^ 4- — c^dt^^) 4- {Idx -^-mdy-k^ndz — cp dt) 

{l^dx -^-mQdy ’\-n^dz- GpQdt\ 

where + it can he used to transform ^ 

electromagnetic field {E^ R) into another electromagnetic field R) 
with an identical relation of the same type as that used in Ex. 1, Ch. vm, 
if the two conditions embodied in the relation 

{mn^ — m(^n) 4- My (nlo - nc^l) 4 - M'^ (ImQ - lom) + ^ (Ipo - kp) 

+ iMy (mpo - -mop) + i {npo n^p ) = 0 

are satisfied. 

Prove that the conditions can also be thrown into the form 
{Uq 4- mmo 4- oiuq 4-ppo) M^ + i (^Po + n^p) My ± ^ (mpo 4- 'rri^p) M^ 

—Iq {lMx’^'^My-\-nM^-\-l 

and similar equations. 

10. In the last example if ^o=^j Wo==n., po= —p, conditions 
are satisfied if Poynting’s vector is in the direction m, n). Show, in 
particular, that the transformation 

z=f{r+ct)+g{r-ct\ ct=f{r+ct)-g{r-ct) 

can be applied to an electromagnetic field in which Poynting^s vector is 
along the radius from the origin and that in the resulting electromagnetic 
field Poynting’s vector is parallel to the axis of z. Apply the transformation 
to the electromagnetic field derived from the functions a, 3 given by 
equations (13), § 5. 

11. If a transformation of coordinates is such that 

ds^ 4 - dy^ 4 - d^ — c^di^ = X^ (dx^^ 4- dy'^ 4 - dz^ ~ 
where X is a function of x, y^z,t\ there is also a relation of type 
+ (y - yo)^ + (2 - - c2 

= XXo [(^ - + (y - yo)^ H- (2 - <o)’']- 

(J. Liouville and S. Lie.) 

12. Prove that the differential equation 



is covariant for a transformation of the type considered in Ex. 11. 

(S. Lie.) 



150 


MISCELLANEOUS THEORIES 


[CH. 


13. Prove that if g, rj, r are functions of o!, y, z, ( such that 
± icrM^ dx + - fi/a- ± icrMy) dy 

4- ir]My. - ^My ± ictM^) dz + ic {^M^-^-riMy + dt 
is an exact differential and c is a quantity whose square may be neglected, 
the value of Jf at the point y'=y+^»?j t' — may 

be calculated by assuming that the integral form 

Mx {dydz — dz by) + My {dz bx — dx bz) + M^ {dx by - dy bx) 

+ icMx {d>x bt ~ dt bx) + icMy {dy bt - dt by) + icM^ {dz bt - dt bz) 
is an invariant for the infinitesimal transformation, it being supposed 
that the function M satisfies equations (10) of § 5. 


14. Let a transformation from the coordinates {x\ y', f) to {x, y, e, t) 
be such that 

dx^^ + dy'^^dz^-^c^dt^^^s/l {dx^+dy^-¥dz^--c^df) 

0' 

4 — -j=:{ldx-\-mdy -\~ndz— c^pdty 

V A 

2^ 

— j=z (Jtdx^mdy-k-ndz — e^p dt) {I'dx -f m*dy 4- n'dz ~ c^p*dt) 

V A 


4 — r=i{Vdx-\-vn!dy 4- ri!dz - c^pdt)\ 
V4 


where 


e == 4* — <^p\ O' ^ 4" n'^ c^p'\ 

^ = cr 4- 4- wtwi' 4- nn' - c^pp\ A = #2 _ 00 * . 

and or, I, m, n, l\ m', p' are functions of Xy y, Zyt^ then if 6 satisfies 
the equations 

,0^ b6.ddt d6j^ 

^0i+”*§5+’^a^+-?’07+^^=®’ 

0 r 




it also satisfies 


W ^ 02^ 1 W 


In the preceding equation we have 
_dR dQ bJL 
0y 02 : 0^ * 

Jbq 'bP W 

dx dy dt ’ 


d^ ^ m 

dz dx'^ dt ^ 
_dZ d_M dJSr 
^ dx dy 02 ’ 


P— vmi T- m!%y R^lm! - Vm, 

L^lp' — I'py mp* '^x^py N^np' — n'p. 

Show, that in certain cases a function of type is a solution of the 
wave-esquation in consequence of the two equations imposed on Discuss 
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the case of the transformation of Ex. 1, p. 139 : also the case of a trans- 
formation which leaves the functions X, Y, Z of Ex. 7, p. 80 unaltered in 
form. 

15. Prove that if a function F satisfies the equation 

' l?!Z 

02/2 02f2 » 

vanishes at infinity and has continuous derivatives except at points of the 
curve or where its normal derivative is discontinuous in such a way that 

the symbol Q being used to denote the coordinates of a point Q, then 

where 

(Levi-Civita, Mtovo Cimento, 1897.) 

16. Prove that 

1 


VCar - + (y -vY + (2 - f 

1 _[” f“ c” 

27r^J _Qo J -00 J ^ 




- dkdfidvj 


and deduce that the integral j d^dy)d( satisfies Poisson’s 


equation A F+ 47rp (^, y, 2 ) =*0. 


(J. Weingarten.) 


17. Prove that the equation A(i) = ^ + 2 ™ is satisfied by 


where (y-yo)‘^H-(«-’2Q)\ 

and Jo (^) Bessers function Jq {i6). Xq, t/q, zq and a*e arbitrary 

(M. Brillouin, Compim JUndus, 1903. 


constants. 


18. Prove that a solution of the dififerential equation 
0jf2 dx^~ 

is given by 

U= J“^e‘^*[i?'(X)cos(i!Vj^) + (?(A) d\ ; 
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02 p- dV d^V 

hence obtfvin a solution of the equation ■0^+^0^°“0^ putting 

o r T 

F= U(x, t) e-‘. If ^7=/ (.ic), ^ =5' («) ** = 0) functions F, G may be 

determined by the equations 

f(_x)=j e^‘‘F{\)dk, ff(x;)=j J‘^G{\)dX, 

with the aid of Fourier’s theorem. 

(H. Poineard, Comptes Rendus, 1893-4.) 


19. Prove that with the conditions of the last example 

where ^ (•^j ^ ; 25, t) = /q V(^ - 

(Laplace (1779) ; E. Picard, Bull soc. math. t. 22 (1894), p. 2.) 


02<J> 00 

20. Prove that a solution <>=e~^F of the equation A'l»=-g^+2 


is given by the following extension of Kirchhoff^s formula : 


477 V {h, ^07 yo, «o)= J 




dV^ dV\ 
dt 07i dn ) 


+ -r1 f 

r- 07^J 


d2 




where = (^ ~ ^0)2 + (y + (2 - zq)\ 6^==(t- tif - ^ ^ 4 

and Iq {6) is the Bessel’s function Jq (id). 

The first integral extends over the volume enclosed by both the sphere 
r=ti and the surface 2; when this sphere cuts the surface the second 
integral extends over the part of the spherical surface inside 2, the last 
two integrals extend over the part of 2 which lies inside the sphere. The 
normal n is supposed to be drawn into the region of integration. If 2 is a 
closed surface and (xq, 3^0, 00) li®s outside, the region of integration is the 
space outside 2 and inside the sphere. 

(M. Brillouin, Comptes Rendus, 1903.) 


21. Let a, /3, <D be defined in terms of y, z by the equations 
y^^sscos acos(j3+^a)), y^y— (josasin(/34-^a>), y20=sin a COS 6), 
where y2= l — sinasin© and ^ is a constant. A solution of Laplace’s 
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equation A^i=0 is then given by u=yF{a,l3) provided F satisfies the 
partial differential equation 


«in2a+l-2cos2a02i?T ^ dF Z ^ ^ ^ 

■^Bo+2cos 2a ^ -F8m2a=0, 

oa^ Sin 2a djS^ da 4 


(U. Amaldi, Rend. Palermo (1902), p. 1.) 


22. Prove that the following transformations of coordinates lead to 
binary potentials, i.e. to solutions of Laplace’s equation of type'ii=i^(§, 77 ) : 

(1) y=77, z=:C; 

(2) ^=^cost, 

(3) ^=|:cosf, y^^sini', 2=r}--mC; 

(4) sin^cos?;, y = ^sin^sin77, s:==feos|; 

(5) ^=^sin77<3^^cosf, y=|sm77<?^^sinf, z=^co8rje^^ ; 


where m is an arbitrary constant. The differential equations satisfied by 
i?"in cases (3) and (5) are 


and 


0^2 + 


m2\^ 1^_ 


02i?^ 1 d^F 1 „ . 0i^ 1 dF 

(1 +^2 C0Sec2,7) ^ - (2 + m2 COSec2;7) cot 


respectively. In the other cases the differential equations are already 
familiar. 

(T. Levi-Civit^, Turin Memoirs, (2) t. 49 (1900).) 


02 7 02 7 02 y 

23. If the differential equation ^ -f == is satisfied by an ex- 
pression of type V=yf{B) where /is an arbitrary function, 6 must satisfy 
the differential equation 


d^J ^{d^J -\dt)‘ 


Prove that if we write fid$— cos adx'^- sin a dj^-r-dt, ss=t cob a +u, 
y— ifsina+v, there is a relation between a, u, v. Discuss the cases in 
which a is a function of 6 and a constant respectively, and obtain the 
general value of y in each case. 


24. If in the last example u=f{a, 6), v—g{a, 6) and we write 

0/ . dg 

^ = rsma, ^=--rCOSa, 

Oa Oa 

0/ dg 

g^=crCOS., ^=<rsm*, 
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d^V dW dW 

the new form of the equation when a, t are taken as 

independent variables is 

dr a- . . dv^ . , T/. ^ 

+ U-<r(«-r)coa(a-.)^=0. 

Prove that this equation can only possess a solution of type F=yi^(^), 
dr 

with -P arbitrary, if ^=0, and in this case 6 is defined by an equation of 

while the most general value of y is 

mL^-Si+mid)y-,,]+7t(6)[t-r-]}i 

where m, n satisfy the relation 


25. Show that wave-functions of type yf(d) may be derived from 
solutions of Laplace^s equation of this type by means of the results given 
on pp. Ill, 114. Hence show that there are wave-functions of type 
yf (^) which are not pariicular cases of a more general wave-function of 
type yf(a, /3) where a, ^ are defined by equations (280). 
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